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Abstract  This  dissertation  deals  with  a  class  of  problems  of  electro¬ 
magnetic  radiation  from  elementary  sources  in  the  presence  of  an  aniso¬ 
tropic  plasma  half-space-  All  of  the  problems  considered  have  the  following 
common  features:  (a)  the  principal  axes  of  the  Sources  are  always  parallel 
to  the  interface,  and  (b)  the  steady  magnetic  field  is  also  always 
parallel  to  the  interface  and  either  normal  to  or  in  the  direction  of  the 
principal  axis  of  the  source* 

The  work  is  divided  in  three  parts*  In  Part  1  we  are  concerned  with 
the  fields  of  a  horizontal  magnetic  dipole  in  and  out  of  the  magnetoplasma 
when  the  steady  magnetic  field  is  along  the  axis  of  the  dipole.  The  rigorous 
formulation  of  the  problem  is  carried  out  to  the  point  where  the  determination 
of  the  pertinent  boundary  coefficients  remains  to  be  a  straight-forward 
but  not  simple  algebraic  process-  Due  to  the  prohibitive  algebraic  complexity 
involved  in  the  explicit  finding  of  the  necessary  boundary  coefficients,  the 
high  frequency  approximation  is  introduced  at  that  point  and  the  approximate 
form  of  the  boundary  coefficients  are  subsequently  found.  The  field  integrals 
are  evaluated  in  the  air-region  for  the  condition  when  the  observation 
point  is  at  a  large  distance  frc®  the  source.  Several  nvmerical  examples 
pertaining  to  the  magnetic  dipole  in  the  ionosphere  are  presented  depicting 
in  graph  form  the  effects  of  the  steady  magnetic  field  on  the  radiation 


pattern  under  several  conditions  of  wave  frequency  and  the  depth  of  the 
source  hurial. 

Part  II  is  devoted  to  the  problems  of  the  electric  current  line  sources 
The  case  of  the  Steady  magnetic  field  norinal  to  the  line  source  is  rigor¬ 
ously  fomulated  and  the  field  integ^ls  are  evaluated  with  no  restriction 
on  the  wave  frequency  or  the  ma^itude  of  the  steady  magnetic  field,  fflie 
subject  of  the  lateriai  waves  is  elaborated  extensively  and  the  graphs  of 
the  radiation  field  patterns  are  presented  for  the  case  when  the  source 
is  immersed  in  the  lower  edges  of  the  ionosphere  for  frequencies  ranging 
from  VIP  to  The  latter  portion  of  Part  II  is  devoted  to  the  formulation 
and  solution  Of  the  problem  when  the  steady  ma^etic  field  is  along  the 
line  Source.  In  this  case  it  is  found  that  the  Steady  magnetic  field  has 
no  effect  on  the  radiation  field  which  could  also  be  predicted  by  physical 
reasoning. 

Finally,  Part  III  is  devoted  to  the  problems  of  the  magnetic  current 
line  sources.  The  case  of  the  steady  magnetic  field  normal  to  the  line 
source  is  rigorously  formulated  and  the  field  integral  evaluated  with  no 
restriction  on  the  wave  frequency  or  the  magnitude  of  the  steady  magnetic 
field.  This  case  is  analogous  to  the  corresponding  ease  of  Part  II  and 
many  conclusions  drawn  in  that  part  also  apply  here.  The  latter  portion 
of  Part  III  is  devoted  to  the  formulation  and  solution  of  the  case  when 
the  steady  magnetic  field  is  along  the  line  souree.  Unlike  the  corre* 
sponding  case  of  the  electric  line  source  in  Part  II,  it  is  foimd  in  this 
case  that  the  radiation  field  is  affected  by  the  action  of  the  steady 
magnetic  field  which  manifests  itself  in  the  lack  of  synmetry  of  the 


radiation  pattern. 
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CHAPTER  1 

INTRODUCTION 

This  dlssel^atidn  deals  vdth  a  class  of  problems  of  electromagnetic 
radiation  from  elementary  sources  in  the  presence  of  an  anisotropic  plasma 
half«space.  AH  of  the  problems  considered  hare  the  foiloviing  ccmnon 
featt^s:  (a)  the  principal  axes  of  the  sources  are  always  parallel  to  the 
interface*  and  (b)  the  Steady  magnetic  field  is  also  always  parallel  to  the 
interface  and  either  normal  to  or  in  the  direction  of  the  principal  axis  of 
the  source. 

1.1  OBJiCTIViiS  iil.J  THE  STATEME21T  OF  THE  PROBLIM 

We  sheill  set  forth  the  following  as  objectives  of  this  dissertation.  We 
shall  attempt  to  understand  and  explain  analytically  the  action  and  the  effeCte 
produced  by  a  steady  magnetic  field  on  the  radiation  of  electromagnetic  waves 
from  sources  in  the  presence  of  a  magnetopOLasma  half-space.  Consistent  with 
these  objectives*  we  shall  formulate  and  solve  certain  basic  boundary  value 
problems  of  eleotromagnetlc  wave  propagation  whose  isotropic  equivalents  are 
vell-ta)own.  By  comparing  our  solutions  with  the  corresponding  solutions  deal«- 
ing  with  isotropic  half-spaces*  we  shall  be  able  to  draw  conclusions  about 
effects  produced  ly  the  steady  magnetic  field. 

Having  thus  decided  on  the  basic  objectives  of  tim.s  thesis*  we  are  now 
faced  with  the  definition'  of  the  problem.  In  order  to  have  as  complete  a 
picture  as  ppsslbls  about  the  effects  of  the  steady  masnetic  field  on  the 
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various  wave  polarizations  and  the  aspects  and  directions  of  the  steady  mag¬ 
netic  field,  it  was  decided  in  the  early  stages  of  t'nis  research  to  investi¬ 
gate  the  fields  produced  by  a  horizontal  magnetic  dipole  antenna  ^th  the 
steady  magnetic  field  in  the  direction  of  the  dipole.  The  choice  of  a  hori¬ 
zontal  magnetic  dipole  as  the  source  of  the  electromagnetic  viaves  was  moti¬ 
vated  by  the  fact  that  such  an  anteniia  in  the  presence  of  an  interface  effec¬ 
tively  produces  both  vertically  and  horizontally  polarized  waves  (20),  thus, 
both  wave  polarizations  Could  be  studied  in  a  single  problem.  i4oreover, since 
the  electric  field  lines  form  circular  loops  about  the  axis  of  a  magnetic  di¬ 
pole  and  the  steady  magnetic  field  is  along  the  axis  of  the  same,  then  the 
electric  field  is  normal  to  the  steady  ma^etic  field  and,  thus,  maximum 
interaction  of  the  eleetroinagnetic  wave  with  the  Steady  magnetic  field  could 
be  expected. 

In  the  process  of  the  investigation  into  the  solution  of  a  magnetic  di¬ 
pole  problem,  it  will  soon  become  evident  that  a  rigorous  solution,  if  at  all 
possible,  would  be  extremely  complex  mathematically  and  the  physios  of  the 
problem  necessarily  obscure  and  difficult  to  Interpret.  It  will  become  ap¬ 
parent,  however,  that  certain  useful  information  could  still  be  obtained  from 
the  formulation  of  this  problem  at  the  high  frequency  limit  vdiere  the  ratio  of 
the  cyclotron  frequency  to  wave  frequency  i;;  small.  Ue  shell  explQlt  this 
possibility  a  great  deed  and  li^ed  we  shall  obtain  an  approximate  solution 
to  the  horizontal  magnetle  dipole  in  or  put  of  the  magnetpidasma  valid  over 
the  entire  hemispbere  excepting  a  smeU  region  about  a  critical  angle  corre¬ 
sponding  to  the  angle  of  total  internal  reflection  in  the  plasma  in  the  ab¬ 
sence  of  the  steady  magnetic  field. 

Since  the  investigation  of  the  vertically  and  horizontally  polarized 
waves  as  produced  by  a  horizontal  dipole  Is  not  expedient,  we  shall  Ipph  into 
the  possibility  of  study^  ^se  two  polarizatdons  separately*  Tb  this  end 
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VB  propose  to  in'/estigate  the  solutions  to  the  folloidng  tuD»dlinenslonal 
problems: 

(a)  the  electrie  ourrent  line  source  with  the  steady  magnetic  field 
normal, 

(b)  the  electric  current  line  source  with  the  steac^  magnetic  field 
parallel, 

(c)  the  magnetic  current  line  source  vdth  the  steady  magnetic  field 
normal, 

(d)  the  magnetic  current  line  source  vdth  the  steady  magnetic  field 
parallel. 

In  the  plane  nomal  to  the  axis  of  an  electric  or  a  magnetic  dipole,  the 
resulting  fields  are  identical  in  fonn,  except  for  the  space  dependence,  to 
the  corresponding  cases  of  the  line  sources.  Thus,  by  solying  the  abore  two-^ 
dimensional  problems,  we  obtain  correct  expressions  for  the  dipole  fields  in 
a  plane  normal  to  the  dipole  axis  for  the  foUovd-ng  eases  of  the  dipole^ 
steady  magnetic  field  configurations: 

(a)  steady  magnetic  field  uonnal  to  an  electric  dipole, 

(b)  steady  magnetic  field  parallel  to  an  eleotrio  dipole, 

(c)  steady  magnetic  field  nozmal  to  a  magnetic  dipole, 

(d)  steady  magnetic  field  parallel  to  a  magnetic  dipole. 

In  ^at  follows,  we  shall  have  the  opportunity  to  study  the  electromag!- 
netlc  field  structure  due  to  each  one  of  the  excitations  proposed.  In  partlo* 
ulaT)  we  will  be  able  to  elaborate  extensively  on  the  salient  features  pro* 
duoed  by  the  steady  magnetic  field  In  both  the  radiation  and  lateral  fields. 

1 .  2  HISTORICAL  BACKGROUI^  MW  RgLATED  PROBLEM 

The  problem  of  radiation  from  sources  in  the  presence  of  isotropic  interv. 
faces  is  one  of  the  oldest  aiid  nust  honored  boundary  value  problems  in  the 
history  of  electromagnetic  theory.  It  originated  with  Sosmerfeld  (20),  >dx>  in 
1909  published  the  formulation  and  the  solu'Uon  to  the  preblem  of  an 
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esolHatlng  vertical  dipole  in  the  presence  of  a  conducting,  homogeneous  and 
Isotropic  half-space.  Later,  In  1926,  Sommerfeld  fornulated  and  solved  the 
problems  of  horizontal  electric  and  magnetic  dipole  in  the  presence  of  a  con¬ 
ducting,  homogeneous  and  Isotropie  half-space.  Since  that  time  a  great  deal 
of  effort  has  been  eocpanded  by  various  workers  in  the  field,  mainiy  on  the 
subject  of  approximate  evaluation  of  the  definite  integrals  that  appeared  in 
SofiB&erfeld' a  solution.  Of  all  these  works,  the  writer  found  the  work  of 
Banos  and  Wesley  (2,3)  to  be  the  most  eomprehensive  one. 

The  phenometK>n  of  lateral  waves,  to  which  we  shall  have  the  occasion  to 
refer  a  great  deal  in  this  thesis,  was  first  reported  by  Schmidt  (19)  in  1938. 
He  performed  a  series  of  experifients  vdth  pressure  waves  produced  by  a  spark 
discharge  in  the  presence  of  an  interface  between  two  liquid  half-spaces  of 
sodium  Chloride  (sound  velocity  about  1600  meters  per  second)  and  xylol 
(sound  velocity  1175  meters  per  second.)  The  series  of  schlieren-iXiotograi^s 
obtained  by  him  clearly  exhibit  the  mechanism  by  idiioh  these  lateral  waves  are 
produced.  Subsequently,  in  1939*  Joes  and  teltov  (12)  pointed  out  that  the 
transient  phenomena  observed  by  Schmidt  can  be  readily  carried  over  to  the 
statloaaxy  fields  in  eleotrodynajuics  and  as  suoh,  those  phenomena  are  indeed 
already  included  in  the  original  Scnmerfeld’ s  solution  (20  ).  Finally,  in 
1942,  Ott  (16)  reviewed  Sommerfeld's  work  and  rigorously  showed  how  the 
oral  field  fpllovs  from  that  solution.  Later,  in  1953,  Gerjuoy  (8)  also 
elaborated  on  the  same  subject.  A  very  comprehensive  treatment  of  the  aob- 
Jept  pf  the  lateral  waves  may  be  found  in  Breldiovsklldi  book  (5  ). 

Ibe  subject  of  eleotrcwiagnetiQ  radiation  from  sources  located  in  mag- 
netopUsna  did  not,  surprisingly  enough,  attract  any  attention  for  a  long 
time.  It  was  not  until  1957  that  the  first  paper  on  this  subjeot  was  pub¬ 
lished  by  Huuldja  (6  ).  He  found  a  general  solution  to  the  problem  of  the 
field  Ffodoeed  by  *  distrlbatloo  of  external  ciurents  in  an  infinite 
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homogeneous  medium  having  arbitrary  anisotropy.  For  the  ease  of  a  magneto<^ 
active  modium»  he  found  the  multipole  expansion  of  the  radiation  field  and 
discussed  some  aspects  of  the  solution  to  the  problem  of  an  oscillating  elep^ 
tzlc  dipole.  He  found  explicit  ^pressions  for  the  dipole  radiation  field 
only  v^en  the  observation  point  was  either  on  the  dipole  axis  or  on  a  line 
normal  to  it^  while  the  dipole  itself  vas  either  along  the  direction  of  the 
Steady  magnetic  field  or  perpendicular  to  it.  The  nature  of  the  difficulties 
that  prevented  Bunkin  from  obtaining  oomplete  eocpressions  for  the  radiation 
fields  was  the  fact  that  the  saddle  point  which  is  given  by  a  solution  to  a 
transcendental  equation  could  not,  in  general,  be  evaluated  explicitly. 

The  same  problem  was  then  considered  by  other  workers.  Kogelnik  (13) 
formulated  the  electric  dipole  problem  by  a  somewhat  different  method  and 
found  the  expression  for  the  total  power  radiated  in  a  lossless  plasma  for 
the  oases  When  a  single  dipole  is  oriented  edpng  or  normal  to  the  steady  ma^ 
netie  field  and  two  crx)S3ed  dipoles  in  a  plane  normal  to  the  same.  Arbel  (1), 
Kuehl  (14) I  and  Meecham  (17)  also  considered  the  same  problem. 

The  problem  of  radiation  in  free  apaoe  from  sources  Immersed  in  a  mag¬ 
netoplasma  half-space  was  first  considered  by  Barsidcov  (4)  in  19^9*  He  ob¬ 
tained  expressions  for  the  Pelting  vector  in  the  air  in  the  case  when  the 
steady  magnetic  field  was  normal  to  the  boundary  and  an  electric  dipole  was 
either  normal  or  parallel  to  the  boundary.  He  also  presented  a  numerioal 
example  for  a  vertical  dipole  at  the  boundary  and  for  the  particular  set  of 
conditions  considered,  he  foimd  that  the  p.a3ma  anisotropy  manifested  itself 
strongly  in  the  directivity  of  the  radiation  pattern. 

Finally,  we  note  that  the  problem  of  a  vertical  dipole  imersed  in  a 
magnetoplasma  half-space  with  the  steady  magnetic  field  normal  to  the  bound- 
ary  was  also  oonsidered  by  ^bel  (1)  in  I960. 
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1.3  coN'raiauTiONs  of  this  thesis 


In  this  vrofk  we  tfeat  for  the  first  time  the  problems  of  electroaagnetiG 
radiation  from  sourees  in  the  presence  of  a  matpietoplasma  half«space  when  the 
steady  magnetic  field  is  p^allel  to  the  separation  boundafy.  We  investigate 
in  detail  the  effects  of  the  plasma  anisotro^  on  both  vertically  and  hori«* 
zontally  polarized  electroma^etic  ve./es. 

We  also,  for  the  first  time,  elaborate  on  the  subject  of  the  lateral 
tAves  in  eonnection  with  the  magneto plasma  interfaces.  In  particular,  we 
find  that  such  waves  are  induced  under  most  of  the  conditions  and  th^  fora 
an  Important  contribution  to  the  total  field  especially  in  the  proximity  of 
the  interface. 


1.4  THE  rEfiiOD  OF  ATTACK 


Unlike  in  the  previous  work  on  a  related  subject  (4)»  we  Shall  use  in 
this  thesis  the  powerful  and  rigorous  method  of  Fourier  transforms  in 
Ceurtesian  coordinates  to  formulate  the  problems.  Such  representation  neoes> 
sarily  yields  a  complete  and  unique  solution  and  furthermore  after  two  inte¬ 
grations,  allows  unequivocal  determination  of  the  path  of  integration  in  the 
complex  plane  of  the  third  transform  variable  irtilcb  must  lie  within  a  oertaln 
strip  of  analyticity  (23,  p.  44* ) 

Vie  divide  this  thesis  in  three  parts.  In  Part  1»  which  Includes  Chapters 
3»  4»  5f  6  and  7,  we  are  concerned  with  the  problem  of  a  horizontal  ma^ 
netic  dipole  in  and  out  of  the  magnetoplasma  ^eu  the  steady  magnetic  field 
is  along  the  axis  of  the  dipole*  The  rigorous  formulation  of  the  problem  is 
carried  out  to  the  point  where  the  determinatipn  of  the  pertinent  boundaxy 
ooefflQients  remains  to  be  a  straight-forwBx^if  but  not  simple,  algebraiq 
prooess*  poe  to  the  prohibitive  mlgebraio  compleKity  involved  in  the  eocplipit 
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finding  of  the  neoessaiy  boundary  coefficients,  the  high  frequency  approxi- 
aation  is  introduced  at  that  point  and  the  boundaiy  coefficients  are  subse»* 
guently  found.  The  field  integrals  are  evaluated  for  the  air-region  for  the 
condition  when  the  Observatioti  point  is  at  a  large  distance  from  the  source. 
Several  numerical  examples  pertaining  to  the  magnetic  dipole  in  the  ionosphere 
are  presented  depicting  in  graph  fora  the  effects  of  the  steady  magnetic  field 
on  the  radiation  pattern  under  several  conditions  of  the  wave  frequent  and 
the  depth  of  the  source  burial. 

Part  II  which  Consists  of  Chapters  S,  9»  and  10,  is  devoted  to  the  pro^ 
blems  of  the  electric  eurrent  line  sources.  The  case  of  the  steady  magnetic 
field  normal  to  the  line  source  is  rigorously  fomulated  in  Chapter  S  and  the 
integrals  are  evaluated  in  Chapter  9  with  no  restrictions  on  frequency.  The 
subject  of  lateral  waves  is  elaborated  on  extensively  and  the  graphs  of  the 
radiation  field  pattern  are  presented  for  the  case  when  the  source  is  immersed 
in  the  lower  edges  of  the  ionosphere  for  frequencies  ranging  from  VLF  to  HP. 
Chapter  10  is  devoted  to  the  fonmilation  and  solution  of  the  problem  when  the 
steady  magnetic  field  is  along  the  line  source.  In  this  case  it  is  found 
that  the  steady  magnetic  field  has  no  effect  on  the  radiation  field  wi'iioh 
could  also  be  predicted  by  physical  reasoiiing. 

Finally,  Part  III  consisting  of  Chapters  11,  12,  and  13,  is  devoted  to 
the  problems  of  the  magnetic  current  line  sources.  The  case  of  the  steadj'- 
magnetic  field  noirjal  to  the  line  source  is  rigorously  formulated  in  Chapter 
11  and  the  integrals  are  evaluated  in  Chapter  12  with  no  restriction  on  fre- 
quen<^.  This  case  is  analogous  to  the  corresponding  case  of  Fart  II  and  many 
conclusions  drawn  in  that  part  apply  here  also.  Chapter  13  is  devoted  to  the 
formulation  and  solution  of  the  case  vdien  the  steady  magnetic  field  is  along 
the  line  spurce.  Unlike  the  corresponding  case  of  the  electric  line  source 
in  Part  II,  it  is  found  that  here  the  radiation  field  Is  affected  by  the 
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action  of  the  steady  magnetic  field  idiich  manifests  Itself  in  the  lack  of 
symmetiy  of  the  radiation  pattern. 


RiLOlAtlON  FROM  A  HORIZONTAL  MAGNETIC  DIPOLE 
IH  THE  PRESENCE  OF  A  MAGNEIOPLASMA  HALF-SPACE 


CHAj^M  2 

MGOROUS  FOflI'nJLATION  OF  TllE  PaDBLa4  FOR  A  MAGNETIC  OlPCtt.E  SOURCE 

In  this  chapter  vie  shall  be  concerned  \^th  the  finding  of  appropriate  in¬ 
tegral  representations  for  the  Cartesian  components  of  the  field  vectors  for 
the  magneto plasma-  and  the  air-half-  spaces.  This  will  entail  the  descrip¬ 
tion  of  the  magnetoplasma's  properties  on  the  macroscopic  Scale,  review  of 
the  fundamental  field  equations  and  the  definition  of  the  source  of  electro- 
magnetlo  waves. 


2.1  STATU^EKT  OF  THE  PROBLEM 


The  geometzy  of  the  problem  is  shown  in  Figure  2.1.  The  horizontal 
plane  z  »  0  coincides  with  the  interface  between  the  anisotropio  homogeneous 
plasma  and  air.  For  convenience  we  shall  call  the  plasma  medium  (1)  and  the 
air  medium  (0).  The  plasmsi  in  addition  to  its  anisotropy^  may  be  lossy 
while  the  conduptivlty  of  the  air  is  zero.  Moreover,  we  assume  that  both 
media  have  the  same  magnetic  Inductive  oapacily  of  free  space, .  The 
factor  causing  the  anisotropy  of  the  plasma,  the  steady  magnetic  field  is 
oriented  in  the  positive  m-direction  as  well  as  the  source  of  the  electro¬ 
magnetic  waves,  the  magnetic  dipole. 

Since  we  shall  be  mainly  interested  in  the  fields  everywhere  in  the  air 
(the  upper  hemisphere),  the  final  results  will  be  eoEpressed  in  spherical  co¬ 
ordinates.  To  this  end  ve  define  the  spherical  angles  and  Q  as  shown  in 


t^c  fieldi  the  electrons  foUov  the  dii^ctlon  of  the  electric  field  between 
the  eolllelons  with  the  other  constituents  of  the  gas.  A  single  collision 
can  be  plct'ored  as  a  sudden  reduction  of  the  momentun  nu  of  a  free  electron 
to  zero.  I^ls  change  In  the  raomentun  of  the  electron  takes  place  within  an 
average  tine  Interval  being  the  collision  frequency.  ^US|  the  oolli> 

Sion  effects  are  equivalent  to  an  additional  force  «mvu  and  the  total  force 
exerted  on  the  electron  becomes 


5  +  i/  U  s  (f  e-‘“^+  U  X  )  (2.1 ) 

where  the  field  of  the  lAve  Is  assumed  to  vary  as  e^***^^ .  In  wrltli^  (2.1)  we 
assumed  that  the  magnitude  of  the  magnetostatic  field  will  always  be  much 
greater  than  the  magnitude  of  the  alternating  magnetic  field  so  that  the 
latter  one  can  be  neglected. 

If  there  IS  N  electrons  per  unit  volume »  then  the  associated  oonductlon 
current  density  Is  given  by 

e-'-fe  =  -  lei  (2.2) 

assuming  that  all  electrons  have  the  same  velocity.  (Otherwise,  we  would 
have  to  integrate  over  the  velocity  distribution.)  Combining  (2.1)  and  (2.2) 
one  obtains 


X  =  iC(je.J  p’  [-ie^la  X  +(1+lq)  (^i,)]"'  ]  I"  (2.3) 

where  ^  is  a  unit  vector  in  the  direction  of  the  steady  magnetic  field,  Si^ 


is  the  Xxonecker  delta  and 


^  S  JUtL  s. 

w 


(2,4) 


lei  Bo 


^e  total  current  in  the  plaana  7^  is  the  sum  of  the  conduction  current  X 
and  the  displacement  current  ^  »  1  oo  Thus  ve  write 
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Figure  2.1. 


2.2  FtlNQAl^IAL  EQUATIONS 


Ibe  definition  of  the  present  tvo^media  problem  Implies  the  solution  to 
Mlxuall's  equation  eubjeot  to  the  usual  boundary  Oonditions  at  the  interface 
and  proper  behavior  at  Infinity.  In  the  alrareglon  it  is  convenient  to  in« 
troduos  the  Hertsian  vector  potential  of  the  mapietlc  type  whereas  in  the 
plaa9aa~region  it  is  fotmd  to  be  more  convenient  to  work  vith  actual  field 
components.  The  imposition  of  the  requirement  of  the  radiation  condition  at 
infinity  in  addition  to  the  boundary  condition  at  the  interface  furxiishes  the 
nacesaaiy  relationships  which  render  the  solution  determinate  and  unique. 

2.2a  MsffnftnpjLAinna  descriotion^lhe  term  *  plasma"  describes  an  elec¬ 
trically  neutral  gaseous  medium  consisting  of  an  equal  number  of  free  elec^ 
trons  and  positive  ions.  We  shall  use  the  term  "magnetoplaama''  to  describe 
plasma  with  an  impressed  magnetostatic  field. 

In  the  absence  of  any  external  electric  or  magnetic  fields,  the  motion 
of  the  free  electrons  and  ions  is  random  because  it  is  caused  by  the  thennal 
agitation  alone.  Since  the  mass  of  an  ion  is  much  greater  than  that  of  the 
electron,  the  velocity  of  the  former  is  much  lover  than  that  of  the  latter. 
Thus  when  the  wave  frequency  is  high  enou^,  the  current  contribution  due  to 
the  motion  of  the  positive  ions  becomes  negligible.  In  the  presence  of  a 
steady  magnetic  field  the  electrons  describe  spiral  paths  between  oplll- 
ilons  about  the  lines  of  force.  This  motion  is  unifoxn  and  its  projection  on 
a  plane  nonnal  to  the  line  of  force  is  a  uniform  circular  motion  with  an  ang¬ 
ular  velooity  oorreepondlng  to  the  qyclotron  frequency  ,  where 

e  ^  »  1 .6  X  TO*'*  Goulcuab  is  the  electric  Phsrge  and  m  »  9*1  x  10*^'  kilogram 


is  the  alaotrpQio  smss.  IS  sow,  la  additipa,  there  is  a  tins  varying  eleo- 
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Jt  =  -  (<fi-  )  ig-%  X-(1+  iq)(5ij  )]“  }  ¥=^iw6,F2;  *  (2.5) 

We  shall  henceforth  call  (T)  the  perfaltti'/ity  tensor  of  the  ma^etoplaama. 
(24}  p.  297.)  One  eaQ}  for  con’/enience}  also  introduce  the  electric  dls>^ 
placement  sector  defined  by 

D  =  e„  €  E  i  (2.6) 

Since  the  medium  is  assumed  to  be  electrically  neutral  (there  are  as  many 
negative  charges  as  there  are  positive  charges)}  we  note  that 


^.15=0.  (2.7) 

The  components  of  the  permittivity  tensor  which  we  shall  need  later  have  been 
found  elsevdiere  (24}  p.  297)  for  any  orientation  of  the  steady  ma^etic  field 
.  When  the  steady  magnetic  field  is  oriented  in  the  positive  x-direetlon, 
idiieh  is  the  case  under  consideration,  the  tezrsor  (?*)  is 


(2.8) 


Since  'T  is  non-singular}  its  In/erse  can  be  found  to  ha/e  the  form 


(2.9) 


%bere 


and 


K  S  jm7 


<  =  5'+  i?" 

6  =  €'  +  i€  " 

7  -  7'+  I7" 


(2.10) 


(2,11) 


Moreover 
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1 


1  +  q* 


r= 


+  q< 


4.  s  .  pMi-<^\-^a 

l(«r  +  -0*  +  qi]{  y  Y  ^  q*] 

*  B  qp^(H-g^-*-  q^)  ^ 

«'  a  _ _ <yp^  Ji^  ) 

”  *  [(er+  1)»  +  q2j[  (ff  -  1)*  +  q»] 

'  [(^+  D*  +  q^lt  (ff  »  1)*  +  q*]  . 


>  (2.12b) 


\  (2J2e) 


When  the  oolllslona  ax^  emali  and  can  be  negleoted»  the  above  parameters  have 


simpler  foms  as  folloust 


$  a  1  -  p^ 


7 


1  ^9*  • 


(2.13a) 

(2.13b) 

(2.130) 


FiurthexBore,  idien  the  magnitude  of  the  magnetoBtatio  field  is  amall  or  the 
wvs  Arsqaeno^  is  suffleiently  high  so  'Uiat 


ws  oah  lirite  for  a  lossless  magnetoplaaDa 

^  a  1  ^  p»  (2.15a) 

(2.15b) 

*  ©•  p*  ,  (2.15c) 


As  a  ppasequanoef  the  rations  (2.10)  pas  be  wittan 
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1 


&  P* 
* 


(2.16) 


Ue  note  that  for  k  to  remain  small  p^<^  1  or,  in  other  vords,  the  wave  freo 
quenpy  imist  be  much  larger  than  the  plasma  frequency. 

Vie  shall  Introduce  still  a  different  notation  vdiich  vdll  prove  useful 
later.  Namely,  we  define  an  index  of  refraction  of  the  plassaa  in  absence  of 


the  ma^etostatic  field 

n^  =  1  -  p^  ,  (2*17) 

in  terms  of  which  equations  (2.15)  and  (2.16)  become 

<  =  n^  (2.18a) 

e  =  n^  (2.18b) 

7  -  n*)  (2.13c) 

X^n^  (2.18d) 

K—  -  &r  (JL|JJ!^)  (2.18e) 


2.2b  The  nature  of  the  3ource~-ji*or  the  purpose  of  the  present  problem 
it  will  be  assiimed  that  the  source  of  the  electromagnetic  waves  consists  of  a 
small  wire  loop  carrying  an  alternating  current  Such  a  small  loop  can 


be  fed  with  a  standard  coaxial  cable  trananlssion  line  with  a  proper  matching 
section  in-between.  If  the  loop  is  smadl  enough  its  electromagnetic  effects 
can  be  adequately  represented  by  its  equivalent  magnetic  dipole  s^ment 
(7,  p.  291) 

!■=  is'e'‘‘^^  (2,19) 

>d)ere  ^  is  the  surface  vector  of  area  enclosed  by  the  loop.  The  magnetic  di¬ 
pole  moment  can,  in  turn,  be  thought  of  as  a  product  of  ’^magnetic  current"  K 
sgid  the  length  of  the  dipole  "1",  thus 

K  1  •  t.,  -  i"  .  (2.20) 

Ip  localize  the  source  properly  we  write  for  the  "magnetic  current  density" 


(2.21) 
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1^-  I  (x)  ^  (y)  S  (z+h) 
i^ere  S  is  the  Dirac  delta  function  (7,  p.  43) > 

2.2.e,  jrh&_field_:eaiuatlQns-*We  ha/e  sho\m  that  our  soturce  of  electro- 
inagnetlc  waves  may  be  regarded  as  a  magnetic  dipole  singularity  located  at  a 
point  (0,0>-h).  Then  for  the  magnetoplaSma  region  the  appropriate  form  of 
Maxwell's  equations  is 

0  X  f ,  =  -  a  iT (x)  ^(y)  S(z‘^h)  1^ 

^  X  H,  =  £, 

(2»22) 

^  *  d]  -  0 

V* I  = 


where  is  the  fictitious  "magnetic  charge  density"  related  to  the  "magnetic 
current"  by  the  continuity  equation 


j  =  « 

"  7>  t 


(2.23) 


Since?  is  non-singular  the  second  equation  of  (2.22)  can  be  written 


T'vxh[  =-ico€oT.  .  (2.24) 

Now  we  operate  with  the  curl  on  the  last  equation  and  substitute  the  first 
equation  of  (2.22)  to  obtain 

-  V  X  ^  X  H,  +  H,  -  -  i<y6.m^(x)  ^(y)^(z  +  h)  .  (2.25) 

Using  the  inverse  of  the  permittivity  tensor  of  equation  (2.9)  and  carrying 
out  the  necessary  algebraic  operations,  one  obtains  a  set  of  simultaneous 
equations  as  follows: 


*>1 

o 

*■ 

w 

'  1 

0 

(2.26) 


19 


l^ch  we  slmll  leave  in  this  foftn  for  the  time  beih^;.  From  the  second  equa¬ 
tion  of  (2.22)  we  note  that  the  electrio  field  can  be  expressed  in  terms  of 
the  magnetic  field  as  follows: 


d 

k  a  1 

W./ 

c 

_  __  j  _ 

yV  A 

H,. 

fcyi 

(2,29) 

dy^LK^Z 

-a* 

tK<9« 

Mi, 

In  the  air»region  it  vail  be  convenient  to  use  the  Kertaian  vector  po¬ 
tential  of  the  magnetic  type  defined  by"** 

^  xTa  (2. 2d) 

and  satis^ing  the  honogeneous  vector  wave  equation 

(vH  k^jfFo  =  0.  (2.29) 

the  magnetic  field  is  expressed  in  tenos  of  this  Hertzian  vector  as  foll,om: 

B,  *  k*f.  +  ^.  (  ^  )  (2,30) 


2.3  FOURIER  INTEGRAL  REFRE^jC^TATIOi^  IN  CARTESIA.N  COORJlIUTES 


The  formulation  of  the  present  boundary  value  problem  can  be  simplified 
a  great  deal  by  expressing  the  field  components  in  the  magnetpplasDraa  the 
components  of  the  Hertzian  vector  in  the  air  in  terms  of  their  triple  Fourier 
integral  representation  in  Cartesian  coordinates  in  the  transfonn  space  as 
well  as  in  the  configuration  space.  Such  representation  necessarily  yields  a 
complete  and  unique  solution,  and  furthermore,  after  performing  two  integra¬ 
tions,  it  allows  the  vnequlvocal  determination  of  the  path  of  integration  in 

In  the  literatiure  on  electrCTnagnetic  theory  the  Hertrian  vector  of  the  elec¬ 
tric  type  is  denoted  by  if  and  the  one  of  the  magnetic  type  by  if  •{22,p.29), 

In  ^8  thesis  we  use  only  the  Hertzian  vector  of  the  magnetic  type  and  we 
shall  denpte  it  by  T  reserving  (*}  to  denote  the  comi^ex  conjugate. 
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the  domplex  plane  of  the  third  transform  Variable  uhich  mist  He  vLthin  a  eel^ 
tain  strip  of  analytlclty  (23»  p*  44)* 

ThuS)  ve  Introdace  a  triple  Fourier  transform  pair  defined  by 

and 


(2.32) 


In  idiat  follows  ue  Shall  also  need  the  transforms  of  the  derivatives.  These 
can  be  obtained  Integrating  by  parts.  ThuSf  for  eKsaple,  oonslder  the 
integral 


4* 


■f  c^UjiF  etKclycIz. 


(2.33) 


and 


«• 


(2.34) 


JJ 


The  vanishing  of  the  integrated  part  of  equation  (2.33)  and  (2.34)  at  the 
upper  and  lower  lirilts  is  assured  for  all  real  values  of  the  transform  yarl* 

able  pt, ,  providing  we  obey  the  radiation  condltloni  l.e.,  with  y  and  s  fixed 
wo  should  have  as  a  solution 
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F 


iLklJtl 

757 


(2.35) 


for  Im  0  . 

IhuSf  assuming  that  our  solution  satisfies  the  radiation  sondltion,  ve  es** 
tahlish  the  follovd^g  eorresfiondences: 


3^f  <=^  i«.f 

l«,f  (2.36) 

F  i«j? 

2. 3a.  Kie  parti oular  integral  corresponding  to  the  sourae—To  transform 
the  Inhomogeneous  ^Stem  of  simultaneous  equations  (2.26),  one  multiplies 
both  sides  by  the  cdcp  [  and  Integrates  vdth  respect  to  tiie 

real  Variables  x,  y,  and  i  between  •>  and  +  .  The  right-hand  side  of 

(2.26)  yields  at  once 


•‘♦«»y***)dkdydx  «  e 


.'ojK. 


(2.37) 


and  the  left-hand  side  transforas  according  to  aq|uations  (2.31)  and  (2.36). 


Thoa,  one  gets 


x: 

gca.k 

0 

0 

The  above  system  of  algebraio  equations  can  be  solved  using  Cramer's  rule. 

The  detersdnant  of  the  coeficients  of  the  square  matrix  on  the  left-hand  side 
can  be  found  to  be 
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A  =  -  [(X  i) k\  -  zott]oc;- 


(2*39) 


We  feeognize  the  above  expression  as  a  fourtn  order  equation  in  «|  .  It  can 
be  readily  put  in  tlie  fora 

^  (<V3^  s.^  )  («/-  s*  )  (2.40) 

\diere  and  a^  are  ^i/sn  fc/* 

}  (2*41) 

Cio  result  of  equation  (2.41)  is  not  Surprising.  It  is  typical  of  what  could 
be  expected  from  analysis  of  a  birefringent  mediam  like  magnetoplasma*  We 
found  two  charaGteristic  modes  of  propagation,  s,  and  s^  ,  with  vdiich  terms 
"ordinary"  and  "extraordinary"  waves  are  often  associated.  We  note,  moreover, 
that  the  determinant  of  equation  (2.39)  contains  all  possible  Information 
about  the  inodes  of  propagation  within  the  magnetoplasma.  These  inodes  are 
given  by  the  zeroes  of  A  ,  vdiich  are  actually  the  poles  of  the  field  compo¬ 
nents.  We  note,  for  Instance,  that  if  o<3  =  0  in  equation  (2.39),  vMch 
is  equivalent  to  saying  that  the  wave  does  vary  along  the  y  and  z  eoprdinatas, 
we  obtain  for  the  charactc  ristic  modes 

(«•  ).,z  -  (2.42) 

idiiph  can  be  immediately  recogt^^zed  as  the  plane  wave  propagation  constants 
in  the  direction  of  the  steady  magnetic  field  (24,  p.  297)*  Alternately,  if 
we  set  Qt,~  =  0  and  solve  for  o<^  we  obtain 
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(«*).  = 

(«i)i  =  k,/r 


(2.43) 


\Moh  Gkgaln  dan  be  recognized  as  plane  uave  propagation  constants  transverse 
to  the  steady  magnetic  field  (24»  p*  297.)  Setting  <x.=  0  and  solving 
for  d()  gives  the  same  results  as  in  (2.43)  as  it  should. 

We  now  proceed  vdth  formally  solving  the  Eystem  of  algebraic  equations 
(2.3S).  Using  Cramer's  rule  and  the  results  of  (2.40)  aiid  (2.41)*  ve  obtain 


Luie.mi 

- 


Sr»a 

(2if^x  A 


(2.44) 


idiere 


(2.45) 


^ere  nov  ve  assume  to  be  compleqic.  The  integrand  has  two  pairs  of  poles 
In  the  pompleac  ot^wplaiie  which  can  be  either  real  or  owplex.  If  they  are 
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I 


flgiir*  2*3  »  Xh» 
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(2.43) 


dosiiAaK  vm  assume  that 

I  m  {  s,  ^  0 

I  m  {  Si  J  0  * 

Then  for  s  -t-  h  >  0  we  close  the  contour  In  the  upper  half  plane  as  shown  in 
Figure  2.4  and  write 


-  It," 

e.c, 


Res. 


-R  e, 

On  C,  we  set  o<jS  Re*’^ «  then  for  R  we  obtain  for  the  second 

on  the  right 

f _ _  i  f^^Rt'*).^^**^"”***-^*-^**^*  r-f  _ _ O 

J  Fi— *.4  J 


(2.49) 


(2.50) 


slnoe  the  order  of  the  numerator  is  nerer  larger  than  three  (see  (2.44)  nod 
(2«45).)  Thus 


(2.51) 


Nov 


and  the  integral  in  (2.51)  beTOoes 


ri?5 


(2.53) 


For  s<f  b<0  we  close  the  contour  in  the  lover  half  plane  and  by  an  analogous 
argument  obtain 
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6.Cs»-s‘)  s^CS.'-sl)  j 


(2*54) 


then  the  transforms  of  the  magnetic  field  components  vihen  in'/erted  vdth  re» 
speet  to  the  oij  transform  '/ariahle,  become 


«  <x*.  s.  )€  ^*‘**'*^*  ^  #.  a*  )!g N 

“♦TxwSl  e,cs.^4i)  s»a.*“a/) 


} 


(2.55b) 


^  0ii.^x(«,j «t.  S* J  T-t K-X^klst  V 

^  J 


(2* 55b) 


**“xo  4TX>4  I  S,(S»^st) 


SafS.^S/)  ^  ]• 


(2.55c) 


Finally,  ye  invert  with  respect  to  ot,  and  obtain  the  desired  represen¬ 

tations 

•• 

§f^\lrsV(s.*-«J)  %c«f-s?)  J  ■  'dm.dcfz 


uM  .  CW€.>'»X  3,  rr fo<>^«C»)  yfcKXifcSs,  gCz.la^Kl 

**K^o  ef*^  ^  JJl  s.<s.*= 


(2*56b) 
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L-jW  ^  tujjs»i<yv_a  ritS.  #kCs,)^tK  Xr!*l«i 


j  t ,^.lS2££.!£,»ilS5  g  fc'Sfci «»k I  le<-'(«,iii  ♦or^wj  j  j 


(2.56c) 


2.3fa  The  GQffljjletnc.ii'tar-/  field„in  J.h.o„plasma^ln  tlie  pre'/ious  section  we 
obtained  the  " particular  integrals'*  of  the  systen  of  equations  (2. >3)  which 
represent  the  primary  excitation  due  to  the  source.  To  satisfy  the  boundary 
conditions  of  the  problem>  v/e  shall  need  an  appropriate  Conplemontdry  Solu* 
tion  of  the  homogeneous  system  of  (2.33).  jj*roin  the  theory  of  linear  differ* 
ential  equations  (10|  p.  145)  it  follows  that  the  magnetic  field  components 
satisfy  the  equation 


(2.57) 


\4iere  oc,  in  this  equation  is  considered  aS  the,  difi’erentiol  operator  with  ro* 
spect  to  the  z-  coordinate.  The  solutions  to  (2.57)  cr.n  bo  written  down 
immediately 


m  • 

— 

i 

C., 

e.^ 

_  _ 

Q-r^S.Z. 

H^: 

= 

til 

C3. 

(2.58) 


where  we  discarded  solutions  with  positive  exponentials  since  we  can  have 


only  waves  going  away  from  the  interface  upon  reflection.  The  coefficients 
C^j  are  not  all  independent.  For  since  ^  has  to  satisfy  the  given  system^ 
the  coefficients  must  be  related  by 


F„  (s-s, )  F,i  (*s, )  F,3  (»s,  ) 

'c,.’ 

F^,  (^s, )  Fjj  (’3, )  Fii  («s,  ) 

Cz, 

(-S. )  Fji  (••s.  )  Fjj  (..s,  ) 

(2.59) 


(2.60) 
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p 

*11 

(-Sa) 

tz  (-32) 

P 

‘13 

(-Sa) 

F 

'^21 

(*Sa) 

Faa  (^Sa) 

^23 

(•“•Sa) 

F 

*3l 

(“‘Sj) 

F32  (-"32) 

?33 

(-S2) 

C.. 

C« 

C32 


vtfiere  the  square  matrix  elements  (s.^j)  correspond  to  the  matrix  of  the 
original  system  (2.33)  with  o(j=  -s,  or  oc^as  -s^  ,  Now  ve  express  wyi  00* 
efficients  C^j  in  terms  of  C„  and  C,2  *  Ihe  j^oedwe  is  Stralght<*forward  and 
w  obtain  as  the  results 


Cii 


• 


and 


^22 

Cji-  *  C.a 


(2.61) 


(2.62) 


For  convenience  in  what  follows  wa  shall  nomallae  the  coefficients  C,,  and  C,,^ 
as  follows t 


C 


II 


C4J  6,0 

ATT  X  Ini 


k, 


Wxk^ 


^2 


(2.63) 


Now,  inverting  with  respect  to  pi,  and  0(^,  we  obtain  the  representation  of  the 
eomplementaxy  field  in  the  piiiiiiia 


C2.64b) 
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”*'  8Tr‘Xkl  JJ I  #.Vs.)  ■  ® 


i|a.<3ld(2 


C^64c) 


As  a  partial  checic  on  the  above  results  we  Investigate  \diether  »  o  as 


it  should.  Indeed,  we  find  that  the  abovre  is  true,  providing 

#.(3,,2)  -  (<^x+  sf,^)  #i(s,,j)  (2*65) 

Which  follow  frora  (2.45)  and  (2*41).  In  this  connection  we  note  the  follow* 
ing  useful  and  not  too  obvious  relationships 

X  (sf  +  s^  )  =  Id  (xko  oi)  -<(1-K*)(v*»Xa?  (2.66a) 

‘Xajt  *  4*i(3|,z  )  +  X  i  %'<l  -  a**  0^)  (2.66b) 

X(sf*s|'  )*  #2(22)*  #4(3.)  (2.66c) 

^z(ai,t  )  ®  (‘^*+  ®.f*  *  Xko  )  (2.66d) 


2.3o  The  field  in  the  ali^As  we  noted  in  section  (2.2c)  the  fields  in 
the  air  are  derivable  from  the  Hertzian  vector  potential  of  the  magnetic  type 
satisfying  the  homogeneous  vector  wave  equation 

(v^+ k|)Tfp  s  0  . 


The  foxin  of  solution  inonediately  suggests  itself  from  examination  of  the  c^oh 
pleoentary  field  in  the  plasma.  Ue  write 


'MO 


% 


20’ 


8ir 

C06, 

8Tr 


ff  B,  e ' d*. doTj 

‘Xki  JJ 

%• 

JJ  B*  d«.d«2 


(2.67a) 


where 


(2.67b) 


30 


So  =  yj  -  af-  Oil 
Im  ^3o}  0 


(a.6a) 


Here  we  assuaed  that  the  vector  TT  has  onlj’-  two  non-vaJiishing  coropOhents 
WhLch  are  sofficieht  to  derive  all  of  the  electric  and  magnetic  field  Compo^ 
hehts  and  to  satisfy  the  boimdary  conditions. 


2.4  THE  BOlfHDARy  CONllTlONS 


In  the  preceding  sections  we  have  found  field  comixjnents  in  the  plasma 
and  aira regions  that  are  solutions  to  Maxwell's  equations.  Moreover,  in  sol¬ 
ving  Maxwell's  equations,  we  have  Chosen  such  solutions  for  field  representa¬ 
tion  that  have  proper  behavior  at  infinity  by  requirihi:  that  the  imaginary 
part  of  the  pertinent  exponents  be  non«negative.  In  addition,  the  field  Com¬ 
ponents  contain  certain,  thus  far,  undetermined  coefficients  yhldh  upon  im¬ 
position  of  the  boundary  will  be  determined  and  thus  render  the  solution 
unique. 

The  boundaiy  conditions  to  be  satisfied  by  the  Cartesian  components  of 
the  field  vectors  require  continuity  of  the  tangential  components  of  the  eleo- 
trio  and  magnetic  fields  at  the  interface  s  -  0.  This  implies  the  foUowingi 


Hxo  -  H,,. 

Hyp  -  Hy. 

Hxo  “ 

Hyo  ®  ®yi 


(2.69) 


Ue  rewrite  the  above  equations  using  the  Hertzian  vector  representation  for 
the  fields  in  the  air 


I^TTco  *  ^Z%a)  *  W^. 

^^o^y^Zo  ~  ^y^Zi  ~^z^yi 


(2.70) 


The  abdv«  system  of  equations  can  be  somevAlat  simplified.  Oiie  can  show  that 


r,i^  ,  Q<._.  y 

"xb  ^‘Xi  ^yi 

Klho  3^  Q(^Q^  Hyi 

k^^  =  -L  ilia.  +  X  H 


(2.71) 


k*ir  -  ^  r2>  +  Mi/J*  J.  -  <\u  +  ^? 

*o  *'xo  ~  ^  „  (T«  "  v7”yi^  if  V  s  “*i 


■“  as  X'  >• 


yhefe  we  have  used  s/  to  denote  the  field  e^cpfessions  vdth  the  integral  signs 

V  V 

removed.  Elimination  of  and  71^^,  gives 


0  *  o(z  (a,-  iX  s^-  *caQK,  -  cz,(^  4  -  iXSo  -  i  O'  «2(^  1 )  Hj, 

0  =  io<.o<*i;H^,  +  [-  Sb^z  +  (iCo  -«^)  ]  H^,  +  i«zsX. 


(2.72) 


Upon  performing  the  necessary  algebraic  operations,  one  obtains  tuo  simal<* 
taneous  equations  in  tuo  iudtnouns»  A,  and  .  We  write  the  results  in  a 
Symbolle  form 


(2.73) 


where  the  various  coefficients  are  given  by 
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(fSj^cKOC^kfs. 


^yv«>v^/V  ^  ^\i-S,  4>,.CS.)  *  UXikjcC^'] 


(2.74a) 


b.  afi^iCs.-  XSo-t^cKoti)  #.a.) 

-C«t(s**Xs.^i  Ko<*)'i>,Cs.)  -*-a.»(^  s^-Xs^  -"tK  X^I<Js.) 


«.*«,(t‘  'X«»4(Si)+*:KXf  k,V)l 


(2.74b) 


(2.74c) 


-1 1.^  ~  ><:KXk?Ks.cs.WJ-»-5.<:f.^.H0b^)1 

~  - - 


(2.75a) 


#.csO  " 


(2.75b) 


^  =  -{«*C#.CS,}-Cs#*0«^)#»Cs.>]+CKX  -gf-f*.*^ 

Jsc<f,*^sf) 

Once  the  ooefricienta  A,  and  A^  are  found,  the  tup  remaining  coeffipienta  B, 
and  Bj  can  be  detennined  from  equationa  (2.71.)  Thua,  ve  find 
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a*B,  =  { «a [ #1  C».)  ^  XjTkSs.a*  j 

+{6(i[#,CsO‘^«.*‘i*i<fs*)]  ?  kaSaa.*} 


and 


2£»^Bi={az  [#,  fs.)  -  X  ? 


_A, 

t'fs.CS.) 

CsJ-*at)j 

s.(rs»-sz‘; 

(2.76) 


(2.77) 


The  explicit  detsimnation  of  the  four  bouadarj;  coefficients  A,,  A^,  B,,  and 
Bj  is  nov  a  straight-forward  although  a  tedious  matter.  The  results  would 
necessarily  be  lengthy  and  probably  not  very  useful.  We,  therefore,  forego 
their  evaluation  in  this  form,  in  the  next  chapter,  we  shall  Introduce  a 
high  frequent  approximation  vhich  will  simplify  the  results  a  great  deal 
without  the  loss  of  the  basic  ingredients  of  the  problem. 

2.5  O.OSURE 

In  the  foregoing  chapter  we  have  rigorously  formulated  the  problem  of  a 
horizontal  magnetic  dipole  in  a  magnetoplasma  with  a  separation  boundary.  We 
oarrled  the  analysis  up  to  the  point  where  the  determination  of  the  pertinent 
boundtny  coefficients  remained  to  be  a  straight-forward  (but  not  simple) 
algebraic  process.  Due  to  the  £p,gebralc  complexity  involved  in  explicit 
finding  of  the  desired  boundary  cpefflcients  It  was  cpnoluded  that  their 
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final  fonns,  if  they  could  be  found,  would  be  of  limited  ralue  in  interpreta* 
tion  of  their  physical  significance.  Therefore,  their  explicit  evaluation 
was  postponed  until  the  next  chapter  in  vMch  ve  anticipate  the  introduction 
of  an  approximation  that  would  make  the  results  algebraically  manageable  and 
their  physical  interpretation  feasible. 


CHAPm  3 

HIGH  FR^UEllCX  APPROXBt^TIOll  FOR  DIPOLE  IN  MAGNEtOPLASI'A 

In  the  preceding  chapter  ve  have  fonBilated  the  present  boundary  value 
problea  finding  the  integral  eKpressions  for  the  fields  In  the  plasma  and 
alx^regions.  these  integral  expressions  contain  certain  boundary  coefficients 
vhieh  are  expressible  in  terns  of  the  media  properties,  source  strength  and 
frequency.  Even  though  in  theory  these  coefficients  are  known,  their  useful¬ 
ness  in  reducing  the  field  integrals  to  form  suitable  for  numerical  calcula¬ 
tions  la  of  Very  limited  value  due  to  their  complexity. 

In  this  chapter  we  shall  be  concerned  with  fiixding  appropriate  approxima¬ 
tions  that  would  put  the  previously  derived  results  in  manageable  form  which 
could  be  Intexpreted  piysically. 


3.1  THE  NATURE  OF  THE  APPROXIMATION 


Aiy  approximation  that  we  may  Introduce  to  simplify  the  field  components 
is  necessarily  contingent  upon  simplifi  cation  of  the  propagation  factors 
corresponding  to  the  a-dlrection,  i.e*,  s,  and  which  are  given  by  equation 
(E.41)*  ^ba^  the  most  useful  approximation  will  be  one  mentioned  briefly 
in  Section  and  based  on  the  assumption 

I  I  ^  1  (3.1) 

I  m  I 

wbioh  can  be  brought  about  by  either  a  weak  steady  magnetic  field,  ,  or  by 
a  sufficiently  hi^  wave  frequency,  to  .  As  a  consequence,  vdien  the  losses  in 
the  are  nai^gible,  the  idasaa  parameters  can  be  written  (sea  Section 
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6  ~  ft* 

^  -  ef  (1-n*) 

X  -- 

i&ere  *a*  la  the  Indax  of  reff^otlon  of  the  plasma  la  absenoe  of  the  steady 
magaetio  fieldi  as  ^  folloin  ue  shall  require  k  to  be  flitall* 

This  tdll  thea  exolude  oases  la  \Ailoh  "a^  Is  close  to  aero. 

We  aote  that)  oonslsteat  uLtii  the  apprcacloatloas  of  equation  (3>^}>  the 
permittivity  teaser  (2.8)  and  Its  Inverse  (2.9)  beoome  to  the  first  order  InK 


approximations  of  equatloas 


(3«2)  are  nov  Introduced  Into  the  equation  (2.41)  giving  apprucimate  foxms  for 


i,  and  as  foUoest 


k,*ak,  • 

If  ofldy  first  order  terms  in  k  are  kept*  the  above 

9*t  •  ±  Kk,«,  • 


(3.5a) 


(3.5b) 

(3.5o) 

aioplifies  to 
(3.6) 
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Pinally,  In  regions  vhero  |s^|  »  |  k.]c,6<,|  one  can  further  approKinate  tgr 


uniting 


a 


111. 


_  4  Kk.Ot, 

®  - 


(3.7) 


In  uhat  foUove  ue  shall  ooncl  uie  that  the  last  approximation  Is  Indeed 
the  only  useful  one  uhen  oaleulatlng  the  amplitudes  of  the  field  oomponents  In 
the  air.  The  usefulness  of  It  ls»  however,  reduced  by  the  fact  that  It  Is  not 
vaUd  In  the  vlolnlty  of  s  »  0  which  oorresponds  to  the  angle  of  the  total 
internal  reflection  in  the  plaaoa. 

3.1  b  aporoxlaflte,  forinff  and  «.-ainr»a  both  <§,  and  4,,oontaln  only 

squares  of  a,  and  s^  >  their  approximate  form  need  not  be  restricted  by  the 
requirement  that  |8*|  >>  |tck,o(.|  .  Thus,  Using  (3.6)  only,  we  rewite 

(2.65)  and  (2.66) 

<it(ai,i)  =  ?  n*KK,«,  (3.8) 

and 

‘i».(s,,z)  *  7  n*Kk,6c.  (s*+a^i  Kk,o<,  )  .  (3.9) 

Ue  also  note  for  future  reference 

*  4  h*K  kfoc.  (k,  !  K«,)  .  (3.10) 

3.2  iUPPRQXmTE  FORKS  OF  THE  BOUNDARI  COEFFICIENTS 

In  the  preoe  ding  section  we  have  derived  oertaln  approximiatlona  that  are 
valid  when  the  steady  magnstle  field  Is  small  enough  or  the  wave  frequent  Is 
high  enough.  AppHoatlon  of  this  approximation  to  the  expression  for  the 
propagation  factors  s,  and  s^  resulted  In  considerable  slmpUfleation.  Sub., 
saquently  we  fouhd  It  neoeesary  to  Introduoe  another  slmpUfleation  iddcdi 
resulted  In  limiting  the  range  of  validity  of  the  apprraimation  for  the  regloa 
where  Is^  I  »  Ikjc.oi,  |  . 

In  this  seotlon  we  shall  apiiy  thla  approximation  to  obtain  simplified 
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expressions  for  the  field  coaponents  in  the  air. 

We  apply  equations  (3*6}  and  (3.8)  through  (3*10)  to  the:  matrix  coeffl* 
cleats  of  equations  (2.74)  and  (2.75)  to  obtain 


a„s  e  -  *cd 
a, 2=  0*  +  Kd* 
a22*  e  *  Kf 
a4,s  6“  +  Kf* 

vhere 


(3.11) 

(3.12a) 

(3.12b) 

(3.120) 

(3.12d) 


Moreover, 

bi  B  (3.13a) 

bj  * ’^{(ixM«zsVi^[swC5Wx)Whs^ri^s(6l<..pr^) *s.(s*..«/))]je‘*H  (3.13b) 

We  alio  revrlte  equations  (2«76)  and  (2.77)  in  the  atm  approiclJBate  font 

*Ka.*Ja.or2(A.*A*)^^(s^-<)^(A,  +Ai)]  J 


(3.14) 
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B*  »  *Ai)  » Ck,s(s^-^^iX^^*^i ) (A. *  Az ) 

♦Kk,ot,fs*-»>a{)f*.a4(A.“A»)-^(^6*^«^K^^  ‘►Aa)]  } 


Next  \fd  eif^uate  (A,  ^  k^)  and  (A,  •  A^  )  \itier6  A,  afid  A^  are  giren 
equation  (2.73) •  ^or  the  determinant  of  the  square  matrix  we  find 


(3.15) 


(3.16) 


7  ^ 


Hovy  using  GrainerU  rule  we  find  from  equation  (2.73) 
*  «  B|  *  Six  62 


«  A 


toploylng  the  notation  of  equation  (3.11)  one  can  show  that 
(A,+  Aj)4=  2  [^iln(e)  -  KRe(f)3  6,+  2  [lim(c)  -KRe(d)] 


(3.17) 


(3*18) 


(3.19a) 


(A,-  Ai)A=  2  [Be(e)«iK3a(f|b.«2 


•  Ik 


^ch  amounts  to 


(3.19b) 


(3.20a) 


ss.Ca'^^fcX****  «•*)]  j- 
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CUiU 

(A,-AJA  ^  j'4C>a*-0ar,s» 

[•" ^Kt«*0^C^ s  SoC^s Vafe*j)  (3,20b) 

^(sVcy»)^('ix*^i),^,S*^k‘o^*)+ss.CriVi)('s^«»*;j^  ^ . 

Substituting  the  above  results  into  equations  (3*14)  and  (3*15)  and  also  using 
equations  (3<16)  and  (3* 17)  one  obtains  finally 


®t  •  •  « It.  {(;»f 

_ ^ _ IIU"'-. 

V.a*.S<>-K.)  «»sCs»'a»s,)/J  J  ^ 


We  note  that  in  the  above  boundaxy  coefficients  we  have  clearly  separated  the 
components  that  are  dependent  on  the  steady  magnetic  field  (pireceded  by  k  ) 
from  those  that  are  not.  I<k>reover,  we  note  that  the  components  of  B,  and 
which  are  Independent  of  the  steady  magnetic  field  are  the  proper  boundary 
coefficients  for  the  analogous  problsn  of  the  magnetic  dipole  in  an  isotropic 
half-space  (2|  p,  26-28).'*' 

The  components  of  B,  and  produced  by  the  steady  magnetic  field  display 
some  heretofore  not  encountered  denominators  >dilch,  as  we  shall  see  lateTi  are 
responsible  for  the  unusual  effects  produced  by  the  plasma*  s  snisotrppf. 


The  reference  quoted  contalzis  soluUon  to  an  electric  dipole.  The  ecBperl- 
son  with  our  resists  can  still  be  made  by  ohanglng  the  appropriate  tnns- 
idssien  ooeffioiento  to  the  magnotio  typo. 
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3.3  THE  COITONEIIIS  OP  THE  HERtZlAN  VECTOR  IN  THE  AIR 


In  the  previoua  Motions  ve  found  the  approadmate  forms  of  the  boundary 
ooefflolents  that  deternine  the  amplitudes  of  the  eomponents  of  the  Hertslan 
▼eetor  In  the  air.  The  pertinent  eocpresslons  are  neoessarily  long  even  in 
their  simplest  form.  In  this  seetlon  we  shall  endeavor  to  arrange  the  results 
for  the  airareglon  in  a  suitalde  form  for  approjElaate  evaluation  of  the 
nsoessa^  integrals^ 

3.3a  The  Hertslan  vector  in  Cartesian  ooordlnateS‘*^We  substitute  the 
results  of  equation  (3.21)  into  the  fonmilation  of  the  Hertzian  vector  in 
equation  (2.67)  and  obtain 


and 


_ _ 

«»Ks*s,Xa-^n^v) 


Vslac  ^  identities 


$♦0. 


(3<£3») 
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azid 


oc,  4= 


S6 


^  -  i^z 


s  4,3^, 

we  c&fl  recast  the  above  inteiJfals  In  a  siraplef  fofin  as  foliows: 


TLa- 


*cu-e^a.*Ha  i  / 1 _ 2 

iTTlut 


_ _  ..  Z -  ^  ^  C  ^  ^  J _ J  ,. 

8Tr»k5  ^ 

-  *o 

C«9 

«« 

C9 


and 


*n^Q= 


(3.23b) 


t.eldtz 


(3.a4a) 


€09 


S<Si-  ri*s4 


g  CC*. K * or^y  >  S.X +SK)^^ 


•!-C0 

,^^(4.i,)|r-^e*<««— .ys.*»K)rf,,st«, 

c<» 


(3.24l>) 
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In  the  above  expressions  vre  have  alenrly  separated  the  contribations  of 
the  steady  magnetic  field  in  the  plasma  by  enelosing  them  in  square  brackets, 
the  first  tern  in  each  of  the  above  Hertzian  vector  components  can  be  recog^ 


nized  as  the  proper  expression  for  the  problem  of  a  magnetic  dipole  in  a  dl- 
electric  with  an  index  of  refraction  n. 


3.3  b  Definition  of  the  fundamental  integrals** To  facilitate  the  ap“ 


prc»Eimate  evaluation  of  the  vazloas  integrals  in  equations  (3i24a}  and  (3^^^) 
ve  shall  define  certain  fundamental  integrals  from  which  aH  others  could  be 
dexlved  by  differentiation.  We  define  for  oonvenienee 


tt. 


L 


(3*25a) 


*9 


(3*25b) 

(3.250) 


Itae  above  integrals  are  not  independent.  Indeed^  tbqy  ere  related  as  follovst 


u,  =  - 

*0 

2ia^Ut  «  (3k+  U, 


(3.26) 


idiat  follovBf  we  shall  avail  ourselves  of  the  above  relation  to  cheek  '^le 
dil^erentiabili^  of  our  results  after  integration  is  performed. 

ibQiloTing  the  above  definitions  of  U,  »  Uz  >  and  0^  we  oan  rewrite  the 
Hertilsn  veptor  oomppnents  in  tens  of  then  and  obtedin 


’"p 


,  »nym  k»  f  1 . 


^  (3.27a) 


and 
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4irkoU»/i«i. 


(3.27b) 


lijidyrloal  _6bQr(Uiiates  in  coo^^ 


tianafonn  apacea*»the  above  Integral  representatlonB  are  surface  Integrals 
over  the  entire  at,  -  or^  plane  Of  a  fom  ^lat  lends  Itself  readily  to  a  trana^ 
formation  to  cylindrical  coordinates  in  the  configuration  space  as  veil  as  in 
the  transfoxnatlon  space.  ThuSt  we  daploj  the  transfoxnatlon 


X  ^  f  cos  (f 
y  =  f  sln^ 


(3.28) 


fs  /x^  +  y* 
for  the  space  coordinates  and 

Ot,  =  1  008(3 

ot^si  ^  sln/3 


(3.29) 


^  =  V  +  Of** 


for  the  transforD  variables.  As  the  oonsaquenoe  of  the  above  transfonsatloni* 
ve  note  the  following: 


«iX  +  fl*j.y  ^  ^foo3(/3-^) 
#9  ^  /  k|  •  V 
s  s  /kf Txi 
doc.  dP(z  »  A  dl  d/3  . 


(3.3Q) 


The  Integrals  appearing  In  the  {oavlpus  seotions  are  of  two  different  types. 
The  first  one  Is 
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1 1  • 

-«o 

C®  T 

a  Jr  C>‘)  Ad  A  Je  e|(S 


(301) 


mm 

■2tjRa*)J.rAf)AdA 

o 

i4i«n  WB  used  the  uel^oHooua  x^pi^sehtatlon  of  Bessei  functions  {22,  p.  273) 

T 

.  (3.32) 


The  aeoond  dlstlnot  Integral  la  of  the  fom 

*9 

(3.33) 

•  ^T 


It  Is  show  In  Appendix  A  that 
t 

4¥c  f ; .  (3.34) 

jV  *'•* 

Substituting  this  Identity  Into  equation  (3.33) »  we  obtain  the  desired  result 

—  r 

h  ^  m  !§I(rirsin(2v*l)<p  j  ,  (3.35) 

O 

In  vhat  foUovs,  It  vUl  be  opn^enlent  to  have  the  range  of  Integration  from 
9  CO  to  ‘t’  ^  rather  than  0  to  .  Tb  this  wd  ve  Intrpduoe  the  follow^  vell- 
knoun  z'elatlon  betveen  the  Bessel  and  Haxdcel  funotipns  (26»  p.  75) 


(3,36) 
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Vhioh  eziables  us  to  reuilte  (3*31)  and  (3.35)  in  the  totm 

I,  sTtJ  F,  H‘;’{Xf)AdA 
and 


li  =  zm 


y  ^l)<^ 


^a*)wiU/Xf)dA . 


(3.37a) 

(3.37b) 


Ve  are  now  ready  to  rewrite  the  fundamentaL  integrals  U, 


^llndrioal  coordinates. 

I  t  o*‘^aK4SaZ) 

"  tjkmr 


We  obtain 

Hraf)Adx 


U2  t  W3  in  the 


(3.38a) 


vmo 


«<• 


CsK.-^a«x.) 


and 


(3.38b) 


(3.38c) 


The  integral  U,  can  be  recognized  as  one  originally  found  by  Sonmerfeld  in  his 
investigation  of  the  field  of  a  vertical  dipole  above  earth  (20).  The  Inte^ 
grals  U2  and  U3  are  someiduit  related  to  vhe  Sonmerfeld  integrals.  The  dlf«^ 
ferenoe  comes  about  from  the  presence  of  an  "s"  multiplying  the  denominator  in 
both  82  and  and  the  factor  (s  ^  n^  s„)  in  the  denominator  of  U^. 

Ve  also  note  for  futmre  reference  the  foUowinv  transfoznatlon  of  the 


partial  derivatives: 

fsoacf^  ^  ^ 

■  8in<p3if  +  ij.  3i^ 


(3.39) 


and  the  transfornation  of  the  components  of  the  fields 


(3.40) 
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(  )y  *  “•  (  )x  slntp  +  (  )y  oosf 
(  -  (  )h  ooa<f  +  (  )y  sintf  . 

3.3-1  Tt^anafoirntion^to  spherical  coordinates  In  confiL-aration  and 
Arahsforn  aPaces^For  oonvehlence  in  evaluation  of  the  fundamental  inte^rala 
hy  as^ptotic  methoda  ue  ahall  now  transfona  the  previous  expressions  to 
spherical  coordinates  in  both  tile  configLiration  ^d  transfoxii  spaces.  fhus> 
ue  employ  the  transformation 
z  -  r  cos  4 

f  =  r  sin  ®  (3.41) 

r  =  \/pnrp 

for  the  space  coordinates  and 

k^sin^  (3.42) 

for  the  transform  variable.  As  the  conaetjaence  of  the  above  transfoiiaation 
ue  note  the  foUouingt 


So  *  ko008(3 

s  =  ko  /n^  -sin*(i 
So  8  A  ^  =  kor  oos(|5  -  9) 
dl  ^  kg  cos(t  d/5  . 

These  transformations  enable  us  to  recast  U, 
forms: 


(3.43) 


U2  and  Uj  in  the  follouing 


9.  B  r  ei-n.flgoeeHSVk.ftw/S>e 


r: 

Uj  «  (rty'9Ln(2v*i)^ 

iimO 


(3.44a) 


(3.44b) 
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and 


J 

P 


(3.440) 


idiere  we  extracted  the  exponential  beha^or  of  the  Hankel  function  and  ve 
denoted 

H^(lc„^sin(S)  =  (3.45) 


The  original  path  of  integration  that  ran  from  *c^to  along  the  real 
axis>  has  now  been  transfonned  aocording  to  equation  (3.42)  into  the  path  17 
as  shovm  in  Fig.  3*1. 


We  also  note  for  future  referenee  the  following  transfornation  of  the 


partial  derivatives > 

00s  Q^r-  a, 

3^  *5  sin  6  +  J?gg-9  n. 

and  the  transfonoation  of  the  components  of  the  fields 


(3.46) 


(  )«  =  -  (  )j  sin  0  ♦  (  )-  00s  0 

’  f 

(  )^  -  (  )^  008  9  +  (  aln  9  • 


(3.47) 


Combining  the  results  of  (3.46)  with  those  of  (3*39)  ve  obtain  the  transfoma^ 


tion  of  the  partial  derivatives  from  the  Cartesian  system  to  the  spherical 
system  as  folio vs; 


2k~  OOSCP  sinOS^  +  ^ 

r  -  r  sin  9  ^ 

ev.  *Ln  0  a.  * 


(3.4S) 


cos  9  3r  «  3# 


Sinllarlyi  combining  the  results  of  equations  (3.47)  and  (3.40)  we  obtain  the 
transfoxmation  of  the  vector  components 
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(  =  coscp  cos  ©  (  )^  +  sin<f>  cos  ©  (  )y  »  sin  9  {  )^ 

(  )y  -  "  sincp  (  +  oos<p(  )y  (3«49) 

(  )y‘  ^  coscpsln  ©  (  -f  sln<psln  ©  (  )y  cos  ©  (  )£  • 

3.4  CLOSURE 

In  the  foregoing  chapter  we  appiUed  a  high  frequency  appro^d^tlon  to  the 
rigorous  formulation  of  the  present  boimdary  value  problesi  of  Chapter  Z*  This 
appro^djiiation  is  valid  for  all  \t&ve  frequencies  appreciably  higher  than  the 
qyclotron  frequency  of  the  plasma«  In  the  coafiguration  space  the  same  ap** 
proximation  is  ei^pected  to  give  good  results  for  ©11  polar  angles  except  the 
one  Gorresponding  to  the  angle  of  the  total  interaal  rejection  in  the  plasma. 

Introduction  of  the  high  frequency  approximation  into  the  expressions 
for  the  boundary  coefficients  resulted  in  a  great  simplification  of  the  latten 
Subsequently,  after  the  determination  of  Hertzian  vector  components  in  the  aii; 
it  was  found  possible  to  separate  the  contributions  of  the  plasma  anisotropy 
explicitly.  These  extra  contributions  display  some  heretofore  unencountered 
denominators  which  will  be  responsible  for  some  of  the  unusual  phenomena  that 
do  not  have  their  counterparts  vdien  dealing  vdth  Isotropic  interfaces. 


CHAPTER  4 


IVALUAHON  OF  THE  FlJNDAMMfAL  INTEGRALS  FOR  THE  OIPQLE  IN  MAQNETOPLASMA 

In  the  preceding  chapter  tie  found  the  approximate  formal  solutions  for 
the  fields  in  the  air  In  the  form  of  certain  definite  integrals^  In  order  to 
have  the  results  In  a  form  suitable  for  numerical  oalculations  these  Integimls 
have  to  be  evaluated^ 

The  subject  integrals  do  not  lend  themselves  readily  to  rigorous  evalua<!- 
tlon  and  in  order  to  obtain  any  useful  Information  from  them  one  has  to  in^ 
dvltably  resort  to  approximate  methods  of  their  evaluation  <  Fortunately 
enoughi  idien  the  field  points  of  interest  are  at  a  large  distance  from  the 
source  there  Is  an  approximate  method  available  that  is  knotin  to  yield  good 
results.  This  method  Is  referred  to  as  the  "saddle^point  method^”  or  the 
"method  of  steepest  desoents"  and  it  tias  first  Introduced  by  Debye  for  obtain¬ 
ing  aayiiq>totlo  expansion  of  the  Hankel  function  (Ilf  p.  503).  The  basic  pro¬ 
cedure  in  applying  the  saddle-point  method  is  as  follotis.  One  first  locates 
the  statlonaiy  or  the  saddle  point  of  the  Integrand  under  consideration.  Nexl^ 
the  original  path  of  integration  is  deformed  to  a  path  going  through  the  sad¬ 
dle  point.  A  nev  path  of  Integration  is  defined  on  tihioh  the  integrand  de¬ 
creases  most  rapidly  \dien  moving  avray  from  the  saddle  point.  When  this  is 
aecompUshed  the  integral  has  been  recast  in  a  fozm  to  uhiCh  the  Watson's 
Lemma  is  applLoable  (Ilf  p.  301).  One  then  proceeds  expanding  the  integrand 
in  Taylor's  series  about  the  saddle  point  and  integrates  tern  by  term.  The 
result  is  an  asymptotic  series  in  teras  of  the  Inveree  distanoe  from  the  souroa 

As  tie  remarlEed  beforef  the  application  of  the  saddle-point  method  is  con¬ 
tingent  upon  being  able  to  contlnuoui^  deform  the  original  path  of  Integra¬ 
tion  to  the  path  of  steepest  descent  passing  through  the  saddle  point.  This 
ISf  of  ooursef  alvsys  legltlmste  by  Cauchy's  theprem  providing  there  are  no 
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singulfiorltles  between  the  two  paths.  If  there  are  singularities  between  the 
two  paths y  a  correction  to  this  transformation  mast  be  made  by  adding  the  con¬ 
tribution  from  the  poles  and  branch  points  of  the  integrand. 

The  integrauids  occumlng  In  the  electromagnetic  problems  usually  contain 
a  Hankel  function  of  one  kind  or  euiother.  Before  the  saddle-point  method  can 
be  applied  one  must  extract  the  exponential  behavior  from  the  Hankel  function 
to  determine  the  location  of  the  saddle  point.  This  is  oftentimes  done  by 
substituting  for  the  Hankel  function  Its  asymptotic  formula.  Whereas  the 
procedure  Is  basically  soundy  it  introduces  an  unnecessaiy  singalarlty  and 
makes  the  results  invalid  in  the  neighborhood  of  small  polar  angles.  The  way 
to  avoid  this  was  devised  by  Banos  and  Wesley  (2,  Oh.  5)y  tdioy  Instead  of 
using  the  asymptotic  formula  of  the  Hankel  function,  make  use  of  Its  integral 
representation  which  yields  a  double  Integral  to  evaluate  rather  than  a  single 
one.  By  formal  extension  of  the  Watson's  Lemma  to  double  integrals,  they 
formed  a  rigorous  basis  for  the  use  of  the  saddle-point  method  for  doubla 
Integration. 

In  this  thesis  we  shsdl  avail  ourselves  of  the  saddle-point  method  for 
double  Integration  of  Banos  and  Wesley  (2,  Cb.  5)  not  only  for  the  reason 
mentioned  above  but  for  the  following.  In  our  case  the  basic  integral 
contains  aa  Infinite  sum  of  Hankel  functions.  Whereas  the  asymptotic  formulas 
for  the  Hankel  functions  of  order  aoaller  than  the  argiUDent  are  simple,  the 
sane  is  not  true  in  the  reverse  pass.  Using  the  integral  representation  of 
the  Hankel  function  in  our  case  avoids  then  these  additional  natbenatioal 


4.1  SINGULARITIES  IN  TH£  IS -PLANE 

The  evaluation  of  the  basic  integrals  U,  ,  U2  >  and  U3  by  the  saddl*- 
polnt  BSthod  of  integration  is  neoessarlly  oontingsnt  being  ahls  to 
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deform  the  original  path  of  Integration  T  to  the  path  of  steepest  descent  f* 
In  Fig.  4*2«  In  the  process  of  the  continuous  path  deformation  we  must  be 
sure  that  any  singularities  between  the  two  paths  are  properly  accounted  for. 


To  this  end  we  shall  need  to  know  the  exact  location  of  the  -/arlouS  slngularl** 


ties  of  the  integrands  of  U, ,  In  the  complex  -plane. 


4.1  a  The  location  of  the  Doles— »The  examination  of  the  integrands  of 
U,  i  U2  and  O3  In  equations  (3.44)  re\reals  a  possibility  of  a  pole  at  a  point 
where  the  denominator  vanishes.  This  occurs  where  and  satisfies 

the  equation 

/n*  -  si*^  +  n*  ooa^f  *  0  »  (4«1) 

The  above  equation  can  be  readily  solved  for  |Sp  giving 

sin /ip  a 

/I  +  n* 

To  determine  the  location  of  this  pole  In  the  complex  /t -plane  we  assume  that 
n  is  complex  with  a  small  Imaginary  part.  We  put 
n  B  n'  +  i  n*  }  0  <  n"  <  n' 

or  (4*3) 

ns  Ini  e'-f  ;  0<<f 

then 

yTTnTs  M,e^®  o<e<«f  (4«4) 


For  convenience  in  \diat  follows  we  shall  Introduce  at  this  point  the  concept 
of  memaiua  sheete*  Vie  shall  call  the  Riemann  sheet  on  ^eh  ln^(n^-sin^^)^j><> 
the  upper  sheet)  and  the  one  on  >Mch  Im  |(n^-8ln-fe  }^|<P  the  lower  sheet* 
Coasequentlyy  we  write  for  the  upper  sheet 


and 


/n*  -»8in% 


M, 


(4.5a) 
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VOa*  *ain%  =  ^^ssssr-  — I— . 

✓l  r  n*  M, 


(4.5b) 


for  the  lover  sheet,  thus,  in  order  to  ha'/e  a  pole  on  the  upper  or  lover 
sheets,  the  equation  (4.I)  must  be  satisfied  together  with  condition  (4.5a) 
or  (4.5b)  for  upper  or  lower  sheet  respectively.  We  obtain 


im*  e'-''’  oos/3p=  0 

(4.6a) 

for  the  upper  sheet  poles,  and 

-  +  Ini*  e*"*^  cos^p  »  6 

(4.6h) 

for  the  lover  sheet  poles,  the  above  equation  can  be  put  in  a  more  convenient 
form  if  ve  denote  ‘f  1  .  this  gives  for  the  upper  sheet 


dos  /ip'  cosh  /ip"  =  *  -1%*  --  00s  Q  (4.7s) 

M 


and 


sin  ^p '  sinh/ip" 


sin  d 


(4.7b) 


From  equation  (4* 7a)  it  foUovs  that  "^/Z  <  (i*  <  ^/Z  .  Using  this  fact  to¬ 


gether  vith  equation  (4.7b)  we  arrive  at  the  location  of  the  upper  sheet 


poles  as  follovsi 


p,  I  ■"■/2  <  (s;  <'rf  ^"<0 

Pj  j  IT  <  (s;  <  ^/z  (t"  >  0 

For  the  lover  sheet  ve  obtain  fr<nn  (4.6b) 

oos  3p'  coi^  (ip*  »  -jl”"  ^  ^ 

\n|^ 

sin  (S^'  sinh  /ip"  »  ® 

vliiph  gives  the  follovlng  for  the  location  of  the  lover  sheet  poles 
P,  I  0  <  |3'  <  ’>72  (S^  >  0 

P4  I  *  Tr/2  ^  /d*  <  0  lb*  <»  0 


(4.8) 


(4.9) 
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the  location  of  the  upper  and  lower  sheet  poles  Is  depicted  In  Pig.  4.I  lAien 
lnl<|  which  Is  of  primary  interest  to  us*  We  shall  ha'/e  the  occasion  to 
return  to  the  subject  of  pole  location  later  in  this  thesis  idien  discussing 
the  evaluation  of  the  integrals  U,  y  Wj »  • 

4.1  b  „Ihe  branch  tx)lnt»‘*the  integrands  of  and  U3  all  Contain  the 

radical  (n^  «aia^^  as  a  result  of  which  the  point  vhete 

±  arc  sin  n  (4*1^) 

idU  be  a  branch  point.  At  each  point  /S  the  Integrands  in  U, ,  and  U3  can 
take  one  two  values  depending  on  vdiich  si^  we  choose  for  the  radical. 

As  before  I  it  is  convenient  here  to  talk  about  two  sheets  of  the  (S-plEUie 
(formed  by  a  two-^sheeted  Memann  surface)  on  which  each  Integrand  is  single 
valuedi  then  according  to  our  previous  convention  we  call  the  upper  sheet 
one  on  which  Im  |(n^-sin^/3  \  >  0  and  the  lower  Sheet  one  on  which 

Im  £(n*’-sln^(t  1  convergooce  of  the  integrals  is  assured  if  the 

path  of  Integration  at  least  begins  and  ends  on  the  upper  sheet.  The  two 
sheets  will  be  joined  along  the  lines 

Im  ^(n*-sin*|ft  )i'  j  »  0  (4*11) 

starting  at  the  branch  points.  These  lines  are  shown  as  dashed  lines  in 
Fig.  4.1.  It  can  be  shown  that  the  In  {(n^»8in^^  )  ^  ^  ^  ^ 

dashed  line  facing  the  origin  and  ^  £(n--!>sin^(}  )  ^  <  0  on  the  opposite  side. 

In  the  following  analysis  it  will  be  necessary  for  us  to  determine  the 
location  of  the  poles  P3  and  more  accurately.  In  particular  it  will  be 
essential  to  know  whether  the  poles  P,  and  p^  lie  below  or  above  the  outs 
issuing  from  the  branch  points  B,  and  B2  respectively.  For  this  purpose^  in 
addition  to  the  cut  represented  by  the  dashed  lines  and  corresponding  to 
Im  ^(n^<«sin*(t  }i  ^  s  0,  we  constraot  lines  Be  |(n^..sin^(S  0*  The  latter 
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lines  are  showi  dotted  In  Fig.  4*1  •  signs  of  Id  ^(n^-sin^/%  )^j  and 
Re  I  (n^-sln‘/%  }  j  ohange  only  uhen  one  of  these  lines  is  crossed.  It  follows 
from  equation  (4* 5b)  that  both  signs  are  negatl/e  on  the  lover  sheet.  Tbere^ 
forSf  the  poles  and  P4  are  on  the  same  side  of  the  cuts  as  the  origin. 

4^2  FORMUL\flGN  OF  THE  C0RfRlBUtl0.'ld  TO  THE  FUNDAMfUTAL  iN'fSGRALS 

In  the  pre^^ous  section  we  detexmined  the  singularities  of  the  Integrands 
of  U, »  (It)  and  U3.  In  particular  we  found  a  branch  point  and  poles  on  both 
Rlemann  sheets.  In  this  section  we  Shan  determine  the  various  contributions 
to  the  fundamental  integrals  as  affected  hy  these  singularities. 

Any  one  of  the  basic  integrals  can  be  written  in  the  form 

(I  =  I  F  ((3  )  Hi,”  (kof  aia$  )  e  U•^2) 

<7 

and  we  wish  to  evaluate  it  approximately  when  r  is  large.  The  saddle  point 
of  the  integrand  occurs  %dien  the  derivative  of  the  exponent  vanishes*  l.e.* 
^en 

<3  s:  0  (4.13) 

The  path  of  the  steepest  descents  is  found  fToa  the  equation 

cos  ((S-0)  s  1  +  ir*  (4.14) 

where  r  is  a  parameter.  It  can  be  readily  shown  that  this  path  will  intersect 
the  real  axia  at  the  point  9*  the  saddle  point*  at  an  angle  and  it 
vUl  oorrespond  to  the  path  T  as  shown  in  Fig.  4*2.  Ue  now  investigate  the 
possibility  of  replacing  the  original  path  of  integration  C  with  the  path  P 
oorresponding  to  the  path  of  the  steepest  descent  tbrou^  the  saddle  point 
~  Since  the  poles  P3  and  F4  lie  on  the  lower  sheet  they  do  not  have  to 
be  taken  into  aeopunt  when  shifting  ^e  path  of  integration  from  H  to  r  • 

Ihs  upper  absst  poles  F|  and  P.  need  not  be  talpsn  Into  aoeovut  either  sinoe 
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they  \)ill  never  be  crossed  lying  beyond  the  region  of  interest,- ^2  <  fi  <  ^/2. 
Thus,  we  are  free  to  replaoe  the  original  path  of  integration  H  by  the  path 
of  l^e  Steepest  descent  P  . 

Nov  ve  focus  our  attention  on  the  foUouing.  As  ve  already  remarked 
earlier  the  integrands  of  U, ,  U^f  and  U3  are  double-valued  since  they  contain 
the  radical  (n^-sin^/t  )!'.  Tke  original  path  of  integration  P  passes  over  the 
upper  sheet  of  the  corresponding  Riemann  surface  and  can  be  deformed  into  the 
path  of  the  steepest  descent  P  only  when  at  least  the  beginning  and  the  end 
of  it  lie  on  ^s  sheet,  in  the  case  when  the  angle  of  incidmce  d  does  not 
exceed  ^  angle  of  the  total  intend  reflection  65  ,  l.e.,  tSien 

-  Og  <  e  <<  Og  (4.15) 

then  in  the  complex  -plane  we  have  a  picture  as  shown  in  Fig.  4*2.  Here 
the  path  of  the  steepest  descent  (nrosses  the  branch  cut  twice  and  the  transi¬ 
tion  from  the  path  p  to  the  path  P  is  aocomplisbed  without  complications. 
£)coept  for  the  part  shown  dotted  in  Fig.  4>2,  the  path  P  will  lie  in  the 
upper  sheet.  In  this  case  the  only  contribution  to  the  integrals  U|  ^  and 
U}  will  be  the  one  from  the  saddle  point. 

The  situation  is  different  when  the  angle  of  incidence,  0,  exceeds  that 
of  ^  total  internal  reflection,  •  Here  the  path  of  the  steepest  descent 
*411  cross  the  branch  cut  only  once,  crossing  from  the  upper  sheet  to  the 
lower  sheet  and  going  to  infinity  along  the  lower  sheet,  without  connecting 
anywhere  with  the  path  P  which  is  inadiidssible.  He  can,  however,  construct  a 
sore  complicated  path  of  integration,  supplanenting  the  path  of  the  steepest 
descent  P  by  a  contour  encompassing  the  out  in  such  a  way  that  the  beginning 
and  the  end  of  the  more  cosqilicated  path  will  lie  on  the  upper  sheet. 

To  this  end  ve  introduce  "cuts’'  in  the  complex  (5 -plane  along  the  branch 
cut,  i.e.>  along  the  lino  on  whi^  1*  |(n*-ain^|a  )i’j-  Q  wbieh  ere  oarked 
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dotted  in  Fig.  4.1  and  4*2.  The  equation  of  the  outs  uLll  be 


n^  -  sin^/S  s=x^  (4*l6) 

and  aecording  to  the  previous  analysis  In  {(n^>sln^if3>  )^j  >  0  on  the  sides  of 
the  Outs  facing  the  origin  and  Im  ^  (n^-sln^^  <  0  on  the  other  side. 

For  0  >  d  this  nev  path  of  integration  vhleh  vas  first  proposed  Ott 
(13,  p.  451)  runs  along  the  line  path  of  the  steepest  descents  r  then  along 
the  borders  of  the  cut,  Intersects  the  branch  cut  once  thus  entering  the  lover 
sheet.  It  then  continues  along  the  path  r  and  intersecting  the  branch  cut 
once  acre  it  ends  on  the  upper  sheet,  fhe  remainder  of  the  path  is  along  the 
path  r  on  the  upper  sheet,  the  situation  is  depleted  in  Fig.  4*3.  For  9  <  0 
the  situation  Is  analogous. 

As  a  result  the  complete  expression  for  any  of  the  Integrals  0, ,  0^ ,  or 
vLU  ooMlst  of  two  parts  and  can  be  written 

|j  s  U»»+  o«u  (0  -  ds  )  (4^17) 

where  U^*^vilL  stand  for  the  contribution  from  the  saddle  point  and  U^^vUl 
stand  for  the  contribution  from  the  integration  about  the  borders  of  the 
brandi  cut.  Since  the  contribution  from  the  branch  cut  does  not  appear  until 
|id|  >  |dfl  >  we  hare  anployed  a  unit  step  function  u(0  -  Ob  )  to  denote 
this. 


4.2  a  Fomuistion  of  the  contrlbu^on  from  the  saddle  tx>lnWAs  we  re^ 
narked  earlier  each  one  of  the  fundamental  integrals  o«a  be  witten  in  the 
fom 

y  *  F  ( (3  )  Ht’  (k,9  8in(5  )  e  d/S  (4.18) 

and  the  saddle  point  of  ^e  integrand  is  located  where  the  derlrative  of  the 
sKponent  yanisheaf  l.e.,  ^re  0.  First  vs  put 

w  -  (3  -  «  (4.19) 

whioh  transfers  the  saddle  point  to  the  point  v  ^  Q  in  the  wmplac  v  <»plane. 


I 


Figuov  4*3  ^  P^onwd  Pith  of  |nt«praUon  For  |0|  >  Ih,) 


1 
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Following  the  procedure  of  Banoa  end  Uealey  (2«  p.  97)  we  naJca  the  next  trans- 


foimetlon  ea  follovsi 


ik^r  (  1  -  ooa  v) 


(4«20) 


noting  that 


d  z  i  kozain  w 


(4.2i: 


tte  alao  note  that  the  right^Btuid  aide  of  aquation  (4*20)  oan  he  expanded  into 


a  power  aerlea  aa  followa: 


=  W^  {  +  O^  +  0^  +  '  • 


}  (4.22) 

}. 


the  above  aerlea  can  be  Inverted  (2,  p.  134)»  ^vlng 


(i^o  r)*  24  (ir)* 


(4.23) 


Fran  the  above  we  note  that  w  la  lan  odd  funotlon  of  x,  w(x)  ^  -  v(-x).  lhaa» 
aa  z  obangea  sign,  ao  doea  w. 

Bov  we  oan  tranaform  the  Integral  In  equation  (4*10)  to 


Ufffm  9ck,r-  I  Q  d  35 


(4.24) 


where 


G(x)  ^{pCW'+ejMjI’Ck.f sin^w+«))  + 


(4.25) 


Patting  thia  expqresalon  under  the  integral  algn  in  (4.24)  reaulta  In  two 
integrala  with  the  Integration  range  froD  0  to  ^  •  In  the  aeoond  Integral 
we  Bake  the  aubatitutlon  x  s  .  x'  and  reverae  the  llBltat  One  to  the  fact 
that  0v/dx  la  an  even  function  of  x>  we  can  add'  theee  two  integrals  to  fom  a 


single  one  as  followsi 


Jp(w)HvCh«^stia(>y^O))€’'*^  dx 


(4.26) 
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vhere 


sin  V 


(4.27) 


Hdv  vs  nake  use  of  the  integxttl  irepTOsentation  of  the  Hankel  funetlon  (I^,p.ld6) 


“v  («) 


vhioh  m\i  we  put  into  equation  (4.26)  to  obtain 


•4#  4 


(4.29) 


i^ere 


A  /«  «)  s  xy*>Frwy.9)l4<;k40St>i,^w^^ 


iw-Ji 


(4.30) 


Ue  recognize  the  integral  (4.29)  as  being  analogous  in  forn  to  single  integral 
of  the  type  to  which  Viatson^s  Lemma  can  be  applied.  A  theoretical  basis  ex¬ 
tending  the  Watson's  Leutna  to  double  integrals  has  been  provided  by  Banos  and 
Wesley  (2«  Cb.5).  Thus,  at  this  point  all  that  remains  to  be  done  to  find  the 
saddle-point  contribution  to  the  integral  (J,  is  to  expand  the  integrand  Q(x»y) 
in  a  double  power  series  and  integrate  term  by  term.  When  perfpndng  such  an 
integration  we  shall  a/ail  oinrsel^res  of  the  following  well-koown  integrals 
(9»  P.  64) 


(4.31a} 


and  in  general 


find  the  eontxlbutlen  to  the  Integral  (4*12)  due  to  the  Integration  about  the 
borders  of  the  cutS|  ve  first  oonsider  the  oaaa  of  6  >  0^  >  0  oorraapoiuUng 
to  the  upper  out.  We  inrite 

Jp  (/S)  Hi’  (k,p  eln/d)  e  d/e  (4*32) 

end  integrate  forBaUy  on  both  sides  of  the  out,  getting 
a«=  jp4  (^)  (k.^sia(i) 

+  |p4A)  H|;*(k.^8in|9)  a  tWv-cpso-e)  dyi  (4.33) 

e. 

vhere  F't-((t )  is  the  value  of  the  function  P  {(i)  on  the  left  side  of  the  out 
and  F.(  (3 )  is  the  value  of  the  same  function  on  the  right  side  of  the  out. 
Interchanging  the  lisdts  of  integration  in  tile  first  integral  vs  obtain 

y<*%  |r  (fl  )  Hi.»  (k.y  8in(a)  d/8  (4*34) 

a. 

vhere 

R  (/a)  «  F,(/5)  •?+(/&)  ,  (4.35) 

We  note  that  since  the  imaginary  pcyrt  of  (n^«*  8ln^y5)'$’  in  sere  on  the  cut  and 
the  real  part  has  different  eigne  on  the  two  borders  of  the  out,  then  F-f ( /3 ) 
end  F.(  0 )  will  differ  fron  one  another  only  in  the  sign  of  this  square  root. 

In  idiat  follows  it  will  be  oonvani^t  to  deform  the  path  of  integration 
in  sttdi  a  way  that  it  goes  f^  the  branefa  point  B,  along  the  line  on  vhi^ 


65 


th«  exponent  in  the  Integrand  decreases  most  rapidly.  This  id.ll  be  the  line 
on  iMeh 

Re^eoe(/&-  0)|  B  const  (4*36) 

Then  it  is  naeessaxy  that 

In  ^co8(/&- 9)  ^  >  0  (4*37) 

Ue  assume  for  convenience  that  n  is  real.  Then  setting  ^  s  in  equation 
(4*36)  giires 

const  B  cos  (  9^  -  8)  (4*36) 

Mov  since 

008(0  *  9)  *  cos  (0'-  9)  oosh/3  ■  +  i  Sin(9  -  ft>')  sinh "  (4*39) 

we  obtain  for  the  path  of  integration  from  equations  (4.36)  and  (4*38) 

oos(  (3  -  9)  cosh  (3  •  =  0O8(9a  -  9)  ,  (4*40) 

This  path  is  shown  in  Fig.  4.4  as  a  solid  line.  The  condition  (4*37)  is  also 
satisfied  on  this  path.  The  defonnatlon  from  the  path  Q  to  r^'  presents  no 
dlffiottlly  since  there  are  no  singularities  between  tbem.  In  particular,  the 
pole  P)  llee  below  and  not  above  the  line  as  we  remarked  earlier  in 
Section  4«1« 

In  view  of  the  results  of  equation  (4>40)  we  rewrite  equation  (4*34) 

I  )  go*  (4.41) 

Now  using  equation  (4.4O)  ve  express  the  entire  Integrand  in  terms  of  ** 
which  we  treat  as  a  new  variable  of  integration.  The  new  Hjiits  of  integral 
tion  are  from  0  to  <>o  .  For  k,r  large  on  the  path  of  the  steepest  descent 
only  small  values  of  "  win  be  important.  Moreover,  we  can  set  '  b  0 
on  the  initial  part  of  the  path.  Thus, 


ud 


Sin(9  9(3')sinh(3  *  aln(9  «  9g  ) 


(4.42a) 
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d(3  =  i  d3  * 
Finally,  we  put 


and  rewrite  (4^41) 


#J=  i e  ^  R ( (3  )  sin/S ) 


(4* 42b) 

(4.43) 

(4.U) 


l&ere  ft  in  the  abo\re  equation  stands  for 

,3  =  9,^  (4.45) 

Next  we  perforiB  the  Integration  around  the  borders  of  the  branch  ait  origina¬ 
ting  at  pertinent  to  angles  e  such  that  -  ^/2  <  ®  .  Integrating 

formally  we  get 
B* 

|f.(  ,5 )  H  sin^a ) 

-c*o 

+  I  F+( (i  )  h';’  sin (3 )  6  <:k.rcos(|a-e)d^  (4.^) 

where  again  F^((3  )  is  the  value  of  the  function  F  (13)  on  the  left  side  of  the 
cut  ^  F«(  (3  )  is  the  value  of  the  same  function  on  the  right  side  of  the  cut. 
Interchanging  the  limits  of  integration  in  the  first  integral,  we  obtain 

-.e'eo 

U^=  -  J  R((3  )  (kp^  3in(3)  (4.47) 

i^ere  R(^  )  is  the  same  as  in  equation  (4.35)  •  Proceeding  as  before,  we  de¬ 
fine  a  new  path  of  integration  on  which  Re  |^cos(|t-  0)j  -  const.  But  tliis 
Use  we  have  to  set  ^  -  0^  in  equation  (4.36)  to  obtain 

const  ss  cos(9b  t  9)  (4.48) 

and  the  new  path  of  integration  is  expressed  by 
Pos(  /&'  •  9)  cosh  (5"  =  CQs(9g  +  9)  , 


(4.49) 


68 

This  path  is  shown  in  Fig.  4*5  as  a  solid  lino.  In  viav  of  tha  aboye  resti^lis 
we  rewlte  equation  (4.47)  as  follovrat 

^  I  a(  ,4  )  Ht’(k^P  sinis)  (4,50) 

<i‘ 

iiow  using  equation  (4.49)  we  express  the  entire  integrand  in  tefas  Of  /3  " 
tdiioh  we  treat  as  a  new  variable  of  integration.  The  new  Malts  of  integra¬ 
tion  are  from  0  to  -  ^  .  Moreover,  we  set  /s  •  ^  9^  d&  the  <ni».iar  path, 

nius, 

sin(o  -  /3‘)  sinh/s"-^  sinidg  +  ©)(&"  (4.5la) 

and 


d/i*  id(3  "  , 


Nov  ve  rewrite  (4,50) 

a<M=  -  iaikorcoicee*e)| R(^  )H<;>  Sln^ )  ^ 


(4.51b) 


(4*52) 


It  will  be  more  convenient  to  have  the  range  of  integration  from  0  to  <f 
rather  than,  from  0  to  —  &9  }  to  ^lis  end  we  put 


/5" 


(4.53) 


which  gives 

i,ck.»*cpsO#i.e^  I  (4.54) 

o 

i^ere(h  etande  for  the  same  as  in  equation  (4.0).  Noting  that 

?eko®  (4.55) 

ws  can  include  (4*41)  and  (4.55)  in  a  singie  equation  as  foMovst 

i*^****"®®*^'*'”®^  I  *R(t/5)HW(k,r^|«l  ein(S) 

p 


(4.56) 
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-/3“ 


j  Figure  4.5  Sefonned  Path  For  Branoti'>Gat  Integration  VIhen  »  ^2  <  0  <  9p 

i 

I 

i 

j 


I 
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Nov  ve  Insert  tNe  integral  representation  of  the  Henkel  funotion  and  obtain 
finally 

ll(»>  2  - g- — - 

Vo  /T  z*^  f(v*yt) 


R(±A)L4Lki^r6uix.\^\Sri*^^  ‘^  y^2 
Vk#  ^5  16 1 S  tVfciil*' 


«''^y^4*57) 


We  pause  for  a  motrant  to  think  about  the  slgnlfieanee  of  the  above  re« 
suits. 

In  the  foriaulation  of  the  braneh^Mt  integration  ve  introdueed  a  change 
of  the  viable  in  the  foxn 

/S  *  «6  +  -1^  (4.58) 

and  subsequently  concluded  that  only  small  values  of  x  vould  contribute  to  the 
integral.  While  this  statement  in  itself  is  true,  it  produces  an  undesirable 
oomplicatlon.  For  if  fi  is  given  by  (4*58)  then  ve  are  forced  to  Integrate  in 
a  small  neighborhood  about  s  =  0  vhi^  is  Just  the  point  ^ere  the  high  fre^ 
quenoy  approximation  (3.7)  is  not  ralid. 

Whether  the  results  produced  by  such  an  integration  would  still  make 
sense,  it  is  difficult  to  appraise  at  this  time.  However,  in  a  more  rigorous 
analysis  which  ve  are  able  to  oarxy  out  in  the  case  of  the  line  source  pro» 
biens.  Section  9*6,  Indieate  that  the  results  of  such  an  Integration  are  inf*- 
deed  not  valid  since  they  produce  expressions  substantially  different  Cr<m 
those  that  are  obtained  when  the  high  frequeni^  epprpadJaa^n  is  introduced 
after  integration. 

We  shall,  henceforth,  neglect  the  contribution  from  the  branch-cut  inte^ 
gratlon  altogether  In  the  present  paroblsm*  First,  becsnse  ve  ere  not  abla  to 
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effeetlirely  earzy  It  out  in  the  present  formulation  euid  second*  as  ve  shall 
See  later  In  Section  9*6,  such  a  braneh<>eut  Integration  shows  that  the  lead¬ 
ing  terms  of  the  lateral  field  are  of  second  order  only. 

4^3  E^ALtliflOIi  OF  Tm  FmmimtAL  int&qrals  at  the  saddle  point 

In  Section  4.2  %ie  derlv^ed  a  general  formula  for  evaluation  of  the  saddle- 
point  contribution  to  the  fundamental  Integrals.  In  this  section  we  shall 
apply  these  formulas  to  find  the  a^jfnptotlc  expressions  for  ll,  *  *  and  U3  • 

4. 3  a  The  inteeral  J  ■ _ — Apolvlng  equations  (4*29)  and  (4*30)  to  the 

SKpression  for  U,  m  equation  (3.44a}  we  obtain 

so  <*• 

®  '  ~  J  /  (**y)  *  d  y  (4.59) 


where 


«  /.  H . 

iinw  s  Crv  Cw  -  e)  -y^  n,*-  s  i  rv*(  w  +  e )  ■  +  la^CpS  C  W'*’  S )  ] 

and  if(x)  is  given  by  equation  (4.23)*  Ue  find  for  the  leading  term 


Q  vl  4  •  ^  ' 

a  n.^cos9  *• 

Patting  this  into  (4*  $9)  wd  using  the  Integrals  in  (4*31)  we  obtain 
'  k#v“  *  "0’**7»x^?j55^3r 


(4.60) 


(4.61) 


(4.62) 


4.3  b  The  integral  Ux  — Applviog  equations  (4.29)  and 
eocpression  for  U2  in  equation  (3*446)  we  obtain 


to  the 


(^n''sL».rzv*09e 


»mp 


2*»*rvri  )**^'  r(Zv 


T^ffckC- 


(4.63) 
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*  ’"  sin  »vfsu-tCw^®)3**'*Vn^5t^^^®r^^^*?^^5e?+  rt^ebsC^w+Oj} 


and  v(x)  la  given  by  aquation  (4«23)>  Ue  find  for  the  leading  texin 


Qi(x,y)- 


/K^-Si/vi,*©  ^✓»l*-Stia*d  ‘►n^osd) 


Putting  this  into  equation  (4*^3)  gives 


0 


UmO 


2*"  V<^i<o»*5<i^<9  **»»1i*COSa) 

•J  f y^z  ^>^<dy  , 

— **o 

Using  (4*31)  we  note 

•  r(2U^  ^ 

e-^3'2  dx  .y^ 


/- 


iddOh  upon  substitution  into  (4*66)  gives 


_ 2  eo»Qg*'*‘’**^*^”*"*^*^  e*^****^ 

^  5tKi  9  /n^-sta^a^v^^^^rt^G^  +  Ki^cps©) 


0,  *  ^ 


VbQ 


(4.65) 


(4*66) 


(4.67) 


sin(Zv-rl)(f  .  (4*68) 


The  sari.es  on  the  extraiM  right  in  the  above  equatioB  oan  be  smswd. 
wite  foxmlly 


He 


^  8in(2y»1  )<f  ^  Ip  ^e‘'‘f  ^  j 


y-p 


and  recognise  the  series  on  the  right  as  the  binoalal  series  aossisting  of 
taxps  (e^^^)’'.  Thusi  «e  oan  svai  ipsediateiy  fay  vriting 


1 


2  rihy  • 


(4.69) 


5^ 
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InMrting  th#  above  result  in  equation  (4«68)  we  obtain  finally 
^  k*r  Si*\B  Stn<f  7??-3S??J^^^'s£r?e'^^*CO»0  ) 


(4.70) 


4.3  o  the  integral  J  a  ~ApDlyiDg  equations  (4*29)  and  (4.30)  to  the 
expression  for  U3  in  equation  (3.44e)  ve  obtain 


I  .“».*>'?»  <l*«f 


(4.71) 


•  w#  6 


idiere 


js  /  \  , X  s£n Cv</ >  Sjcgs/ w-»  <9) [4fck«SLiaf w*  e)*  y  ** *** ^  /^  ~  *_*' ^**^ /  ,  m\ 

*  si*xvJ 

where  v(x)  is  given  by  equation  (4*23).  We  obtain  for  the  leading  ten 


O  /v  v1  ~  — 


i  k.k  y  vib^ski^e 


(4.73) 


Substituting  the  last  result  into  (4*71)  and  aaking  use  of  the  integrals  in 
(4.31 )«  we  obtain 

co56g^****^C^***^ 


0» 


(4.74) 


4.4  OIFFERENTlAfilLITX  OF  IHB  iSXMPTDTlC  FORMS 

The  usefi^sess  of  the  asgraptotio  fons  for  U, »  U^f  and  0,  found  in  the 
previotts  seotion  is  necessarily  oontingent  upon  whether  these  fons  are 
differentiable  or  not.  In  this  section  we  shell; show  that  our  asyspiotie 
fons  are  differentiable  by  bhe  way  of  showing  that  the  differential  relations 
ssistlng  aaong  their  integnl  representations  ere  also  satisfied  by  Uieir 
asyiiptotio  fons* 

It  can  be  idiowD  that  the  following  differential  relationships  are 


<0 
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satisfied  bgr  the  integral  representations  of  U, »  »  and  U3 1 


0,  3^3,0^ 

*0 


(4.75a) 


2  i  3y  Oi  s  (  +  n^cie)  0 4  .  (4.75b) 

In  vbat  follovs  ue  shall  show  that  the  above  relationships  are  satisfied  bgr 
the  asgpBptotio  foxu  of  the  oonti^butlons  from  the  saddle  point* 

te  dinonstrate  the  differentiability  of  0,  1  02  »  and  0^  1  we  first  note 
using  (3.40)  and  the  results  of  Section  4*3} 


)  X-  ikooost^sin  ©  (  ) 

(4.%a) 

3y(  )  x.*  ikealn<fsin  ©  (  ) 

(4.76b) 

i^(  )  x.«  i  koCOS  ©  (  ) 

(4.760) 

3l^(  )  x.-  i  ko\/a*  *iln*  ©"  (  ) 

(4.76d) 

where  the  first  three  relationships  are  valid  only  to  the  first  inverse  order 

of  r.  Applying  these  to  (4.75a)  and  (4.75b),  we  obtein 

0,  B  2iBin<f  sin  ©  /n^*  sin^  ©  ©z 

(4.77a) 

2  i  sin<f  sin  ©  ©2  *  (  /n^  -^^©  <f  n'ooe  9)  0,  . 

(4.77b) 

An  examination  of  (4.^)  and  (4*68)  reveals  that  the  relationship  (4*77a)  is 
satisfied.  Similarly,  an  exandnation  of  (4.68)  and  (4*74)  reveals  that  the 
relationship  (4.77b)  is  satisfied. 

This  Qonoludee  the  proof  that  0,  •  0^  »  end  03  are  differiintlable  to  the 
first  inverse  order  in  r. 


4*6  CLOSURE 


Iq  the  foregoing  chapter  ve  evaluated  approximately  the  fundamental  Inte'^ 
grals  U, ,  Jj.)  ^3  •  ^  process  of  their  evaluation  ue  found  that  the 

as^ptotlc  representation  of  each  of  these  integrals  consists  of  two  basic 
parts*  The  first  part  comes  from  the  contribution  of  the  saddle  point  of  the 
integrand  and  represents  the  well*»laiovm  radiation  field  vdth  -  the  leading  term 
eharacterized  by  the  factor  r'^'  exp  (ik^r)*  The  second  part  which  appears  for 
d  >  @^  ^ere  corresponds  to  the  angle  of  the  total  intex^al  reflection  in 
the  plasma,  comes  from  the  integration  along  the  borders  of  the  branch  cut  and 
represents  the  lateral  field  (5,  p*  270).  while  the  contribution  of  the 
branch**cut  integration  was  fonaulated,  it  was  not  used  to  find  the  actual  coSk* 
ponents  of  the  lateral  field.  The  reason  for  this  was  the  fact  that  the  main 
contribution  to  the  lateral  field  appeared  to  have  come  from  the  point  where 
the  earlier  introduced  high  frequency  approximation  was  not  valid.  The  sub<* 
ject  of  lateral  waves  will  be  discussed  later,  however,  when  dealing  with  the 
line  sources  where  rigorous  evaluation  of  the  lateral  field  is  possible. 

Finally,  we  remark  that  the  aeymptotle  representations  of  the  fundamental 
Integrals  U,  ,  Uj,  and  Uj  are  differentiable  and  thus  are  valid  representations 
and  as  such  they  can  be  used  to  find  the  asymptotic  form  of  the  components  of 
the  .Hertsian  vector  and  subsequently  the  field  components  in  the  radiation 


CHAPfm  5 

FllLOS  AND  POm  FLOW  IN  THE  AIR  FOR  THE  OlFOLS  IN  MAGNETOPUSMA 

1q  the  preceding  chapter  ve  found  tlie  leading  ten^  of  the  a^ptotie 
representation  for  the  fundajBente^  integMls  U,  »  Nj  *  ^3  *  ^  addition 

ve  hare  dei^nstrated  that  the  aeo^ptotic  foxDs  were  differentiable  to  the  de« 
sired  order  and  as  such  they  could  be  used  to  find  first«  asQuivtotio  fomi  of 
the  components  of  the  Hertsian  rector  and  secondi  the  asgrmptotlo  foxns  of  the 
field  components  themselres* 

In  idxat  follows  we  shall  use  these  facts  to  arrlre  at  the  results  of 
this  diapter. 


5.1  THE  HiRlZlAN  VECTOR  AND  THE  FIISJS 


In  this  section  ve  first  find  the  aagrmptotie  forms  for  the  components  of 
the  Hertzian  rector  in  siRkerical  coordinatee.  Following  that,  we  express  the 
spherical  components  of  the  fields  in  terms  of  the  components  of  the  Hertzian 
vector  and  finally  find  the  explicit  forms  of  the  various  field  components  in 
SpheriofLl  coordinates. 

5.1  a  The  Hertjd^  vector^ApuIying  the  expressions  for  the  partial 
derivativaa  from  the  results  of  (4.76},  we  reerite  (3.274)  and  (3.27b) 


*«;cosV*vn.s^  cQse^;^1gJic  ^ 
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and 


U.  -^KptStrKi^StnaUj 

^  ^i;o.Y^n,*f --cos^e  ^J  ^Lk^(cos^sU*e  ^  j3^g£^;2^SS^)uJi>  • 
eos0  ■♦  /tx*- st>i.*0  *  V  '  Cos  &  /  *JJ 

Nov  Inoox^rating  the  z^aults  of  (4«62)»  (4«70)  aM  (4*74}  ve  obtain  the 
deali^  Insult 


(S.lb) 


likaK 


2lfisu  fA^l*\*c&s&  *  st«^*e 


li' 


I  iiiCKSmdl 


f. 


l*X*£05© 


^  - jiMjatarr”  ^  tltah^^Sc  l  I 

2  Sbn^v'n*«  St.ix*©'(Cd«9  Sen.*©'  ZSkn^  J  J 


and 


fc»x^ 


"T  m  COSLe  r  /n^Si. 

*®  2 


tK -  f — Sit^<seMn^  ©tKi5»aw^»cg^._ 

StVi  9  Scn5p!?crs??^[  nxoael/^^!.  £(c^0*  ) 


klc«r 


(5^2a) 


(5.2b) 


5.1  b  The  eleotrlc  field  in  apberioal  ooordlnateo^ln  Section  2.2  ve 
noted  that  the  eleetrlo  field  vas  deilvable  f^<ni  the  Hertidan  veotpr  a© 
followii 

Ip  ^  l0i)|*p^¥lfo  . 

Then  eadi  Curteelan  oonponent  of  the  eleotrio  field  is  giv^  hgr 

N^p  “ 

®yo  ®  “  ^^©)  (5*3) 

Njjp  “  *  • 

fo  find  the  ponponants  of  the  eleeUlo  field  in  spherioal  opprdinates  ve  use 
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(3.49)  and  get 

ico^o|(eos  Q  sin<f>3ji  +  sin  9  -  eos  9(sin<f 9^  aobi^piTig |  (5.3a) 

Eyo=  {  eosifdjX^-^  (cos<if9*+  eincf  3^  )%o  ^  (5.3b) 

Eyo-  *  |^(obs  9  9^  -  sin<f  sin  6  9^)li„  +  sin  9(sinY^-eos<p3pilio^  •(5.3e) 

le  now  use  the  results  of  equation  (4.76)  for  the  partial  deriTatlves  and 
Obtain  for  the  components  of  in  terms  of  as  follows: 

®eo  =  •  (5.4a) 

E^e  ^  -  k^c<Jfu„(eda<fa3a  B^o  *  sin  91^^^)  (5.4b) 

S.o  -  0  • 

We  note  that  to  the  first  order  in  r  the  electric  field  is  purely  transverse 
as  it  should  be  at  a  far  distance  from  the  source  (7>  p«  67). 

5.1  c  The  magnetic  field  ia  spherical  coordinate».^ln  Section  2.2  we 
noted  that  the  magnetic  field  in  the  air  was  derivable  from  the  Hertzian 


vector  as  follows: 

H,  =  klW,  +  v(^’ ir,  )  .  (5.5) 

Then  the  Cartesian  components  of  the  vector  are 

Bxo  -  isoTTko  +  a,(  +  a,T,^)  (5.6a) 

Hyo  -  a,(a*Tr,o  a,Tr«)  (5.6b) 

+  azTTzp)  .  (5.6c) 


To  find  the  components  of  the  magnetio  field  in  sidierioal  coordinates  we  use 
(3.49)  and  get 

®ep  *  {cssyoos  9(k^  +  4f)  ♦  ^(sin<fCOs  99^  -  sin  9  az  )  jlT*© 

(5.7( 

+  ^OOS  9  9,(008^^^  sin<f  ay)  »  sin  9  (k^  t  2^)  jH^p 
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®  i  *  ainsf  (k^  +  ^  )}T^Ke 

♦  ^g{Qoaef>hy»  8ln<^a^)TQg 


(5.7b) 


fi^e  *  9  (ko  +  3*  )  +  slQif  Bia  ♦  ooB  ®  ^  3^  l^xo 

(5.7c) 

+  |cos  ®(k^  ^)  +  ®  Bln^sin  0  9y)  jlT^o 

Nov  vs  uBe  the  fSBults  of  equation  (4>76)  for  the  ^tlal  deritratiireB  and 
obtain  for  the  eomponents  of  ^  la  terms  of  aB  fOlioVB: 


®  J'o  [  OOBy.  008  *  Bin  ®T*6  ] 

s  «  Bin5(sTK<s 

H,«=  0 


(5.8a) 

(5.8b) 

(5.8e) 


Ue  note  in  passing  that  to  the  first  Order  the  na|pietio  field  Is  porely 


transverse  as  it  should  be. 


Ve  note  the  foUoving  useful  relatloasbipB  Bnbng  transverse  electric  and 


■agnetio  fieldst 


(5.9) 


5.1  d  ^e  Taxation  field^We  shall  npv  vrlte  the  erpUolt  foni»  of  the 
field  oomppnentB.  Qsing  the  results  of  (5.4)  as  veil  as  those  of  (5.2a)  end 
(5.2b)  f  ve  obtain  for  the  coaponents  of  the  radiation  field 


>o 


2tf  (n»cos  9  > 


i<-d^>siln9 


.  co«Veo59 


(5.1Qi) 
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aad 


ZTr(n*cose*y'n*^sin^0  )  (  '  / 


(.>7Sl.»1  6^ 


case  y  iSui*^ 


'6  L  «.*co5  ©  Z(eos3  * 


(5.10b) 


SB  tk.k  CyL\eos& 


2(Cd6@ 


V6  Qdte  that  la  general  the  field  patterns  id^H  not  be  synmetrlo  in  the  0 «  as 
well  as  la  the  cf  •planes* 

Ue  note  that  sons  of  the  ooBpoaeats  of  the  above  field  expressions  be^ 
ooM  Infinite  idien  sin  0  s  n*  This  is  a  direct  consequenoe  of  the  approxlna* 
tlon  we  have  introduced  earlier  in  equation  (3.7)  which,  as  we  stated,  was  not 
valid  for  s  s  0  corresponding  to  sin  ®  =  n  in  the  above  results.  The  results 
for  the  radiation  field  are  valid,  however,  everywhere  SKoept  in  the  neighbor¬ 
hood  of  the  point  Sin  0*0  and,  of  course,  within  the  limitations  of  the  ap- 
proxlBation  itself.  This  stateisent  we  shall  prove  fomaJLly  in  the  latter  part 
of  this  work  when  we  discuss  the  solutions  to  the  line  sources  in  Parts  II  and 
Ill. 


As  a  partial  check  on  the  above  results  we  shali  go  to  the  lindt  as  n-=-*l, 
¥e  should,  in  this  case,  obtain  the  weU^lBiown  solution  for  an  isolated  nagne^ 
tlo  dipole.  To  this  end  we  note  laraedlately  frea  (5.Aa)  and  (5.4b}  that  when 


h  8  0  then 


¥  *  ^ 
0 


r 


(5.11) 


vtaiA  is  in  agreenent  with  the  well-known  result  (Z5,  p.  184). 
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5.2  TiiE  POWER  FLOW 


In  the  preceding  section  we  found  the  electric  and  magnetic  field  co:st- 
ponents  in  the  air.  The  form  of  these  comoonents  rev’eals  certain  additions 
that  are  due  to  the  anisotropy  of  the  plasria.  In  this  section  v;e  shall  deter* 
nine  the  effect  of  these  additions  on  the  direction  and  magnitude  of  the 
energy  carried  by  the  electroma^etic  wave. 

The  tlr.e*a7eraged  energy  density  carried  by  an  electromagnetic  wave  is 
given  by  the  Footing  vector 


?=^Re[Ixl»}  (5.12) 

which  in  our  ease  vdll  have  only  a  single  component  in  the  radial  direction 
given  by 

S^o  =  i  Re  I  (5*13) 

which  in  view  of  the  relationships  in  equation  (5.9)  can  be  written 

Sfo*  {  lEeJ"  +  lE^J"  ).  (5.14) 

Row  substituting  for  Eeo  end  from  (5.10a)  and  (5.10b)  respectively,  we 
obtain 


ytiere 


ko  m^  cos*  ® 


n^cosOv  /n*<»  sin*0p 


oosV  f'*’ I 


(5.15) 


+  goe^y  °QS  9 


(Q,<p)an^8iny»»  jy  sin  0 

Ln*cosO+  /n**  8in*Q'  2  /n^*'Si&( cos^  /^^35) 
+  4k,hooaV  j  (5.16a) 


and 
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J'-P  )  =  /n*-  ali^O  (cos  fl  /n^-  aLn^d  sln^O) 

_  (5.16b) 

+  j-yg1-bfl8j4<g  f  etM-  VaSslJ^^  ain*'e  _  1  -f  Ik^hCfl^eoo^  \/drj!jit*'6)  /j  ^ 

[  n^oosO^  /nr‘-Sla^'®'  2(oo8  6  +  /n*-  sin*®  )  J 


Ibe  power  pattern  repreBeatad  by  equation  (5*15)  ms  conpated  and  plotted 
in  Figures  5«2  -  5*4  Tor  the  dipole  situated  at  the  lower  edge  of  the  lonn^ 
sphere.  For  the  purpose  of  this  eocaaple  the  lower  edge  of  the  Ionosphere  was 
assumed  to  be  homogeneous  and  sharply  bounded^  having  an  eleetron  density 
H  B  750  eleetrons  per  cubic  oentlneter  and  the  earth’s  magnetic  field  .4 
gauss.  Since  the  eyclotron  resonance  associated  with  the  assumed  earth's  maga 
netlc  field  is  1.12  x  10^  oyeles  per  second^  to  stay  within  the  Hadts  of  the 
traUdlty  of  the  high  frequent^  appro^Lnationi  the  lomst  frequenqr  considered 
was  3  X  10^  oyeles  per  second  giiring  s  ,373  vtneh  was  assmaed  to  be 
adequate. 

Ramination  of  the  Figures  5.2  >  5*4  reveals  that  the  correction  due  to 
the  earth*  s  magnetic  field  manifests  itself  most  strongly  in  the  plane  s  0 
and  <f  B  'n'/4  where  is  azimuthal  angle  measured  from  the  au^s  of  the  di¬ 
pole.  Even  here  this  oorreotion  is  not»  however,  very  large  at  these  frequen¬ 
cies.  It  O8U1  be  seen  ^lat  for  moderate  polar  angles  this  oorreotion  is  oo9« 
pletely  negligible  and  it  begins  to  be  noticeable  for  6  >  60°.  Furthermoref 
this  correction  becomes  infinite  at  the  critical  angle  ®  idiere  oinr  high  fre¬ 
quency  apprcccimatlon  doea  not  hold. 


5.3  CLOSORS 

In  this  chapter  we  have  ocapleted  the  approximate  high  frequency  solution 
to  the  problen  of  a  borlaontiLl  magnetic  dipole  in  magnetoplaaDa.  In  partioulsii 
wi  forad  all  of  the  radiation  fl^d  eomponente  in  the  air  and  were  able  to 
separate  defiaitely  the  opntribations  of  the  pUima's  anisotropy. 


5.2  hmt  Fitteni  in  iir  of  a  Borisontal  Magnatlo  Olppla  in  Ibgnsto^ 
WMm  (r*  0»  ■  •  750  ^aatvona  par  oiMe  oentiaater.  h  =  1000  aatara, 
and  .ifuaa,  f  «  3  x  lO^adaf  par  aaaond)  ^ 


mm 


A 


MMl 
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the  resultd  of  this  analysis  were  applied  to  the  problem  of  a  horizontal 
■agoetle  dipole  situated  in  the  lower  edge  of  the  lonosjphere.  The  power  pat¬ 
terns  in  the  air  were  obV  Ined  for  frequencies  above  3  megacycles.  It  appears 
that  in  these  frequenciea»  the  oorroction  due  to  the  earth's  magnetic  field 
maMfests  itself  most  strongly  in  the  polar  plane  tlnreugh  the  axis  of  the  di¬ 
pole  and  T*74  from  it*  Moreover,  the  patterns  are  symmetrie  and  the  corTeo* 
tion  due  to  the  earth's  magnetic  field  appears  to  be  significant  only  for 
large  polar  anglesi  i*e.*  in  the  regions  close  to  the  interface* 


CHAPTER  6 


PORMtnUTlOH  OP  THE  PROBLEM  POR  A  MAGNETIC  DIPOLE  SOURCE  IN  BE  AIR 

la  this  chafer  ve  shall  be  eoaoenied  with  the  finding  of  appompii^ate 
Integral  representations  for  the  Cartesian  eomponeats  of  the  field  vectors 
for  the  aagnetoplassia^  and  aiz^half  spaces  with  a  source  in  air.  In  manjf  re« 
spects  the  present  boundary  iralue  problem  is  similar  to  the  one  in  Chapter  2. 
Vie  shall,  therefore,  avail  ourselves  of  the  many  results  of  that  Chapter  that 
are  applicable  to  the  present  problem. 

6.1  statemert  op  the  problem 


The  geometry  of  the  problem  is  shown  in  Pig.  6.1.  The  horlsontal  plane, 
8^0,  coincides  with  the  interface  between  the  anisotropic  homogeneous  plaaaa 
and  air.  Again,  we  shall  call  the  plasma  medium  (1)  and  the  air  medium  (0), 
and  assume  that  both  media  have  the  same  magnetic  inductive  capacity  of  free 
space,  .  The  steady  magnetic  field  Hoc  shall  be  oriented  as  before  in  the 
positive  direction. 


6,2  POHDAI^TIU.  EQHATXOBS 

The  definition  of  the  present  boundary  value  problem  Implies  splution  to 
Maxwell's  equations  subject  to  the  usual  bpus4ary  conditions  at  the  interface 
and  proper  behavior  at  infLaity.  As  before,  ve  shall  work  with  the  Hertsian 
vector  of  the  magnetic  type  in  the  air  regips  and  with  the  actual  field 


figsn  6*1  ^  QfMUx  of  FnUIiii  of  a  Higiiotlo  Dlpolo  iJB  Air 
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eoii^nenta  in  the  piaona. 


6.2  a  The  field  eouationsi^in  Section  2.2  it  via  ahovn  that  the  aource 
of  eleetroaagnetic  uatrea  way  be  reganied  as  a  inagnetio  dipole  alngulafity 
iMch  in  the  pz^sent  caae  wili  be  located  at  a  point  (0>  0>  h). 

then  the  Hertiian  sector  in  the  ai^  will  aatiaf^  the  inhonogeneoua  ireetor 
¥aiw  equation 

IS  -  -^|i^  £  (it)  £ (y)  <S  (a  *  h)  (6.1) 

vhtre 

Eo  ®  (6.2a) 

H*  «  k‘fo  +  f  (  0-if*  )  •  (6.2b) 

In  the  plasma  the  magnetic  field  must  satisfy  the  homogeneous  ^atem  of 
equationa  of  (2.26)  %diioh  ve  revrite  heza  tot  oontiaaienoe 


Xk* 


2x) 


H., 

xkl 

(6.3) 


and  the  eleotzlo  field  is  given  in  tezna  of  the  ugnetlo  field  as  follovsi 


S-. 

0  fK  ~f< 

E  yi 

"<.tO€.X 

Hy, 

H*. 

m.  « 

!»  « 

9  S 

6.3  FOORim  lNTBQRiy;i  REFRiSPT4TI0N  IN  GINTKIAN  COOBPlNAtES 


As  befpi^y  ve  shall  emidqy  here  the  triple  Fourier  tranafoni  pair  to 
iiaqillfy  the  fozoulation  of  the  present  problaii. 
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6«3  a  Ihe  oartle^^  eorra^panding  to  the  aouree — Applying  the 

Results  of  Seetien  2^2  to  the  wve  equation  of  (6.1)  ve  readily  Obtain 


“  (  ’0k  ) 

tdiere  ^  -  o< V  atl  as  before*  lie  inreirt  idth  respect  to  the  o/^  transform 
tradable  to  obtain 


*  tt 

4ir 


g  fcSi,|z-Kl 

si 


(6.6) 


Sow  inrertlng  vdth  reipeet  to  the  oi,  and  or^  transform  tariableSf  we  obtain 
the  well-knoua  result  (2^  p.  22) 


--e-f? 

8ir*uj^,  JJ 


-.Ofxy 

... 


dot,  dotj 


(6.7) 


6.3  b  The  oomolementary  integral  in  the  air-^te  saUsfy  the  boundary 
conditions  at  the  interface  ve  shall  need  appropriate  ocasplenentary  solutions 
of  the  differential  equation  (6.1).  these  can  be  written  at  once  as  follovst 

and 

Tf^u  [f  (6.8b) 

9^9 

is  all  of  the  abore  integrals  we  shall  require  ^t  Im  {so}  >  0  as  before. 

^.3  o  The  fields  in  the  nagnetoplaBiaa^The  fields  in  the  sagnetoplasma 
nust  satisfy  the  system  of  the  differential  equations  in  (6.3) •  Their  forms 
will  be  identical  to  ^se  obtained  in  the  previous  problem  in  equations 
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(2.64a) f  (2.64b,  and  (2.64o).  fhua,  v&  vrlte 


"•«'  "  8tr*cot* 


u  g  jgimka,  a  1 1 J 


<a<«.#«^^yxLfaigi  «feK)do(. 


(6*^) 


t  «  g>  f  ff  y^g^t'a.fc 

8t*cujA.  JJl  ^.a> 

^  4kX^|4o<_*  A.e"‘'**^  J ■*■  *»y a<.  dot*. 

^.Ca»>  ^ 


(6.9e) 


vhere,  as  before 


-  • 
T 


A2i^  .S2«j  *\/[(»t>4-  4fK^|<6.10O 


(6.10b} 


(6.10e) 


6.4  in  BODNOaRl  CONSITIOBS 

The  boundaiy  oondltlone  to  be  eatlsfled  In  the  present  problsBi  are 
identlail  to  tboae  of  Cbapter  2.  HovsTer,  due  to  ^e  faot  that  the  souroe  of 
the  eleetrpiiiflnetio  viTes  is  npv  situated  in  the  air  rather  than  in  the 
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■ignetopiaamay  they  oan  be  restated  In  more  eontrenient  form. 


6.4  a  Statement  of  the  boundary,  eondltlona-^The  boundi^  eondltlons  to 
be  satisfied  by  the  Cartesian  eomponents  of  the  field  vectors  require  conti¬ 
nuity  of  the  tebogential  oompcnents  of  the  electric  and  magnetic  fields  at  the 
interface  s  ^  o.  this  was  foxnally  stated  in  equation  (2.70) .  For  the  pur* 
pose  of  the  present  problem  we  put  these  equations  in  the  following  form: 


a  y  V 


V 


ki  hfzo  *  -  i  Hy, 

Ha.  ♦ 


(6.11a) 
(6.11b) 
(6.1 1C) 
(6.11d) 


lAiere  we  used  v  to  denote  the  field  expressions  with  the  integral  signs 
removed. 


Carzying  out  indicated  mathematical  operations^  one  obtains  fotir  alge* 
brale  equations  with  four  unknowns  ^ch  we  shall  write  in  a  eynbolie  form  as 
followst 


•  ^ 

— 

-V  — 

Om  0  ^.5  <»14 

B, 

b. 

® 

Ha 

8 

Q 

0  9  32  ®33  ^34 

A, 

0 

1 - 

t® 

O 

s' 

o 

t 

1 _ 

1: 

N 

_ 1 

(6.12) 


lAiere 


0  B  .  0(a. » CK XCkls.QC* 
'»  ckx^Cs.) 


(6.13a) 

(6.13b) 

(6.1|e) 
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e 

ex, 

Q  3  a^I ^  ■  CsQ  -^Ckl^  -  Ck  %  k»^j<>*a(.*),Si 

24  «i<i».?*V) 


'3t 


^3  Qtt  #,&J 

®S4.  oii#.C&0 


^1  “ 


1 

<s.ft.3/5i  —  C  K  0(t)  ^  A.yx  ^^tKQt«.Ol^x  ^  t's,  ^  ^-fc  K  X_4‘J<i‘‘At3 

So«a?X<#,C«-) 


.  cy.26C-^)[-j5.  ^66.)4-4K^<:k>0(t] 

<X#iC«) 


(6.13d) 


(6.1 3e) 


(6.13f) 


(6.1 3g) 


(6.1 3h) 


(6.131) 
(6.1 3j) 
(6.13k) 


’44 


<0(x  «.yx^^cK«,^)r«t^Cg.)4.c  icX^k*:s«3 

s*«*{:X#.Cs*) 


tl^X^kSofel 


(6.131) 


(6.13b) 

(6.13b) 


6.5  HIGH  FREQUMCX  APPROUKilllON 

The  detexninatlon  of  the  imknowns  B, ,  A, ,  and  fkom  the  e/atem  pf 
eiinations  (6.12)  is  a  straight^forvai^  bat  tedlpus  process.  The  results  would 
neoesssrlly  be  very  lengthy  and  probably  not  too  useful.  Iberefore»  as  In  the 
case  of  the  dipole  in  the  nagnetoplaaii  ve  shall  forego  their  evaluation  in 
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this  form.  Insteadf  as  before^  we  introduce  the  high  frequent  approxiiutleii 
(3*7)  which  ve  rewrite  here  for  eoavenienee 

B  *  0  (6,14) 

'  2  8 

Where  a  «  (Icf  -  af-  6^]^ae  before. 


^1?  t  ga’^laation  of  hi^  frequencgr 

approximBtion  we  find  the  ep^xiiaate  value  of  the  matrix  ooeffioiente  a^-^- 
as  foUovat 


o„  =  1 


’13 


»|4- 


»  n,  /  (^CWxk,  *40<.s)^3SLafA*)-^KQC.*  [  X 

^  I,  (  Cia;».k«i>t‘e<.s)^s^/)j!i*K«.*tC'e*-y^^^  ’t-oc.ottl  ? 


«az  * 


(6.15a) 
(6.15b) 
(6.15c) 
(6.1 5d) 


»a» 


^  *)>k* ot.or.^nVOlO^-^at^)  4 K  k.w, )Cy.<>rx* Ce^o^*) Vza  1 

(6.15e) 


CW,Oft(S^+«^)* 


Oaa** 


c„- 


O(.o<, 

i  k, 

bit.ofx 


O41  *  1 


c,»5«k, 


*44 


(6.1 5f) 

(6.1 5g) 
(6.1 5h) 

(6.151) 

(6.15J) 


■  j  (Ma) 
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Ptu^hennore,  ve  eu  reduce  tbe  system  of  equations  in  (6.12)  to  a  simpler 
e/stem  as  foUeuS: 


*2b, 


Whefe 


^43  *  ®ii 


n*8oO(4 


n^s«o(x 


and  a„  $  a,^y  aj,*  and  a^xare  ^ven  by  (3*11  )• 
Using  the  above  results  ve  obtain  finally 


(6.16) 


(6.17a) 

(6.17b) 

(6.17c) 


(6.l7d) 


-  i  *  sISfe  * 

B,  s  -/  2(»a*-»i)oc.  2y..Q(i,_  _  <x,Css’»-fyC/  ,  .  ] 

*  (<;sVs«Xs-f  rt*s#)  <LsCs-^wi^)  «iS(s-KS.)CS'^«.*a*1  J  i 


(6.1$a) 


(6.10b) 


veetoi>~Subatituting  the  results  of  (6.18a)  and  (6.10b)  into  (6.8a) 
and  (6.8b)  respectively f  one  obtains 


2yt5g^a*<a*K) 

~  s  *  ri^So 


(6.19a) 


f  rfgCSjz^M  _  e«S.(«.-KJ  2n5g^»*<a*K) 

*  8Tr*cu^.  JJ I  S*  ^So 

^ I* w f  j  ■  4  ■»<. _  ^ 

o<x§Ca-*'».Xs*ri^s.)J®  P  waort 
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_  nx 

ill 

f  2ai0m_  _ 1 

,(0«*-S*)(S+rt*Sa)  1 

.S(S-*>lT?S.)*  «,S(S-hS^<S-f  rt?4o)J 

,ecCct,K  ♦  C^y  +  (6.19b) 


U«  note  that  vhen  h  ^  0,  (6.19e)  and  (6.19b)  an  esaotly  the  aane  as  0,i2a) 
and  (3.^)  nspectlTely  as  they  should  be. 

Using  the  matlonshlpe  of  equation  (3*23b)  ve  noast  the  abon  integrals 


in  oore  oontrenient  forms 


■tie 


St*Ui(Xt 


-  g  X  ^  ot^y)  ^ a*,  d a* 

J 

’7L^^/J  J  (s+  rt.*5.)*  (6.20a) 

^  *JJ  oiiS(s-^sJ(s-*‘»a.*^s*)  J 


and 


bGvi  X 


-7M 


doc.da^ 


^  «  n1  ■  -  dot,  da  J, 


’►t 


4«« 

^  1  f 

•  ‘Jj 

(6.20b) 


6.5  a  Dsfinitinn  nf  f?^«YiiiBntAi  iotegrals«»’»To  facilitate  the  eraluation 
of  the  above  definite  integrals  «e  8hall>  as  in  Section  3.3f  define  certain 
fundamental  integrals  from  whiob  all  others  couJUi  be  derived  by  differentiiH 
tion.  TO  this  end  wa  dsfins 


c[ai;K-«-  ort.y-i' Sa^Z'V'K)] 


Y  «xy- 

J  • 


d<x,al«2 


(6.21a) 


(6.21b) 


rr  quLv.m.*  Olt.y  *Sa(t-*k)2 


S(s+rt.*4.)* 


Soi.S^z 


(6.21 o) 


V  S  - - - r-T-  ^  0(  ^  Q(^ 


(6.2ld) 


^  C  Coi ,  K -►  «iy  ti  S.^Z»  k)) 

■S— «■ — « — V  "  eioc.aax 
Oft*  (s  ‘*‘a«)(s  ♦  «.**•  ) 


(6.21  e) 


Bnjploying  the  above  definitions  ve  pan  irtwrita  the  Cartesian  components  of 
the  Hertsian  vector  as  follows: 

{t JJ[  ^  ]  e  - 


.  *  !>  ■  ■%  ^  ]  1 


(6.22a) 


the  trensfornations  of  equations  (3.2$)  and  (3«$9)f 


ve  rewrite  the  fUndanental  integrals  as  follnvst 


V.  e 


i-  fjihli: 

ko  i  B  +  a* 


K) 


(;\f)a  dA 


A  /|1  -  V  - 


f  ^Cs,(Z*h.) 

.“2“--- - -  Mi»(A9)AdX 

J  (8  4  8,)(8^n*So) 


^cs,(as-»K.) 

♦8o)(s 

_ _Cs«_.  ^ 

(*1)  aUi(2y+1)Y  -S - - I d> 

J  8(8  ^  Be)  (S^*0 


i»Md 


(6da3a) 

(6.23b) 

(6.23e) 

(6.23d) 

(6«23e) 


6.S  a  TrimdifQnBAttnn  to  utihayl^l  la  oonfiguratlon  and  traath- 

foaa  8Bae8>— In  the  oonflgtixatlon  Bpaee  we  tnuoBfom  88  8how  In  Kg.  6.1.  Vie 
note 


«  -  h  a  008  «o 

f  a  8ln  Oo 

r,»  /f*+  (e-  h)* 
and 


8  h  a  r  008  0 
^  a  r  sin  0 
r  a  /if*  +(•  ♦  h)^  , 


(6.24) 


(6.25) 


In  ^e  tnnefoni  ipnoe  ve  put  ae  before 

^  a  ke  8in(S  (6.26) 

Mid  88  a  oonaeqinenoe  obtain  the  folloving  reeultBi 


n 


(6.27a) 
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‘  I  {6.27b) 

n 

y  «  f  Q  ck.r<so.03^e;^p  (6.;wc) 

n 

f  a  r  — €  ck.reosr(9*©)d^  (6.27d) 

^  I  (eosib*  vVt*-  Sbn.*/s)(v^*£os^ /rt*-  su»-i3jS  ) 

*#  ^  .  A.../  .  ^^GOSAHiv^iCWe*^^  _ .>  V 

V«0  n 


6b6  CLOSURE 

In  the  foregoing  chapter  ve  first  rigorously  formulatcKl  the  problem  of  a 
horlsontal  magnetic  dipole  in  t^e  air  in  the  presence  of  a  magnetoplaaina  half- 
space*  FoUoidng  the  rigorous  formulation  we  then  applied  the  high  frequeniy 
approximation  as  in  Chapter  3*  This  enabled  us  to  find  the  explicit  form  of 


the  boundary  coefficients  and  the  approxiouite  inte^pal  representations  of  the 


Cartesian  components  of  the  Hertzian  rector. 


Purthemore,  as  a  partial  check  on  the  results  of  this  chapter  as  well  as 
those  of  Chapter  3,  ve  hare  shown  that  the  integral  representations  of  the 
Hertzian  rector  In  the  cases  of  the  source  in  the  magnetoplasma  and  In  the  air 
are  Identical  vdien  the  depth  of  the  source  burial  Is  serof  that  ls»  ^en  the 
Bouroe  Is  situated  right  on  the  boundaxy. 


CHAPTER  7 


FIELDS  AHO  POUm  FLOW  IN  AIR  FOR  THE  OlPOLi  IN  AIR 

In  the  preceding  chapter  ve  found  the  approxiinate  Integral  expreesiona 
for  the  Cartesian  oomponents  of  the  Hertzian  vector  in  the  air*  We  also  d»> 
fined  certain  ftmdaaental  integrals  vhieh,  vSien  avaluated>  can  be  used  to  find 
the  explicit  form  of  the  Hertsian  vector  and  the  field  ooB^nents  bgr  differem» 
tiation. 

In  this  chapter  ue  first  evaluate  approxiaately  these  fundamental  inte^ 
gralSf  find  the  components  of  the  Hertsian  vector^  and  finally*  the  oomponents 
of  the  radiation  field  and  the  Pointing  vector. 


7.1  EFALOATION  OF  THE  FONOiMTAL  INTEQSALS 

The  fundamental  Integrals  of  the  previous  chapter  are  very  similar  to 
those  of  Chapter  4*  the  only  essential  difference  being  the  absence  of  the 
factor  eoqp  (Ish),  The  reason  for  the  factor  ''h”  not  appearing  In  the  Inte^ 
grals  under  oonslderatlon  Is  the  fact  in  the  present  problem  the  corresponding 
factor  exp  (iSeb)  was  combined  with  exp  (iSes)  hF  naans  of  a  proper  shift  in 
the  coordinate  SFstem*  and  thus*  absorbed  in  the  radial  distance  r  or  re  • 

The  basic  method  that  has  fomiulated  in  ^pter  4  can  be  used  here  effeotivily 
and  the  evaliiation  of  the  present  integrals  can  be  accomplished  almpst  by 
inspection* 

first  we  consider  the  first  two  integrals  in  equaUon  (b«2Ra)*  Ihese 
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fofls  are  vell^knovn  (2,p.22)  aM  the  results  can  be  written  at  once^  Ue  ha\fe 

(7.1a) 


.  "  ■  . r---.,-r-r  d«.dS(l»  “  ITt 


II 
fl^ 


.  - dot,  d«j  =  «  ZTTl  — ^ 


{7.1b) 


Next,  we  foeua  our  attention  on  the  Integral  V,  In  equation  (6.27a)  and 
ooiqpare  it  with  the  integral  U,  In  equation  (3.44).  We  note  that 

V,  sU,  (h«0).  (7.2) 

IBwit  W  (4*61) 

itW,>r 


^  _  \ 

'  ko  I  n*oos  9  +  /n**.  sln^d  /  r 


(7.3a) 


Sli^Jjarly  oonparlng  In  (6.27b)  and  tl.  In  (3.44c)  and  (4.74) »  we  can  write 


* 


xz 


f  . . . . _,co8._6  . .  "I  ,el‘**'' 

L  (n*oos  0  +  iin^isi)  ^  J  r 


(7.3b) 


ihialogouaiy  with  Yj  and  Y4 


_  _  ^  f  _  ooa  «  '] 

kp  L  ^B^-sl^(n*coB0  +  /n*’-  sia*'0)^  J  r 

12  r _ 008  9 _ 1 

kp  I  (ope  0  +  /a*»ain*0)  (n^coa  ®  +  /n*-sin*9)  J  r 


(7.3c) 


(7.3d) 


Finally f  comparing  in  equation  (6.27e)  and  U2  in  equation  (3*44b),  we  ob^ 
tain  using  (4.70) 


kaj.  8in08in<f  ya 


008  9  ) 

i^»idjnF9(  ooi  0  ♦  (n^ooaO+  /n^^ello)  J  ^ 


(7.3e) 
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7.2  THE  HESIZIAN  VECIOE  AND  THE  FIELDS 


Haring  found  the  leading  tezins  of  the  a^ptotlo  eocpanelona  for  V,  throogli 
V5  ve  aball  now  use  these  results  to  vrlte  dovn  the  components  of  the  Hertslan 
reotor  In  explicit  form  and  later  the  oomponents  of  ^e  radiation  field. 


7.2  a _ The  Hertslan  vmotor«>Apolylnf  the  expressions  for  the  partial  de« 

zlratires  from  the  results  of  (4.76),  «e  rewrite  (6.22a)  and  (6.22b)  as 
foilowsi 

4T60M«  t  ra  r 


4T60^«  t  To  r 

r  It^kgSln  0  ^isin<fV^«  sin  0  cos  6  cos^  |  (7.4a) 

•  »3 

+(slnifsln*^©  +  00^6)  7s] I* 

Nov  Inoorporattng  the  results  of  (7.3)»  we  obtain  for  the  radiation  field 


T  j-  *mlt|»OQa<f  sin  6 


1  « 

4T4i)^o 

1'  r. 

r  1 

1*^  sin  < 

}  oos  0 

f  2  aln^ 

(n^ooa^ 

/n^oeijD  ) 

1  n^ooaOf  /n*«i 

2B^e08  9  •* 


r  n^oos  •  +  /n^-  slni‘0 


(7.5a) 


105 


7«2  b  The  widlatlinn  H^fjl  tin  shedl  now  %alte  the  expUeit  forms  of  the 
field  ooi^hsnts.  Using  the  restalts  of  (5*4)  as  veil  as  those  of  (7.5)>  ve 
obtain 

a  i  ain»  / s*^ **«»*;  ^  e^^\^ _ n^goe  0  siny_ _  e^***;*“ 

2T  c  2  \  r  r«  /  n^oosO^f  /n^-  slii^O  r 


I  7  sin  6  w  s  0  f _ _ 

n^ooid  ■*■  /n^tein^g  L  n*oos  0  * 


(7.6a) 


ops  0  cos 


and 

‘f’ 


"lluiwosy  epsO 
4Tr 


2  /n’^Aikn^OC  OPS  0  ‘f  ^ 


^  Atkor 


_ 1  e*’***^  ) 

sin^Q  )  J  ) 


4  ^  /n**’»ia^O(sin^O»eosO_v^»Sl^  )  ,  lilllL 

'  '  n* OPS  0  V  r 


i9siny  sin  0  f  2(oin^0  •»  opsQ  ^n^«i»  oin*0  ) 

‘-sin^O  (n*oo8  9  ♦  /n*-sio^9)  i  n^oes  9  +  /a*-  sin^O 

r 

_1 _ 1  e*^^^  7 

+  /a^-  sln^O  J  r  i 


(7.6b) 


OPS  0  8in*9  J  r 

lb*  aagnetlo  field  opagipnents  aura  related  to  tbs  eleotvie  field  ponponents  as 
la  (5.9)f  is. 


>«<p 


<fO 


H 


(7.7) 


eo 


vbors  Sp  is  the  fkes  ipaes  impedsnes. 
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7.3  THE  POWER  FLOW 


The  radiation  field  oonponents  fotmd  in  the  previous  seotlon  can  now  be 
used  to  obtain  the  tisu^veraged  P^pting  vector  given  by  (5.14) «  Thus*  we 
obtain 

3^^  -  1  ^  1^  ■*■  I  0^(6^cf)l^  j  (7.8) 


Soo  s 

fo 

Sir*  Zo 

Where 

^•0  ” 

/a 

2  ( 

(  JP  To  / 


008 


n^oos  Oi-/i 


^asssBsasr  ^  B^Slnc^ 
/n*«8in*0  2  ^ 


♦  -  t  — - 1  i 

In^oosQV  ^x^-sin^Q  2  ^?-slj?0(co8d+ /?«^aii?0  )  J  r  3 


(7.9a) 


n*  oosO  V  /n*-3in^0 


I  ya*»8ia*0(8ia*0+ coed  /n®-»8in*9) 


ly  alny 


sin  0 


/n’’-  sin*« 


sln^OteosC  /n*^-sli^ 

^  _ 1... - ] 

1  e‘^'‘-»“ 

n^cos  0  V  /n*-sla*© 

2(co8  8+/n*^-  sin^O  j  J 

ST” 

Again  ve  note  that  when  b  =  0  the  Pointing  vector  In  this  case  Is  Identical  to 
that  In  equation  (5.15)  and -corresponding  to  the  case  when  the  dipole  Is  situ-* 
ated  In  the  Bagnetoplaana. 


7.4  CLOSURE 

In  this  chapter  vs  haTe  oonpleted  the  approximate  high  frequent  solution 
to  the  lOoblsB  of  a  horlsontal  magnetic  dipole  In  air  In  the  presence  of  a 
xegnetoplatsaa  half<!>8paee.  In  particular}  ve  found  aH  of  the  radiation  field 
oonponents  In  the  air  and  were  able  to  separate  definlt<^  the  oontrlbutlons 
of  the  plasna's  anisotropy. 

As  a  parUal  eibsek  on  the  results  of  this  chapter}  as  well  as  ^so  of 
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Qlapter  5|  ve  noted  that  the  coxTeepondlng  results  were  Identloal  when  the 
depth  of  the  source's  burial  in  each  ease  was  aUoved  to  approach  seroy  l«e»y 
Idien  the  source  was  situated  at  the  interface^ 

Although  no  nuaerloal  eKample  was  considered  in  this  probleniy  the  results 
can  be  expected  to  be  substantially  the  sane  as  those  of  Chapter  5  and  thera^ 
foray  the  conclusions  drawn  in  that  chapter  also  apply  here. 


PAM  ii 

FIEU)  OF  H.£0tRIC  CUBMNT  LINE  SOURCES  IN  MAQNEIOPLASMA 


UIIH  A  SEPARATION  BOONQARX 


CHAPTIIR  8 


BIQOBOUS  FORMULATION  OF  TH£  PROBLEM  OF  AN  ELECTRIC  CURRENT  LINE  SOURCE  VOiEN 
THE  STEAES  MAGNETIC  FIELD  IS  PERPENDlClfiiAR  TO  IT 

In  this  ohafitor  ve  shall  be  eoneenied  with  the  finding  of  appropriate 
Integral  reprosentatlons  for  the  Cartesian  oonponants  of  the  field  veotors 
for  the  magnetoplasina-  and  alr»half»  spaces.  This  vlll  entail  revlev  of  the 
fOndaaental  field  equations  and  the  definition  of  ths  source  of  eleotro* 
oagnetlo  vaves. 


8.1  STATIMBiT  OF  THE  FROBLBf 

The  geometry  of  the  problem  Is  shown  In  Fig.  8.1.  The  horizontal  plane 
■  s  0  oolncldee  with  the  Interface  between  the  anisotropic  homogenous  plasma 
rad  air.  Again  for  oonveidLenoe,  ve  shall  call  the  plamna  medium  (1)  and  the 
air  medium  (0).  Aa  before»  ve  assume  that  both  madia  haro  the  same  magnetic 
Inductlro  capaol^  of  froe  spaecf  The  steady  magnetlo  field  Hp^  vlll»  In 
iHaSLB  case#  be  oriented  perpendicular  to  the  Una  source  vdiersas  the  plane 
fomed  by  the  line  source  rad  the  direction  of  the  steady  magnetic  field  Is 
parallel  to  the  Interface* 
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f^urt  $«1  -  Qttmtiy  of  th«  FxoUm  of  an  iQnotj^o  Guorwit  I4na  Soaroo  Umb 
tlii  StMdtf  Nugnallo  flald  is  Fsrpfndiealar 


Ill 


B,2  FUUDUmTAL  fiqMATIONS 


The  definition  of  the  present  boundazy  valiie  prohlem  implies  solution 
to  Maxwell's  equations  subject  to  the  usual  boundazy  conditions  at  the 
interface  and  proper  behavior  at  infinity.  The  use  of  the  auxUiazy  vector 
potentials  does  not  seem  to  siapliiy  the  prohlem  any.  Thus«  we  shall 
operate  with  the  field  components  direotiy. 

8.2  a  The  nature  of  the  source^for  the  purpose  of  this  problem  it  will 
be  assumed  that  the  source  of  uhe  electromagnetic  waves  consists  of  a  Very 
thin  straight  %dre  of  Infinine  extent  carzylng  an  alternating  current  • 

To  localiie  the  source  properly  we  write  for  the  electric  current  density 

J  s  I^(x)S'(s  +  h)  (8.1) 


where  S()is  the  Dirac  delta  function  (7|p.43). 


8.2  b  The  field  eQuatione-~We  assumed  previously  that  our  source  of 
electromagnetic  waves  in  this  problem  may  be  regarded  as  an  electric  current 
line  singularity  at  a  point  x  ^  0  and  s  >  •  h .  Then  for  the  megnetoplaame 
region  the  appropriate  form  of  thucwell's  equation  is 

^  X  E,  =  ico^,H| 

^  x^,  «  •  iwe^f  +  l^(x)^(s  ♦  b)  ly 

(8.2 

f  ,  a  0 

Ubarw  f  if  the  electric  etamrge  doneily  related  to  the  electric  current  by  the 
continuity  e^tion 


(8.3) 
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lb  obtain  tha  wva  aquation  we  perfozn  a  curl  openitlon  on  the  first  aqinatien 
of  (d«2)  and  than  aobstituta  for  ^  x  H  firoa  the  second  and  obtain 


-  ^  X  ^  X  f,  +  k“  ^  -  i&>^.l5(x)  i‘(s  +  h)  . 


(8.4) 


Oiing  the  ooaponants  of  the  tensor  as  in  equation  (2.8)  and  eaxxying  out  the 
neoesiarj  algebra  noting  that  one  obtains  a  iqrsten  of 

differential  equatioas 


-f  a/  0 

0 

0  ek*  f 

*y. 

B  -i^/^  S(x)S (sfh) 

1 

0 

whieh  we  shall  leaire  in  this  fona  for  the  tine  being.  Pron  the  first  equation 
of  (8.2)  we  note  that  the  aagnetie  field  is  expressible  in  terns  of  the  ele»> 
trie  field  as  foUovsi 


*  • 

••  «■ 

H*, 

Hy, 

H*. 

0  •bg  0 

3,  0 

»y. 

1.. 

iCJfA, 

0  2k  0 

in  the  alr-r^(ion  the  appropriate  foxn  of  the  Muniall's  equation  is 


(8.6) 


%,  «  ico^,X 

^x  K,  »  *  iu>d,To 
V-  S’,  Q 
**  9 


(8.7) 


and  the  Cartesian  oonponents  of  the  eleptrio  and  nsgoetio  fields 
seetor  nave  equation 


satisiy  the 

(8  .$ 
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8.3  FOURIER  INTEX^  RSPRESEiNIiLTlQt]  U  CARTESIAN  COORDINATES 


The  formulation  of  the  present  boundary  value  problem  can  be  simplified 
a  great  deal  hF  socpresslng  the  field  oomponents  In  magnetoplasEia  and  In  the 
air  In  terms  of  their  double  Fourier  integral  representation  In  Cartesian 
coordinates  In  the  transform  spaoa  as  well  as  In  the  configuration  space. 

Thus*  we  Introduce  a  double  Fourier  transform  pair  defined  by 

»  (<»‘.»ai)  *  ^  J[  ^  (***)  d  X  d  s  (3.9) 

and 


f  Cx»t) 


(«.,  0^3)0  ♦V)  ^ 


(8.10) 


In  tSiat  follows  we  shall  also  need  the  transforms  of  the  derivatives*  These 
can  be  obtained  Integrating  by  parts.  As  we  showed  In  Chapter  2  for  the 
triple  toransform,  the  vanishing  of  the  Integrated  part  at  the  upper  and  lower 
Halts  is  assured  provld^  the  radiation  condition  is  satisfied.  ThuSf  assuoH 
ing  that  the  radiation  oondltion  is  satisfied  we  can  establish  the  following 


correspondences: 


(8.11) 


8.3  a  to  the  sourcOi-^TO  transfoim 


the  inhoBogeneotts  systen  of  simultaneous  differential  equations  (8.5),  one 
multiplies  both  t^des  by  mcp{iiot,x  t  oc^s)^  and  jetegrates  with  respeo;,  to 
the  real  variables  z  and  s  betueen  «  ^  and  4^  cko  •  Tbe  ri^Vhand  side  ef 


($.5)  yields  at  once 


JJ^(x)S  (s  +  h)  e*^(^**'*’^**!?  d  z  d  s  ^  e‘^i^ 


(8.12) 
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and  the  left-hand  side  tx^nsfozms  aeoording  to  equations  (8.9)  ajid  (3.11). 


Thusi  one  obtains 


et}-  cxi 


7^*  « 


(8.13) 


^ch  can  be  liimediately  reduced  to 


O  oi,Qt_ 


^ere  now 


^0.?]%  4K*{k;«*  I  (8.15) 

Ue  compare  (S.13)  ulth  (2.41)  and  obserre  tliat  the  former  is  the  same  as  the 
latter  witb  oc^  -  0.  From  (8.7)  ve  can  vrite  each  transformed  component  of  the 


electric  field  as  foUovs: 


f fo)  3  .  rCakl-«.*)-  6e<^1  CctsW 

Sire  L  Cor3»-s»C«/.i^)  r 

f  ^(^kl^oc/’)  ]  ta,K, 

2T€  L  J  ^ 

-  AHt^M.1  f _ _ 1  c«,k 


(8.16a) 


(8.16b) 


(3.l6e) 


The  inversion  \cLth  respect  to  the  ot^  transform  variable  can  be  performed 
iomediately.  According  to  (2.53)  and  (2.54)  va  obtain 


-^j**”**^”  *^****  ^fek. -ocr)^  6S/^flg.l2*kl  1  #{.  . 

2>/^6  I  s.(a^.s^)  ^  j 
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**1 


5.(4.*- S»*)^  I 

f(P)  *  ^  isi£0*^  f  I 

*'  2H*x^d  2N/2ir  'J  s**s^  s«-Si*  J 


(8.17b) 

(8.17e) 


Finally,  ve  Invert  with  respect  to  the  transform  variable  and  obtain  the 
desired  representations 


®r  -  «:4.iz.-kU 

•i-*# 

«d 

Jg;,  (8.18b) 

f  (8., 8,) 


fa  the  preceding  secUon 

ve  obtained  the  "particiUnr  Integrals”  of  the  system  of  equations  (8.13)  vbloh 
represent  primary  excitation  due  to  the  source.  TO  satisfy  the  boundary  oon> 
dltions  of  the  problem,  ve  shall  need  an  appropqrlate  oompleaentary  solution 
of  the  homogeneous  system  of  (3.13). 

From  (8.14)  It  Is  clear  that  satisfies 

(«»'-.?  )(«.•»  »|  )f«’»0  (8.19) 

vhere  In  the  above  equation  ot^  Is  looked  upon  as  a  differential  operator  vitb 
respect  to  the  iHooordInate.  The  aolutions  to  (8.19)  can  be  vritten 
limdlatialy 


^*1 

"c,. 

Ciz’’ 

•yi 

* 

®a» 

Cgx 

a’-vsiz 

.1S>, 

a» 

_ li 

(8.9)) 
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^^ere  we  discarded  solutions  wlt^  positive  exponentials  since  we  can  have 
onl^  waves  going  away  from  the  Interface  upon  reflection,  the  coefficients 
C^j  can  be  found  by  a  method  Identical  to  the  one  of  Section  2.3.  We  obtain 


C,.  ^ 


C3,  » 


1  y  pt.  jp. 


otj^)  -  esf 

i?  ) 

<  («  k* 


Ci. 


Cl. 


and 


C.i  = 


i  ‘7  _^  l  J,z 


>22 


§  (  €  kf  -  af )  *  6  Si 


C32  =  .  Cji 

C(«kf  -o<?)*68^ 


(8.21 a) 

(d.2lb) 

(8.21  e) 

(8.2ld) 


For  convenience  In  what  follows  we  shall  nozinalize  the  Coefficients  0^-  and  put 


Ca.  « 


2\/2T€ 


d| 


(8.22a) 


Caz  * 


- 


2v^e 


(8.22b) 


Wov  we  Invert  with  respect  to  the  ot,  transform  variable  and  obtain  the  desired 

(8.23a) 
(8.23h) 
(8.23e) 
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5k.  4lf  ^ 


It  oan  be  reaMly  sboin  that  the  aboTe  field  coaponenta  satisfy  the  divefgeaoe 
equation  s  0. 

8.4  the  BOUNQARX  COMDXtlONS 

In  the  preceding  sections  ve  found  field  components  in  the  pHaaraa  and 
air-regions  that  are  solutions  to  Maxwell's  equations.  Moreorer,  in  solving 
Maxwell's  equations,  we  chose  such  solutions  for  field  representation,  that 
have  proper  behavior  at  infini^  by  requiring  that  the  ioaginaxy  part  of  the 
pertinent  exponents  be  non-negative.  In  addition,  the  field  components  con¬ 
tain  certain,  thus  far,  undetensined  coefficients  which  upon  imposition  of 
the  boundary  oonditions  v*  11  be  determined  and  thus  render  the  solntioa 


8.4  a  St«*^»«»««^  C»»*dition8— The  boundary  conditions  to 

be  satisfied  by  the  Cartesian  components  of  the  field  vectors  require  conti¬ 
nuity  of  the  tangential  oomppnents  of  the  eleotrio  and  na^etio  fields  at  the 
interfaee  i  s  o.  Ibis  inpilies  the  foUowingi 


Hko  *  K. 


Hyfl  “ 

^X6  “  ®X| 


Eyo  “  E/>  * 


(8.25) 


Ue  teweitn  the  aboi/e  equatione  In  texna  of  the  eleotrlo  field  on!lj|r| 


getting 


E  -  E 

"xi 


®yo  “  ®yi 
^*Ey^  “  8yi 


*“  ^®zd®  ^®Xi  **  • 


(t.26a) 

(3.26b) 

(8.26e) 

(8.26d) 


These  equations  oan  be  alm{di.fled  to  a  nor#  oonvenlent  font 


0  s  i  ^2  Ey,  +  s*  E^ 


(3.27a) 


^  V  ai  4  V  V 

0  « t  a. t  ♦  «,**, 

•o 


(3.27b) 


idiere  v  denotes  the  field  representations  uLtb  the  Integral  sign  rcsio/ed. 


8.1  b  Apniiaxtina  of  the  boundary  condltlon^Perfonaing  the  neoessary 
algebrale  operation  to  satisfy  equations  (3.27a)  and  (3.27b)  glTes  a  set  of 
sijedtaneeus  algebrale  equations  in  tvp  mtigiouns  as  follovst 


F^CiSi} 


s*^s.  A, 


PaC»d 


pi:s.h.  _  Sa.-Ca«  tfiasK.  (3.23) 


(■..»)  *  <(«  k*  •  a? )  ^  ^  B,%  . 


(8,29) 


Hi  iol?s  the  shove  gfstoi  of  equations  using  Gveser'i  rule.  He  find 


(8.30a) 


119 

.  ^  pK-C$!i)[CSiTSo)C^Sx  ->'4lS.)RitSi)  —  ^Sf»So)^S|-  ^5«)Pa^6«)l  ^  lSiK. 

*•  s.(s.*-St*)Kle 

“ 


and 


= 


'g*r  ,S>)CSj  -jh <S  JPaCSt)  - (S,-*-S^Sz ^c'5*k 
Ssi(5.*-^s^)Ne 


(8.30b) 


ligs  (a,  ♦  8^)(a^  +  ^a.)  Pe(a, )  *  (a^+  So)(a,+  s®)  ♦  (8*31) 


fo  find  B,  and  M  uae  {itTha)  and  (B.l^b)  and  the  above  resulta  and  obtain 


B,  =-26 


Pg  (a.)(aai-<an)e**‘*"  *  Pe(aa)(a.  t,^ao)e*^'^*^  ^ 

3 


(8.32a) 


and 


Bj  =  2  4:  a. 


a,)e^^>‘*  (a.  ♦  s,)e«-^**^ 


1 


(8.32b) 


8.5  FIEU)  CQHPONBNTS  IB  THE  AIR  IN  CIUNDRICAL  COORDINATES 


Having  found  rigorous  expressions  for  the  field  components  In  the 
Cartsslsn  coordinates  in  both  of  the  regions^  we  can  now  transform  to  ^llnK* 
drlcal  opordLnates  in  whl^  the  field  oomponents  will  have  soaetdiat  simpler 
fon. 

He  transfoni  to  cyllndrloiJL  coordinates  in  both  conflgurati^  and  trai^ 
form  spaces  using 

ot.  m  k^sin/S 

8  =  ^  COS  0  (8.33) 

X  =  f  sin  0 

As  a  oonsequmee  of  the  above^  we  observe  that 
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»,x  +  8„a  ®  00B(jS-  0) 

(8*34) 

ftothenaore, 

au^8ine2^ 

•3e 

(8.35) 

®  OOS  0 

•3. 

and 

(  )e  =  •"  (  sin  ©  +  (  )^ 

(  )y  =  (  )2  008  0  V  (  }^ 

Cos  0 

sin  0 

(8*36) 

Usl&g  the  above  x^lationship8>  the  eleetrle  field  doni^nente  oan  be  wltten 


c 

B^^ar  3^J cosiSco»(j&-e)e^'*»y*«»*(i»*»>d/l  (8*37b) 

n 

By^g*  ^ f co*^ SLn(j»-©)ei««^oc>»(h-«>d/5  (8.37o} 

n  * 

^ere  P.  is  the  integration  path  in  the  complex  ^  -plane  as  shown  in  Fig*  4*1  • 
In  the  above  equation  it  is  understood  that  s, ,  and  s^  are  all  functions 
of  15  aocording  to  the  first  equation  of  (8.33)* 

8.6  GLQSURi 

In  the  foregoing  chapter  we  formulated  rigorously  the  problem  of  an  ele^ 
trie  current  line  source  in  a  magnetoplasiia  with  a  separation  boundaxy*  In 
partioulaTf  we  found  the  integral  representation  for  the  field  oaip>n^ts  in 
the  Cartesian  and  cgrllndrioal  ooordULaates  in  the  air<”regioa  and  in  the 
teaian  coordinates  in  the  plnsma  region*  Onlike  in  an  ordinazy  isotrppie  irnrn 
all  of  the  eleetrle  and  nagne^o  fields  are  preswit  in  both  r^lons* 


CMPTER  9 

RiSULlS  fOK  THE  AIA-RBQION  WHEN  THE  StSAOX  MAGNETIC  FIELD  IS  NORMAL 

10  IHE  LINE  SOURCE 


Having  obtained  the  rigorous  fornal  solution  to  the  problem  of  an  eleo» 
trio  cuxrent  line  source  In  a  nagnetopilasBa  ulUi  a  separation  boundaxy  In  the 
fotm  of  oertaln  definite  Integrals^  we  ahadl  endeavor  in  this  chapter  to  V9» 
dttoe  these  Integrals  to  a  form  that  Is  suitable  for  numerioal  calculations. 

As  It  was  remarked  In  Chapter  4>  these  definite  Integrals  do  not  lend 
themselves  to  rigorous  evaluation  and  we  shall,  therefore,  use  the  method  of 
steepest  descents  to  extract  the  desired  Information  from  them.  Ihe  method 
to  be  used  here  vUl  be  somewhat  simpler  than  in  Chapter  4  ^nce  only  e  ein|^ 
integration  is  neoessaxy. 

9*1  SllklULARITlES  IN  IHE  3  -PLANE 

For  reasons  mentioned  in  Section  4*1 »  It  is  mandatozy  to  know  the 
singularities  of  the  integrands  before  their  apprpadjsate  evaluation.  TO  tfeis 
end  wa  shall  now  examine  the  slngularitiei  of  the  integruds  in  (8.36a), 
(8.36b),  and  (8.360). 

These  integrals  whan  tnmsfozmad  to  the  (8  -plme  take  on  the  fon 

cpsfi  (9,1i0 
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f5 

.gtkfco«C(9-»)^p 

*f’*~  ’^m'l'^  «•’«> 

uben 

o; ‘  4:) -  0 ^ )si*i?<S  ^  ^[(x-?)- (i-  |.)si:rt?/5]*+  4k*<&^  j*  (9.2a) 

Jtg  >  ('eiH-  <:os|8X^a^ ^cosA)Pef^)“  C«z*eos(»)C«^  ^  4:co5(a)P,(isi)  (9.2b) 

He  also  note  that  P^  (<r,)  and  P^  (s^)  can  be  written  In  the  form 

Pe(«.>)  =  -  [.^7*-«|"ft-«aiaa?]  (9.3) 

and  that  Pi  and  er^  can  then  be  obtained  from  (9.2c). 

9.1  a  The  location  of  the  tx>le3-.^Ihe  poles  of  the  Integrands  in  (9.1a), 
(9.1b),  and  (9.1c)  are  the  sexoes  of  the  denoadnator  iZe  •  Thns,  to  obtain  the 
location  of  the  poles  ve  most  sol^e 

(s;+  ops/i)(C5.+  Cqo8^)PE(er)^^+  cos/t  )(6r  +  ^  oos|i  )Pe  (^)  (9.4) 

Unfortunately  ^e  analytical  solution  to  the  above  equation  is  very  difficult 
to  tiad,  and  one  nnst  resort  to  graphical  means  to  accomplish  it. 


123 


9.1  b  The  branch  polnta—The  integrands  under  consideration  contain  the 
radicals  associated  with  cr,  and  er^  as  a  result  of  vhloh  the  points  /3  -  O01 
and  B  Og,  at  vhlch  er,  and  respectively  vanish  %iill  be  branch  points. 

At  each  point  ^  these  integrands  can  take  on  four  different  values  depending 
on  vhlch  sign  ve  choose  for  the  radicals.  It  vill  be  convenient  here  to  talk 
about  four  sheets  of  the  /S -plane  (fonned  by  a  four^sheeted  Riemazui  surface) 
on  which  each  integrand  is  Single-valued.  These  four  sheets  will  be  defined 
as  follovst 


Sheet  1 

Sheet  li 

Sheet  III 

Sheet  IV 

lajcr.i 

+ 

- 

+ 

m 

mieii 

♦ 

+ 

- 

- 

The  ^nvergenoS  of  the  integrals  is  assured  if  the  path  of  integration 
at  least  begins  and  ends  on  Sheet  1.  For  convenience  we  shall  introduce  the 
feUowlng  "outs*  in  the  ooni^ex  (i  -pianet 


0  .  (9.5) 

It  is  clear  that  the  equations  for  the  "outs"  met  satisfy 

^  0  ,  (9*6) 

hov  the  branch  points  are  located  at 

Oei  ^  ^  VO  alh  (  /a ¥7  )  (9f7a) 

»e^'  X  Vc  sin  (  »  7 )  ,  (9.7b) 


Ve  recall  that  /e +7  and  ^e  »“7  vn  respoctively  the  indices  of  refrae- 
tion  of  a  rii^ti-  and  left-hand  circularly  polarised  plane  wave  travelling  in 
the  direc^on  of  the  steady  aagnetio  field  in  the  plaaia.  Furth^Raorey  in 
tenn  of  a  losiless  plasiiibPin>M^oFs  (a ^7)  and  (*«'7 )  arc  givsn  by 
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6  +'7  “  1  -  •  (9.t) 

We  obsezve  from  the  above,  the  location  of  the  branch  polnta  Is 
strongly  dependent  on  the  value  of  plasma  and  cyclotron  frequencies  In  the 
plasma.  When  p^  <  1±er  both  branch  points  vlU  be  found  on  the  real  axis 
and  vhen  p^  >  1  ^er  they  will  be  located  on  the  Imaginary  axes  In  the  coo^ 
plex  |3  •plane.  Finally,  vhen  (1-S‘)  <  p  <  (Ho*)  one  branch  point  will  be 
located  on  the  real  axis,  idiereas,  the  other  on  the  Imaginazy  axis.  Figure 
9^1  depicts  one  of  the  above  discussed  situations  together  vdth  the  branch 
cuts  Im  =0  and  Im  { 0^  {  so  ihoiia  dotted  and  denoted  by  i*;  and 
respectively. 

With  the  branch  cuts  defined  ue  can  draw  the  four^sheeted  Rleoann  aiu>o 
face  as  shown  In  Fig.  9.2.  We  observe  that  the  cut  joins  the  Sheets  I  and 
li  and  the  out  fx  joins  the  Sheets  II  and  IV. 

9.2  FORMUUTION  OF  THE  CONTRIBimOIiS  10  THE  FIELD  INTEQRALS 

In  the  previous  section  ve  have  deteznlned  the  singularities  of  the  field 
integrals  and,  in  particular,  ve  found  four  branch  points.  In  this  section  we 
shall  determine  the  various  contributions  to  the  field  Integrals  as  affected 
by  these  singularities. 

All  of  the  field  integrals  are  of  the  fom 

W  a  J  F((3  (9.9) 

f', 

and  ve  vish  to  evaluate  these  integrals  approximately  when  the  distance  f  is 
large*  The  saddle  point  of  the  integrands  occurs  vhen  the  derivative  of  the 
exponant  vaniahea,  i.e.,  vhen  /3  ~  9.  By  the  argument  identloal  to  that  in 
Section  4.2,  the  original  path  of  integration  P,  can  be  defonied  to  the  path  of 
the  steepest  dssoents  P  the  integration  can  be  perfonsd  elong  that  path. 
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Nov  ve  fooua  our  attpontlon  on  the  foUovlng.  As  vs  remarked  earlier^  the 
integrands  under  consideration  are  four^valued  since  they  contain  the  radicals 
of  cr;  and  •  The  original  path  of  Integration  passes  over  Sheet  1  of  the 
four>aheeted  Riemann  surface  and  can  be  deformed  into  the  path  of  the  steep¬ 
est  descents  r  only  ^en  at  least  the  beginning  and  the  end  of  it  lie  on  this 
sheet,  in  the  case  when  the  angle  Q  does  not  exceed  one  of  the  angles  of  the 
total  internal  reflection  d,,  or  e.a.  then  in  the  complex  -plane  we  have  a 
picture  as  shown  in  Fig.  9*1.  Here  the  path  of  the  steepest  descents  crosses 
one  or  both  branch  cuts  twice  euid  thus  ends  on  Sheet  1.  Transition  from  path 
n  to  the  path  r  is  then  accomplished  without  any  oomplicatioas.  Ibccept  for 
the  part  show  dotted  in  Fig.  9.1,  the  path  P  will  lie  on  the  Sheet  1.  in 
this  case  the  only  contribution  to  the  field  integrals  %dll  be  from  the 
saddle  point. 

The  situation  is  quite  different  when  the  angle  4)  exceeds  the  angle  of 
the  total  internal  reflection  of  one  or  both  waves.  Suppose  thet 
Og,  <  6  <  952  •  P  will  cross  the  branch  out  f  only  once  and 

the  branch  cut  twice.  Thus,  the  path  P  will  cross  from  Sheet  i  onto  Sheet 
il  and  then  IV  and  then  back  to  li  Wich  is  inadmissible.  When  9  >  9bx>  O5, 
the  path  P  will  cross  each  cut  y,  and  only  once.  Thus*  the  path  P  will 
cross  fFom  Sheet  i  onto  Sheet  ii  and  then  onto  Sheet  iV  Wlch  again  is  inad-> 
sissible.  Ve  can,  however*  construct  a  more  opmplioated  path  of  integration* 
sttpidsMnting  the  path  of  the  steepest  descents  P  by  a  contour  encompassing 
each  cut  in  such  a  way  that  the  beginning  and  the  end  of  the  more  complicated 
path  win  lie  on  Sheet  i.  This  situation  is  similar  to  that  in  Section  4.2 
only  that  here  two  brwcb  outs  we  involved. 

To  this  end  ve  propose  the  contour  shown  in  Fig.  9.3>  For  9  >  0  this  new 

path  of  integration  runs  along  the  path  of  the  steepest  descents*  P  *  then 
•loiig  the  borders  of  the  cuts*  intersecting  each  branch  out  once  thus  entering 
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ShMt  IV*  It  then  oontlnuee  along  the  path  P  and  Intersecting  each  branch  out 
once  nore  it  oonea  back  on  Sheet  1  ufaere  it  oontinuee  along  the  path  P  •  tw 
0<  0  the  situation  is  analogous. 

As  a  consequence  the  complete  expression  for  asgr  of  the  field  integrals 
vLll  eozisist  of  three  parts  and  can  be  written 


(9.10) 


idiere  V^vill  denote  the  oontrlbutlon  from  the  saddle  point  and  W^^^Wid 
Idll  denote  the  contributions  from  the  borders  of  the  cut  assooiated  with 


teandi  point  B,  and  B2  respectively. 

9i^r  oontribotion  th*  ue  re¬ 

marked  earlieTf  each  of  the  field  Integrals  under  consideration  can  be  uritten 


in  the  form 


U  S 


(9.11) 


The  saddle  point  of  the  integraxid  occurs  vhere  the  derivative  of  the  exponent 
vaniahesy  i.e.,  at  ~  6.  The  leading  term  of  the  asymptotic  expansion  for 

(9.11)  mill  then  be  given  by 


F  (®) 


(9.12) 


9.2  b  Formulation  of  the  eontributione  ftnai  the  branch  cats— First  con¬ 
sider  the  contribution  to  the  integral  (9*11)  due  to  the  integration  along  the 
borders  of  the  cut  corresponding  to  the  branch  point  B,  .  Vie  write 

Vl<^>  ■  [fO)  (9.13) 


and  integrate  formallv  on  both  sides  of  the  cut  getting 


(9.14) 


gj. 

**♦  a,^ 

4*0 


(^)  *  F-’  ((»)  -  F?  O) 


(9*15) 


F^^  ( |3  )  denotes  the  vslae  of  the  function  F  ( |3 )  on  the  left  side  of  the 
cut  and  Fi'^  (^ )  is  the  value  of  the  seine  function  on  the  xlj^t  side. 

Next  us  defoxn  the  path  of  Integration  to  the  path  of  the  steepest  descents 
through  the  point  .  The  procedure  is  Identical  to  that  in  Sectlcn  4.2  b 


and  ve  write 

10 


(9.16) 


/a-ee.^%‘ 


E{]r  analOg7  we  can  write  similar  expressions  for  the  oontrihution  fToa  the 
borders  of  the  cut  )r»'*  associated  vith  the  branch^point  .  Ve  have 

09 

v<^)  s  ie  jx  Rf*>(i(i!)e-««.f8wv(8- <»*«)xyadx  (9.17) 

p 

• 

Nov  we  find  the  opntrlbutions  fruo  the  branch  cuts  idien  the  angle  0  <  0. 
Integrating  along  the  borders  of  the  cut  asspelated  with  the  braDcbi.point  B*  > 
ve  obtain 
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*4^ 

B.* 


•  ^ce 

B,- 


(9.18) 


i&ere  Is  givsn  fay  (9^1 5)*  Rolloulng  the  procedure  of  Section  4«2»  we 

deform  the  path  of  integration  to  a  path  of  steepest  descents  through  the 
point  B,~  and  obtain  as  a  result 


«• 

w*?’)  s  ie Jx  R^A)®****? 


(9.19) 


By  analogy  the  contribution  from  the  branch  cut  associated  Vd.th  the  branch^ 
point  Can  be  wltten  at  once,  vie  have 

y^cec)j-  ci«,pco*o«i-  Bmi)  f X  e*w.p4CK(isi‘  >cj|^  x  .  (9.20) 


Vie  can  combine  the  results  of  (9*16)  and  (9.19)  ns  veil  as  (9.17)  axtd  (9*20) 


and  urite 

a  i«ii«.feo*fiei-e„) 

eno 


f 


xR<tn<5)€“ 


s  ie  ce*0Bi^ei!i) 


eno 


O 


I X  e  **‘»f  (w-«m)  xJi  dx 

O 


(9.21a) 


(9.21b) 


9.3  EirMUATION  OF  THE  FIELO  INIBQRALS  AT  THE  SAPULBi>P01NT 

In  this  section  ve  shall  apply  the  result  of  the  formulation  (9.12)  to 
find  the  saddle-point  eontribation  to  the  field  integrals.  Because  of  the 
oompleacity  of  the  results  ve  shall  only  find  the  leading  term  representing  the 
xmdiation  field. 
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9^3^a  -  ihe  field  compouent  Ei*i  ^iODiicatioa  of  (9*12)  to  the  integral 
expression  for  E^*^in  (9>1a)  gl\res 


(g)ooee5^ 
/Jir  /kTip 


(9.22) 


vbere 


B  ^g<  fe(9)»  4oosQ]  *Pet(0)  $0609]  e*'***^^^^^ 


(9.23a) 


>^*(e)®P6,(»)  [ff.(«)+ooa9]  [ei(®)4<oose]  -Psi(9)  [ei(®)+oos9].[or(«)+i‘oos»] 


(9.23b) 


*‘eoi  (®)  ®  2  («)  +  ✓&*  (0)  +  ^7^sin‘® 


(9.23e) 


^^(er’Sio*9)  *  Pi..,.(®) 


(9.23d) 


R(®)  S  ^  ’»(o..»)ain® 


(9.23e) 


9.3  b  The  field  oomponent  ^Application  of  (9.12)  to  the  integral 
eQqairesaion  for  E^^  in  (9.1b)  giree 


i*^  ^  OJ  I 


P_  (®)aiaOoosO 


lAere 


f  m  ^^(^)'*'0Pb9]  t^Mvq-.CeJ  ,  [e;(®).»ooe9]  e**«.*^g»Cg) 

>^,(9) 


(9.24) 


(9.25) 


•Bd  q:(9)»  ^(9)f  and  ^(9)  vmps  defined  prefionnly. 
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9»3  0  Tha  field  conponent  — Ae  can  be  noted  from  the  integral 

representation  of  In  equation  (9*1  e),  the  integrand  in  aero  at  the  saddle 
point  /S  B  6.  tbuSf 

«  0  (9.26) 


9.4  EVALUATION  OF  THE  FIELD  INTEGRALS  ALONG  THE  BRANCH  CUTS 


In  Section  94  2  ue  have  foxnulated  in  general  terns  the  contributldn  to 


the  field  integrals  due  to  the  Integration  along  the  borders  of  the  branch 


outs.  In  this  section  ve  shall  apply  this  formulation  to  find  the  branch^cut 
oontributlons  to  the  field  Integrals.  These  branch.>aat  contributions  are  the 
aouroes  of  the  lateral  wares  ae  ve  remarked  earlier. 

9.4a  Ihe  field  conponents  EvV  and  sSo*^  —.We  apply  (9.21a)  to  the 
integral  representation  of  £  in  equation  (9.ta)  and  noting  that  the  lnt». 
grand  is  an  even  function  of  /B  ,  we  obtain 


g(e.) 


2 j 


(9.27) 


idiere  by  (9.16)  and  (9.1a) 


‘  \  (9.28) 

eos|8Xe< ■*^«9Sd)^Cei30  J 


(9.29a) 
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Nov  ve  imist  axpresB  In  teroa  of  z.  To  aooonpllBih  this  ve  set 
(i  *  Qj,  erexsrwhere  except  in  o;  %Moih  becomes  eezo  Mhen  /a  *  e,, ,  Ve 
evaluate  o;  then  to  the  next  hl^er  appraadaatlon*  lb  this  end  ve  note 

sln»/3  Sin*  fig,  +  ix*  Sin  fi,,  cos  fig,  (9.30) 

and  we  find  0:  to  be 


o;  a  -  e^^^A 


Forthetnorei  w  note 


^  ain  2  fia, 
C  ^  e  ^  ^ 


a  -  ^  7 

Pg(«;)  a  -  <"1^ 

(«i)  ®  7  (^♦7). 


(9.31) 


(9.32) 


SabsUtating  the  abore  results  into  (9.22)  and  carrying  out  the  necessary 
algebraic  operations,  ve  find  the  leading  term 


vbere 


a  2  C 


(O 

yo 


(9.33) 


.0) 

'yo 


(9.34a) 


^Cl<oKb„*  Pb,  +  0-C)Ca+7)€*^*»«*»''^  >  •coafig,  ] 


and 


Si  ■(^■^7+^)' V^*^rr+77?T+(e  +  7+ ^).y/f::^?’5^  (9.34b) 

<  v^iHaty)  [(U6+7)  f(a+7+t)+(«+7+<^)  ✓HatY)  j. (9.34o) 


Perfox^ng  the  integraUon  ehere  ve  use  the  refiilts  of  (4.31b)  ve  obtaiB, 
finally. 


iC»0  -  .M  e-*:^ 


(9.39) 
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Mow  VO  oToluato  the  mbo  Intognl  along  the  borders  of  the  oeoond  eat.  Ue 
write 

.^■.1  .  " 


yo 


2T 


((»)•- 


kfS4;oei-®»0»‘>iic<x  (9*36) 


iMiere 

a 

*yo 


PP-  _P#rg.)feqk  ^  FWf  egjrer.  ^cos#? 

([cr,^eoc(S)(-< 


-«» -t-Scoft  OJ  ■*■  cos(3)^»i  •♦■^’cosia)  Pgf  fli) 

(JO;  +  cdS|l)^<ra^  <’co5|3i^f^^  «;)  “  C*5 •»■  coifiX^*  9cos|S )R| C^}  f 


(9.37) 


Where 

(9.38a) 

Sin  s  v^€-  ^  .  (9.38b) 

Vs  express  in  tents  of  x  we  set  (S  »  eveiTWbere  exeept  in 
lAloh  bsoomes  sero  when  (3 «  •  Ve  evaluate  then  to  the  next  hif^er 

approxlBatlon,  getting 

^  a  •  \/  -■  X  (9.39) 

fnrthennref  we  note 

($+€«.y) 

fg{^)  ^  ^  (9.40) 

?»(«&)  ^  ^7. 

Sttbstitatiag  the  aboxe  resalts  into  (9*37)  ae  obta;ia  for  the  leading  term 

R^i  ®  ZG^l^  e-^^  X  (9.41) 


idisrs 


136 


(9.42) 


' [4 k.kb.j  -0*1  +  (I -?X«- 7 -CO* ] 


■ad 

»gj»(<+€-^  )  V^  (^♦<-7)  ♦  (  '^*+■-7)  ^l-(«  -  7 )  (9.43a) 

^82*  ?/1-(a-7)  [(Hfi-y)  ^^(t  +  «-7)+(( +*-7)  /i-(e-7)  j  (9.43b) 


Performing  the  Integration  vhere  ue  use  the  reealte  of  (4.31b)  %»  obtain, 
finally, 


•  (9%)  a  *  f.(X)  *<>><« f go* rtei*-o,»jg^t^ 

'•  l/51f  '•  [k.,.lll(l8l  jj't  ■ 


(9.44) 


9.4  b  the  field  ooapwaeatfl  and  BiS*^..*We  mpfij  (9.21a)  to  the 
integral  representation  of  E*^  In  eqnaUon  (9.1b) 

*eo  ■  "^tZV cic.f.cos(t»K e»^ e" *»»f *«•’>» ee.) »»?i  4xj(9.45) 
vhere 

<  6  ♦  €o*  A>r®x^  ?eosp>^  gp*(*)6 

„  Cen»co«6)e^*»*ke;  -  4-ce»fl3gvMteit  ) 

Cw;  ’f  goef5)('^-t’<co»jO*%C’n>»Ca|  ’►cpa^)(’or*<'oiaii5)%C'ii.)  1  (9»4(6) 


a^tlttttlng  the  acpreeaione  for  ^  ,  Pe(e;)  and  ?«(«&)  ftoi  (9.»)  ve  obtain 
for  the  leading  term 

-  *2cy,  ooe(  |QI»  9»,  )o«^V* 

^  7/«+7 


IE 


(9.47) 
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«h«re 


rC')  -  >/r[stKu2e,.l^t  COSest^Z-Kf  ?•»#•»  »7) 

c,.,  ZbVy^-^JV 


(9.40) 


and  tdiara  a,,  la  givm  Iqr  (9.34b}  and  b^,  bj  (9.34o).  Parfoning  tba  iatagra* 
tlon  idiara  ve  use  the  raault  of  (4.31b)  va  obtain 


e-^  i  cosdSi-  )  /g.Qi 


fo  parfpm  the  differentiation  we  use  (0«35)  and  note 


*a.  “■  *  1 

for  the  leading  ten.  Thus*  finaUgr 


(9.50) 


«<a.>  ^ 

••«eIo  ' 


I  ^<0 


£  h*^  Sivx  O*!-  ©ei)] 


(9.51) 


lov  wa  anlaata  tba  ssm  IntagnX  along  tha  bordara  of  the  aaoond  bnuub  out. 


Va  wtta 


9(ht  ^ 

2Tk,a 


abara 


ti^V±£)s  aoa6epsA>M)^  _ 

•oV  P'  ^(e^+ccsaJ&^k-*-  ?eo»il)»VC«) 

*  (<^~«c«a>X«k'^^aoig/fbCa(>^^^  X 


(9.53) 
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tebstitutijag  for  o;  ,  P,  (a;),  and  P((er,),  tha  aaqpraaalon  found  In  (9*40)f  va 
obtain  for  tha  landing  tan 


«(t>  ^  .  2  Gia  ooa(  Iftl  -  86i.)a“^^  - 

••o  ■  ♦  ;,rr-  -m,  * 


(9.M) 


aC»)^  Vr  [a^2i9»  J^4oveBt»y^Ktf-»*  (9.55) 

•[tk.kbgi-ati  I “ ^ )C«- 7)ce»s Oji C 7>  ] 

and  ag2  ia  glvan  bgr  (9«43a)  and  Iqr  (9*i3b).  Parfondng  tba  intagntion 

ataara  aa  uaa  tha  raaulta  of  (4«31b)  aa  obtain 


^Ccu**. 


eos^iai  -  ) 

[kaf  abiaObi*daa)j^  J 


(9.56) 


fa  aany  out  tba  Indioatad  diffaraatlatlon  aa  nota 

for  tba  laaditig  tana.  Ifaaay  aa  obtain  finaUgr 


-ca.)  _  _  -fcoM..X  PiV  0Pi(l«l  -  9„) 

*•0^*  "’  M*  -  - 


(9.57) 


(9.58) 


9.4  0  Iha  fiaid  oonponenta  and  -.-Tb  find  and  aa 
ooold  prpoaad  aa  bafora.  Ihia  la*  boaarar>  not  naeaaaaxy  if  aa  notioa  tha 
adAUjcrity  of  tha  Intagranda  oorraapooding  to  Be^and  B^a  in  aqoaticma  (9.1b) 
and  (9*1o).  Ibua,  B^p^  oan  ba  obtainad  frma  B^p^by  putting  ±  ain(l9l  «  P»i ) 
in  pliioa  of  ooa(  181  ^891)  in  tha  nnarator.  Sinilarlj  B^^  oan  ba  obtainad 
frail  Igr  pitting  x  aln(  I8i  •  8u)  in  plnoa  of  ooa(  191  «•  in  tha 


ttUMrator.  Ihusy  we  oan  write  4*— 


jceo  s  »  .  «ln(|g|  e-^^ 

^  >/5T€  [k«f  ain  ( Iftl  -  «•, )] 

and 

jrt«>  W^ICeV  alnClftl  ■  <>»>) 

^  V^«  Lk,f  aia  (l®l  * 


9.Z,  d  Pbaraical  Interpretation  of  the  lateral  watee—Bie  oontrlbution  of 
the  branob-out  Integration  found  abore  haa  a  well  defined  twinning  In  the 
Idqraloe  of  wave  propagation  In  layered  media  (5»p*270).  It  la  the  nathaaati* 
eel  equivalent  of  the  i^aloal  phenoaenon  of  lateral  wavea  whoae  hlatorloal 
aooount  we  have  given  In  Section  1.2#  Here  we  aball  attenpt  to  eocplaln 
piqraloallgr  the  Mohaalan  of  the  lateral  wave  foiaatlon  on  naipMitopiLaaae 
laterfaoea. 

Vd  thla  end  we  oonalder  the  geonetiy  of  the  Plgt  9«4«  We  note 


X  «  d,  +  g,  Bin  9^,  (9*6la) 

a«g,  ooaOe,  (9*6lb) 

an  a  oonaequwiee  of  whioh  we  obtain 

eoa  (9  *  9*,)  «  +  7  d,  ko  g.  (9*62a) 

^  ain  (9  -  9b,  )  «  /  1  -  U +7)-  d,  (9.62b) 


Oalng  the  above  reaulta  we  oan  reoaat  the  fltU  eonppnenta  due  to  the  bran^ 


* 


4 


flgu^  9«4  Conowming  tbe  of  «  liRtopal  Itevo 


U1 


In  the  above  ve  zvoognlie  that  /e  +y  and  ko  -y  are  the  propaga^ 
tlon  faetora  of  plane  ordinary  and  extraordinary  raya  reapeetively  In  a  laag- 
netoplaama  vhen  the  direction  of  propagation  oolnoldea  with  the  direction  of 
the  ateady  loagnetio  field.  Thus,  the  total  phase  of  each  one  of  the  lateral 
wavea  oonalats  of  two  parts.  The  part  for  exanple,  k  ^  »  tells  us 

that  the  wave  travels  along  the  interface  In  the  nagnetoplaaaia  through  a 
distance  d,  with  the  phase  velocity  of  Cy/ s/e  +  »ji  .  At  point  A^  Fig.  9«4>  the 
wave  changes  Its  direction  and  travels  along  the  path  A  -  B  With  the  phase 
velooitiir  c  $  Which  is  the  phase  velocity  of  free  space,  nue  other  lateral 
wave  associated  with  the  second  branch  eut»  B2  >  follows  a  similar  pattern. 

The  complete  picture  of  boUi  wives  is  showi  in  Fig.  9*5  where  we  assumed 
for  convenience  that  the  source  is  located  at  the  interface.  The  ray  OB  is 
the  direct  ray  from  the  source  representing  the  radiation  field  as  found  fron 
the  saddle-point  integration^  The  lateral  wave  represented  by  the  ray  OA 
pcopagates  along  the  Interface  with  the  phase  velocity  of  the  ordinary  wave 
In  the  aagnetoplaaaa.  Binoe  this  velocity  is  different  than  the  one  Allowed” 
by  the  free  space,  there  is  a  disturbance  created  along  the  interface  in  the 
air  as  a  result  of  which  a  new  wave  is  produced  represented  by  the  ray  AB. 
Sinoe  the  spatial  period  of  this  disturbance  along  the  interface  is  equal  to 
.  the  wavsi-length  of  the  ordinary  wave  in  the  magnetoplaasaf  the  wave* 
length  in  the  air  can  "fit”  this  period  only  if  the  angle  9^,  between  the  nor* 
■si  to  the  boasdai7  and  the  direction  of  propagation  is  such  that  the  wave* 
length  in  the  air  ^^^ain  •  This  is  Just  the  direction  of  the  let* 

enl  wave.  The  seas  argument  holds  for  the  second  lateral  wave  represented 
by  the  ray  OA,  The  resultwt  field  at  a  point  B  is  a  saperpositiott  of  three 
wtves  as  ve  ooneluded  ewlier  in  equation  (9,10}, 

Us  diiiyns  of  fig*  9*5  dtopiots  the  sitaitiott  whan  i  ±er ,  i.e*»  whan 
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(9^64) 


If  the  ebotre  Inequality  la  not  eatlefied,  one  or  both  of  the  lateral  waTee 
uLll  dlaappear.  Physloally  this  ocoure  uhen  either  0,, ,  or  or  both  are 
outaide  the  range  of  interest  0  to  ^2« 


9*9  1^  I^WEB  FLOW 


In  Uxe  preoedlng  seetion  ve  found  the  eleotrlo  field  oomponents  in  the 
air*  the  fore  of  these  oonponenta  tell  a  great  deal  about  the  physios  of  the 
pdrohlea*  Howerer ^  still  a  better  insist  into  the  pvobleei  oan  be  obtained  by 
eonsidering  the  direction  and  aa^tude  of  the  energy  serried  by  the  vavo  re^ 
presented  by  these  field  ooaponents. 

Mow  the  tliae  averaged  energy  density  oarried  by  an  eleetroeagnetio  isnre 
is  given  by  Poynting  vector 


T  (9.65) 

thieh  oan  be  witten  in  coapporot  fore 

^  Re  (I,  -  Sy  H»  )  (9.66a) 

S,  =  Be  (By  H»  •  H*  )  (9.66b) 

Sy  ^  ^  Re  (B^  ^  I,  )  (9.66o) 


Thus*  before  ve  deteniLne  the  ooeponents  of  the  Foynting  veetor  «e  first  east 
find  oonponimts  of  the  aagnetie  field  veetor. 


9.5  a  The 


obtaiiiahle  fras 


field  vsetor  in  the  air  is 


(9.67) 
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In  nliat  foUovi  \m  nball  onlgr  b*  oonMdRiad  with  thn  radiation  fiald  ainoa  the 
lateral  field  Is  of  seoond  order.  Unie^  In  the  radiation  field  0  and 
aoraoter  to  the  first  orisr 


iko 

0 


(9.68) 


litiah  anablea  us  to  find  the  ooi^onents  of  the  lagnetio  field  teotor  as  foyeast 


(9.69a) 

2o 

Hy^«  1^  (9.69h) 

3o 

there  is  the  free.'Spaae  li^edanoe. 

9.5  b  The  cottponenti  orthe  ^  results  of  (9*69) 

together  idth  (9.66)  we  obtain  only  one  eoaponent  of  S 

V  (  i*yol*  +  )  (9.70) 

2  Zo 

Hov  using  the  results  of  (9.Z2)  and  (9*24)  we  obtain  for  the  first  order  temt 
8  "  ^  ♦  sin’olFiol*  )  (9.71) 

there  fy^  is  glren  by  (9.23)  and  tgr  (9*25) .  We  note  that  the  esaantlal 
difference  between  the  abors  result  and  the  result  obtalnod  Croai  a  oprresjf»nd« 
lag  Isotropio  ease  is  the  firesttioe  of  the  seoond  tens. 


power  pattern  represented  by  equation  (9.71) 
was  OQBpjited  and  plotted  In  Figures  9*7  and  9*8  for  le  eleotrlo  ourrent  line 
souroe  situated  at  the  lower  edge  of  the  ionosphere.  For  the  purpose  of  thle 
■»— pT*  the  lower  edge  of  the  ionosphere  was  assuned  to  be  sharply  boundedi 
hating  an  eleotron  oonoentratloa  H  ■  790  ^otrons  per  eUhie  oentlneter  wid  the 
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eiufth’s  fflagnetle  field  ,4  gausa.  For  the  sake  of  eoaparlaon  the  power 

patterns  oorreapondlng  to  R  =  0  we  also  Shovn  in  the  Same  graphs. 

the  examination  of  the  radiation  patterns  In  Figures  9.7  and  9.3  reveals 
that  the  treninisslon  of  energr  from  the  source  through  the  plaana  and  Into 
the  free  space  is  generally  enhanced  by  the  presence  of  the  steady  magnetic 
field  especially  at  wave  frequencies  below  the  plasma  frequency.  At  frequen* 
des  weH  above  the  plasma  frequenqjr,  the  action  of  the  steady  magnetic  field 
becomes  Inslcpilfleant  idilch  is  well  understood.  Another  interesting  phenonH> 
enon  is  the  appearance  of  a  peak  in  the  radiation  patteins  at  angles  ooxres^ 
ponding  to  the  oritloal  angles  (see  Fig.  9.6).  As  we  showed  In  equation  (9.7) 
these  oritloal  angles  O^i  and  related  to  the  Indices  of  refraction  of 

the  ordinary  and  the  extraordinary  plane  waves  In  magnetoplasma  vtien  the  dlr^ 
eotlon  of  propagation  Is  along  the  Steady  magnetic  field.  Since  the  critical 
angles  are  generally  different  (except  at  frequencies  high  enough)  for  a 
plasma  with  and  without  the  steady  magnetic  field  (see  Fig.  9*6)  the  oorres^ 
ponding  patterns  peak  at  different  points. 

Finally,  we  may  mention  that  at  certain  frequencies  two  peaks  in  the 
radiation  pattern  can  be  expected.  In  our  example  this  would  occur,  for  In- 
Stanoe,  at  a  frequent'  of  2  z  10^  ^cles  per  second  >diere  both  oritloal  angles 
0^,  and  9b2  (Fig.  9.6).  We  were  not  able  to  show  this  phenomenon  due 

to  the  unfortunate  ohoice  of  wave  freqaenelee. 

9.6  IHS  VAUItllZ  OF  THE  HIGH  FREQUBICI  APPRQnKATlOH 

In  solving  the  dipole  probleme  we  found  it  neoessuy  to  use  the  high  fTe^ 
qunhcy  appitnCUiatlon  before  carrying  out  the  integration.  Since  no  rigorous 
solution  to  the  dipole  problems  is  available  it  was  somewhat  difficult  to 
•ppraise  the  vall^ty  of  such  an  apprsadmation  aftw  tha  integraUpn  was 


7  «•  Power  Pattenx  in  Air  of  an  ©ectric  Current  Une  Source  in  Magneto^ 
Ha^Noraal  to  the  Line  Source,  N  =  750  eleotrona  per  cubic  eentteeteri 
-^ters,  0,  .4  gauae,  f  a  10*,  10^  cyalaa  per  aeoobd 
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MXTiadl  out* 

la  the  prasont  ca8«)  ho%feTery  it  vaa  not  nooossazy  to  um  angr  appiroxina^ 
tiona  to  arrive  at  the  final  aixaver.  Tfauay  our  preaent  aolutloa  ia  rigorous 
and  it  can  be  effeotiteljr  used  to  test  the  validity  of  the  hig}i  fretiuen^ 
approsdnation. 

to  thia  end  ue  shall  perforn  the  following  testa*  Firsti  we  apply  the 
hi^  fMqoenoy  approxination  before  the  integration  was  earrled  out  as  was 
done  in  the  dipole  problens*  Seoondf  we  apply  the  hi^  frequency  to  the 
rigorous  solution  found  in  this  chapter*  Finally*  eonparing  the  results  of 
the  first  netfaod  with  that  of  the  second*  we  estahliah  the  talidlty  of  this 
approxination* 


9.6  a  IntrodueUon  of  the  aparoxi«i«t^“P  bigh 

frequanoy  approadmation  of  equation  (3*7)  in  the  notation  of  the  present  pro^ 
VLan  eaounts  to  setting 

2  /n*-  8in*/5 

Introducing  this  into  the  expression  for  the  denominator  in  equaticm  (9*^) 
and  neglecting  the  second  order  terms  in  k  we  obtain 

;2Ka?8in|^  (  ♦  cos^  )(  /n*-oi?jS  ♦  n*'oos/a  )  (9*73) 


Furtheznorei  it  can  be  shown  that  the  integrands  in  (9«1a}»  (9«1b)f  and  (9*1c) 
rsdttoe  to  the  foUowingt 


'yo 


% 


(9.74a) 
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’•e 


K6*Jm-«i.I  a  f sLr^6_*^0^ )le *•  *^* ** ^^**'j^**y 
4Tri<o  J  ♦‘CdspX*^**st«.*(S  +  yx.*cos(S  ) 


.g£W,fco»<j»-«; 


(9*74b) 


*fo  = 


TKfaj«i*aI^  ^  I  ^Qsftat 

4Tfk.  ^JTSSsT 

r 


t  ees(9yv^»x'-stKv‘/r  ♦  tl*co»/a  ) 
g<:lc^feos^(4-#>e<^  . 


(9^746) 


the  saddX^polnt  lategration  ei&  be  earrled  out  at  onee.  tlaing  (9*12)  wa 
obtain 


S  c  « 


•yo 


cot4sl 

v/sr  ' 


✓a*-  8ix^9  ooa  0 


(9.75a) 


kW 

•oo 


fcKcOiJuI  /scrtZeCi-  IMUS 

4  v5n  4  T^^s5?3^^%^S3cSS9  +  c6s  0  )(V^5?s7r??^*co«ai 


S^V  *  0  .  (9.75o) 

Va  note  that  the  Integrande  In  (9.74a),  (9.74b),  and  (9.74o)  now  oontaln 
a  branch  point  at  ain  Qp  b  q.  to  pei^osai  the  branah-oat  Integration  «a  naa 
the  fomalation  of  Seetion  9.2b  ag^  obtain 


g<e^  ^  i  y/n  (1^ik,h l/l^)  a  ^**.f  go«C»ei» e,) 

^  ^  ^  9,1  ]  H 


(9*76a) 


T  Kujt*.!  \/^(1*ihoh  )  a^  qoa(  191  *  9a  )t*-**.f^^^^*^*^**) 

2  \/io?ldSlimrT 


(9.7i5b) 
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«<•)  _**»<**J<**I  /n  (1-ikoh  >/ 1-n*  )e‘’^  8in(  191  -  ^  ^ 

■  ^o  "  .  ■  '  ■  ‘  •\y»ToC) 

2/2f(1-o^)^^  lint  101  -0») 


9.6  b  iQttDdactlon  of  the  apn^^nu^ti  after  laid^Uoa^IntroducUon 
of  tbo  hl^  ftroquoncy  approxlBatlon  In  the  foxn 


/n^  -aln^O  ± 


Ka  alj  0 _ 

2  /n*^-r^*« 


(9^77) 


into  the  field  expreselone  of  (9.22)  and  (9.24)  gliree  jrosulte  that  Identi- 
eal  with  those  of  (9.75a).  (9.75b),  and  (9.75o). 

The  appUeatlon  of  the  sane  approodjoation  to  the  components  of  the  later* 
al  field  In  (9.35) i  (9.U)t  (9.51),  (9.56),  (9.59)  and  (9.60)  affects  the 
•HpUtudes  of  these  components  only  bdt  does  not  alter  the  order.  Thus,  all 
of  these  oomponents  remain  to  be  of  the  second  order. 


ia^.g— V^gXop  9f  YftliditT  of  the^hlgh  freonenpiY-  the 

preceding  seotlons  ve  applied  the  hi^  frequency  approximation  before  and 
after  the  Integration  and  obtained  the  corresponding  results.  We  found  that 
In  the  oase  of  the  radiation  field,  vhloh  Is  given  by  the  saddle^polnt  contri¬ 
bution,  the  results  did  not  differ.  Hovever,  In  the  case  of  the  lateral  waves 
glvea  by  the  branch^cut  Integration,  ve  found  that  the  results  differed  sub¬ 
stantially. 

The  reason  for  this  apparent  dlsorepency  Is  as  foUovs.  The  high  fte- 
quenoy  approadbaatlon  Introduced  prior  to  the  Integration  was  of  the  form 

/n*- 8ih^(5  t.  — •  (9.7$) 

2  8in  «|| 

Obvioualy  this  apin’Qxlmation  Is  not  valid  when  the  denoadnator  of  the  second 
tita  vanishes,  l.e.,  ^en  sin  /t-  di  n.  gov  In  the  saddle-pplnt  method  of 


integration^  the  process  of  integration  laerely  changes  |%  to  the  polar  angle  0 
leaving  everything  else  unchanged.  Thus,  the  result  is  valid  everyvdiere  eac> 
cept  vSiere  sin  Q  ti  and  it  Is  independent  of  t^ether  we  apply  this  approxi«> 
^tion  before  or  after  the  integration. 

In  the  Case  of  branch-^cut  integration,  however,  the  situation  is  quite 
different.  Whereas,  in  the  rigorous  form  of  the  integrands  contain  two  branch 
points  at  07  s  0  and  0,  the  approximation  (9*73)  in  effect  coalesces 
these  two  branch  points  into  one,  the  branch  point  of  plasma  in  absence  of  the 
steady  magnetic  field.  The  branch-cut  inte^ation  in  this  ease  is  then  not 
valid  and  its  result  must  be  disregarded. 

9*7  GLOSDBE 


In  this  chapter  we  found  the  field  components  in  the  air-region  due  to  an 
electric  current  line  source  in  a  magnetoplaaima  when  the  steady  magnetic  field 
is  normal  to  the  line  source.  In  particular,  we  found  that  the  field  in  the 
air-region  consists  Of  two  main  components,  the  radiation  field  and  the  later- 
al  field. 

The  radiation  field  contains,  in  addition  to  the  oomponents  present  in  the 
corresponding  isotropic  case,  additional  transverse  oomponents  of  the  electric 
and  magnetic  fields  vdlch  render  the  resulting  fields  to  be  ellptleally  polaz^ 
iaed  everywhere  except  right  above  the  souroe. 

The  lateral  field  is  of  second  order  as  in  the  ooneiponding  iaotroplo 
case,  however,  unlike  in  that  oase  it  may  ooneist  of  two  waves  rather  than 
one.  Ibis  is  due  to  the  fact  that  the  magDetoplasiia  pofseesea  two  distiBOt 
indices  of  refraotlon,  eorrespondlng  to  the  ordinary  and  eictrsordinary  waves, 
eaOh  one  of  these  waves  possessing  its  own  qritlssl  ssgls  wHsii  is  tbs  sngls 
of  the  total  internal  refleotion. 

Finally,  by  introduoing  the  hi(^  fraqusncy  approxiaatlom  la  ths  rigorous 
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ex{iresslons  before  and  after  the  integration  and  coaparing  the  oorresponding 
resulta,  we  were  able  to  detexvlBe  the  region  of  validitjr  of  this  appitadjna- 
tion.  it  appears  that  the  final  expression  for  the  radiation  field  is  not  de¬ 
pendent  on  idiether  this  approximation  is  introdaoed  before  or  after  the  inte¬ 
gration.  Tbusj  the  proeess  of  its  Introduction  prior  to  the  Integration  is 
valid  as  far  as  the  saddle-point  method  of  integration  is  eoneemed.  In  the 
case  of  the  branch-out  integration^  however}  it  is  not  so.  Here  the  corre- 
spondlng  results  differ  significantly  and  therefore  the  branch-cut  Integration 
of  the  apprpxliBate  integrand  is  not  valid* 


w 


CHAFES  10 

TIS  StiAOf  MAONSIIC  fllLO  PARAULSb  TO  THE  UHl  SOUSOB 

Ill  this  chapter  va  shall  be  eoneerned  id.th  finding  the  electrLo  and 
natie  fields  due  to  an  electric  ourrent  line  source  located  in  a  nagnetofOiaiBa 
shea  the  steadgr  aag^tic  field  is  in  the  direction  of  the  line  source* 

10.1  FOBKDLATION  OF  TEE  PBOBIJSM 

the  geoBeti7  of  the  probleiii  is  ahovn  in  Fig.  10.1.  the  horieontal  plane 
ISO  coincides  vlth  the  interface  betveen  the  anisotropic  hoBogeneous  plaaaa 
and  air.  For  convenience  ve  shall  oall  the  pdaSMa  aediuB  (1 )  and  the  air 
BcdiuB  (0).  As  before*  we  assuae  that  both  aedia  hare  the  saae  aagnetic  in^ 
ductire  eapadtj'  of  free  space*  ^  .  the  steads  aagnetic  field  will*  in  this 
case*  be  oriented  along  the  line  source  idiereas  the  line  spuroe  as  well  as  the 
steadar  aagnetic  field  will  be  parallel  to  the  interface. 

10.1  a  Fm^oinnwfail  eQuationpM»Ihe  solution  of  the  present  boundazy  value 
probleB  entails  solution  to  the  Mezwell'B  equations  in  the  following  fora: 


^  a  I,  e  ico^.B,  (10.1a) 

^  z%  ♦  X^(y)^(s+li)  10.1b) 

^  f  (lO.lo) 

«  0  (10*14) 
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^anu^  we  transform  aquation  (10«5)  as  followat 

\/fir 


(10.7) 


iinea  the  vaalahlne  of  the  Integntad  pai^  la  aaaurad  lAen  the  radiation  ooih> 
ditloii  la  aatlaflad.  the  parameter  a  it  gitfan  bgr 


a  *  /  -  ot^ 

Im  ^  a  }  0  . 


(10.8) 


The  aolutlon  to  the  differential  equation  (10.7)  la  weHaknown  and  we  write  it 


at  onoa 


2/i7r  a 


(10.9) 


Inverting  with  raapeet  to  the  «z,  trantfom  variable  we  obtain  the  detired 


raault 


f d«, 

iir  2  a 


(10.10) 


The  foxB  of  the  conplamentaiy  field  in  the  plaaoa  inaMkUataly  auggeats  Itaalft 


we  write 


4^ 

m  ^  f  A,  a^<««)'-  d«» 

ATT  J 


(10.11) 


and  tiadlarly  for  the  air-region 


i  I  A.  a  da* 


(10.12) 


la  {  a,  }  >  0 


(10.15) 


the  boundary  oonditlona  to  be  aa'^afiad  at  the  interfaoa  f  «  Q  raqulra 
oontinuitj  of  tangantial  ooiqionanti  of  the  alaotrlo  and  aagnatie  fl.ald& 


15S 


thii  marnim  to 


lUoh  dotcoidMi  the  odefflelents  and  A,  as  fbllovst 

2  a  *^  _ 


So  a  a 


A,  B  .g^^ga  ..  .  « 

a  'i‘  s#  a 

BuiSi  integ^  foioa  of  the  elaetElo  field  eoaij^anta  are 


(10.14) 


(10.15) 


4ir 


a  +  a.  J  * 


(10.16) 


*.o  *  -  I  i;^-  -  e‘^-“y-*-*)  d(y^  . 


(10,17) 


10.2  IHS  PIBJ)  00MF0NBI13 


Having  found  the  Integral  representations  of  the  field  ooaponents  in  each 
regioni  vs  shall  now  reduee  then  to  the  fom  sui table  for  nuaerioal  ooaput»- 
tions.  In  this  prooeas  vs  shall  avail  ourselves  of  the  teeluiquea  of  saddle- 
point  and  brandi-oat  integration  developed  in  the  provioos  ohapters. 


10.2  a  ilelda  in  tranafon  the  integral  espreasion  (IQ.  16) 

to  dplindrloal  aoordinatea  in  both  transforn  and  oonfigure  I  n  apaeea  using 


Qixm  k,sln|S 
jBf.,*  sin«,,» 
i  -r  h  s  ooB  0, 

a  «  h  «  fv  oos  9a 


(10.18) 


«ad  obtain 
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■«l 


eoft^  ♦  vy^x~  Surr?l^  J 


(10.19) 


Iho  flLrst  two  taim  la  tao  aboiro  IncogTsl  oaa  be  reoogDlzed  lioDed^tely  ae 
IntegtaiL  rapniee&tatlona  or  ruaetloae.  Bittay  we  can  rewrite  (10.19) 

la  tbe  foaa 

1__^  s  „r^ytl.  I  Hi'  (k,  )  •  Bi’  (k,  f,  )l  -  BiT' 

wti«re 

J^r)  ^  coe-.I  r  eoa  .  .  (10.21) 

2Tr  J  ooa/b  +  -  iln*^ 

n 

The  integrand  in  (10.21)  has  do  poles  in  the  ooiqplsK  ^  ^plane  but  it  has  a 
tarsnOh  point  Where  the  radleal  ( *  sin^^  vanishes,  that  is  la^ere 

sia^  *  -  '  (10.22) 

For  |n|  >  1  this  branoh  point  vill  lie  within  the  region  of  interest  and  in 
defoming  the  path  of  Integration  to  the  pa^  of  the  steepest  descent  through 
the  sfljdlT  pointywe  Bust  add  to  the  saddle^point  contribution  the  oontributlon 
fTcB  integration  along  the  borders  of  this  branch  cut.  Following  the  proce¬ 
dure  of  the  previous  chapters  we  write 

!<»•)  m  u  (0  -  0.  )  (10.23) 

•m  “  ■  xi  *  '  S' 

Where 

0^  ■  arc  sin  («|*)  (10.21) 

and  stand  for  the  saddleepoint  oontribution  and  stand  for  the 
hraneh''oat  contribution  to  the  reflected  field.  Using  (9*12)  we  obt^  for 
thf  iTty"g  teni  of  the  ssfln1e~in^int  contribution 


£frs)  ^  I  ^  _  oog  6  ^  e J-^**'/*  ~^) 

\f^  (oOBfi  +  8li*«)  v/JTK 


(10.25) 


fo  evaLuate  the  braneb^eat  oontslbation  wi  use  the  foxoalatloa  (9^21  a)  and 


vxlte 


(10.26) 


uhere 


R  (%a  )  =  IL  (*/a  )  -  (±^  ) 


(10.27) 


Hoif  Ml  find 


R  (=^/S  )  * 


Goap  s/y^*-  -  ' 

T  [o6a*/a  ••(34v*^  ain*/3  )] 


(10.23) 


Ue  can  put  (6  >  0^  erexyidiere  except  in  the  ndleal  ubioh  vaMshes  at  that 
point,  lie  oompute  it  then  to  the  ziext  hii^ier  order  of  a|n>>wiiBation  and 
obtain 


vhere  of  the  tuo  possible  signs  of  the  radioal  ve  have  ^sen  the  one  that 


guarantees  la  (^*  -  sin^^s  )i'  |  >  0.  Va  find  for  the 

B  )  ~  -■■■— » 
ttW  *  1 


(10.30) 


Sobstituting  the  aboTe  result  into  (10.26)  and  perfoxaing  the  integration 
ahere  ve  use  (4«31b)  ve  obtain  finalij 


^(ra)^  «  oJM.  le*^ ^ 

nta  (IW 


(10.31) 


-J-W  K"! 
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and  note  tbat  the  latenl  field  is  of  seoond  order. 

To  evoiMrlie  the  abot^  reanlta  we  write  the  first  order  electzlo  field 
in  the  pl^ff*** 

2-  v^T.  I  /r  /K  ®  ♦  /kf  -  Sis'll  i.  ^ 

the  Biihetie  field  in  the  plama  has  only  one  eonponiAt  idiieh  is  gitren  67 

1.,  (10.3» 


where  Z  o  is  the  free  space  inpedanoe. 


10.2  b  Fi*'»«1*  tha  ftir>~Ua  transfora  the  integral  edcpression  (10.17) 
to  OFlindrioal  ooordinates  in  both  tiie  transfora  and  oonfiguration  spaoes 

»tf 4  fig 

ofjS  sin0 

7  s  ^  sin  0  (10.34) 

s  s  9  oos  0 

and  obtain 


»xp 


2T 


oosfS  e 


t  k*  K.  *vrv*S 


000^  +  /n*  -  sin^ 


(10.35) 


where  n  is  the  appropriate  contow  in  the  (i  •plane. 

The  integrand  contains  the  radioal  (n^  •>  sin^  )^  as  a  resolt  of  whidh 
the  point  P»  9if  where 

0^  B  i:  arc  sin  n  (10.36) 

will  be  the  bransh  point.  Following  the  aethod  of  the  prerious  sections  we 
thsa  write 

»xo  *S  ♦  »4V  u  (®  -  »s  ) 


(10.37) 


162 


where  and  the  saddle^polnt  and  the  braneb-eut  eontrlbutlons 
respectively.  The  saddle>^polnt  Contribution  can  be  found  Unadlataily  ualng 
(9.12).  We  obtain 


**  V^2Tr  cos  6  +  /n*  «  sia*  0 


L  k«k  sirt*e '  ^ 

_ _ _  _ _ _  _  9  . . . 


/kTF 


To  evaluate  the  branch-cut  oontribution  we  use  the  fornulatlon  of  (9* 
wlte 

<>• 

where 


and 


How  we  find 


R  (*fi  )  »  (±fi  )  *  (*fS  ) 


(5  a  ©j  + 


(10.38) 


and 


(10.39) 


(10.40) 


R(±fi)  m  oosh,h/n*-siiiis  -isinkji>/n*-slnj8j 

^  IT  '  V  oosV-  (a*-ein*^  )  J 

Ve  set  |9  a  0,  everywhere  eeceept  In  the  radioal  wtaldh  beooaea  eero  at  that 
point.  We  oonpute  it  then  to  the  neact  higher  appapnartnatiun  and  obtain 


\/a^  «gin*^  'V'  ^  .  y/n  /  1  •  x 


(10.^) 


lAere  of  the  two  possible  choioes  of  the  sign  for  the  radioal  ve  have  chosen 
the  one  that  asyhes  In  |  (n^  «  ain^/8  )^  J  >  0.  Oubstltuting  the  last 
Sion  into  (10.41)  we  obtain  for  the  leadiog  tore 


yi??  ‘ 


X  • 


R  (*/l  )  -v 


(10.43) 
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Osliig  this  result  in  the  Integrand  of  (10.39)  gives  finslly 


\/W 


(1  -  Ik^h) 


g  L  rW.pcesO*l  -  Si)  1 

Sin  (101  -0,  )j^ 


(10.44) 


sM  note  that  this  lateral  field  is  of  second  order. 

The  magnetic  field  components  can  be  found  from  the  appropriate  curl 
equation.  The  on]^  non«vanishiBg  component  is  H  uhich  is  given  by 


H 


•d 


uhere  Zq  is  the  free^space  impedanee. 

The  Poyntlng  vector  has  a  single  component  given  by 


I  cos  e  e^*^**^*^*^*-*^*^*^  1^ 

iTTf  I  cos  ®  V  /n*  •  sin*  0  I 


(10.45) 


(10.46) 


10.3  CLOSURX 


In  this  chapter  ve  have  formulated  and  solved  the  problem  of  an  electric 
current  line  source  situated  in  a  magnetoplasma  vben  the  steady  magnetic  field 
is  along  the  i  jpe  source  and  the  line  sougrce  is  parallel  to  the  boundary*  In 
particular  ve  fowd  the  radiation  in  the  air  and  in  the  magnetoplasma  as  veil 
as  the  lateral  field  in  each  region  vhen  oonditione  are  favorable  for  its  ex'^ 
istsuice.  In  this  problea  the  steady  magne^e  field  has  no  effect  on  the 
plasma  since  the  electric  field  vector  is  alveys  in  the  sane  direction  as  the 
steady  aagnetio  field. 


PART  lU 

PULD  or  MIQNETIC  CURBEHT  LINS  SOURCES  IN  MAONSIOFLASMA 
Ulia  A  atPARATlOB  fiOOHOAlS 


G84PTIS  11 

BIGOBOUS  FOBMIUtIdll  OP  THE  PBCfflLiM  OP  k  MAGNETIC  CURliENT  LINE  SOURCE  WH^ 

THE  STEAQS  MAGNETIC  PIELO  IS  NORMAL  TO  IT 

Xn  thla  dhaptar  ve  ahaVl  be  ooneerbed  vith  the  finding  of  appropxlate 
integral  repreaentationa  for  the  Carteaian  eomponenta  of  the  field  vectors  for 
the  magnetoplaanw*  and  alr^half^apaoea.  This  %dll  entail  a  review  of  the 
fhndaiMntal  field  a^iatlona  and  the  definition  of  the  aouree  of  the  electro^ 

■agnetio  waves. 


11.1  STATEMENT  OP  THE  PROBLEM 


The  geometry  of  the  problem  ia  shown  in  Pigure  11.1.  The  horizontal 
plane  a  ^  o  ooineides  with  the  interface  between  the  anisotropic  homogeneous 
plaema  and  air.  For  oonvenlenoe  we  call  the  plasma  medium  (1)  and  the  air 
medium  (0).  As  beforSf  we  asauize  that  both  inedla  have  the  same  magnetic  in» 
dttotlve  eapacitgr  of  free  apace*  .  The  steady  magnetie  field  H  ^  will*  in 
this  case*  be  oriented  normal  to  the  line  source  whereas  the  plane  formed  by 
the  line  sovurce  and  the  direction  of  the  steady  magnetic  field  will  be  parai^ 
Hel  to  the  interface. 


11. E  FONGAISnTAL  EQUATIONS 

The  definition  of  the  present  boundary  valve  problem  impliea  the  solution 
to  Maxwell's  equations  subject  to  the  usual  boundary  conditions  at  the  inters 
face  and  proper  behavior  at  infinity^  The  use  of  the  auxiliary  vector 
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l^fur*  11,1  - 
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{iotentlals  doeis  not  seem  to  simplify  the  problem  aziy>  IhuSf  ve  shall  work 
%£Lth  the  field  components  directly. 

11.2  a  The  nature  of  the  source^For  the  purpose  of  this  problem  it  will 
be  assumed  that  the  sovuree  of  the  eleetromagnetio  wares  consists  of  a  very 
thin  straight  solenoid  of  infinite  extent  carrying  an  alternating  current 

It  is  wsH^loiown  (25»p.185)  that  suoh  a  solenoid  can  be  replaced  for 
donTenienoe  by  a  line  Oarrying  a  magoetie  ourrent  .  To  localize  the 

souroe  properly  we  shall  write  for  the  magnetio  oitfrent  density 

K  «r  (x)  ^  (s  +  h)  (11.1) 

wbere  £  (  )  is  the  Dirac  delta  funetion. 

11.2  b  T^a  fiai/i  equationa^^Aa  we  remarked  before^  our  source  of  elec» 
tromagnetio  Waves  in  this  irobleim  may  be  regarded  as  a  magnetic  current  line 
Singularity  at  a  point  x  »  0»  s  »  «  h.  Then  for  the  magpetoplasma  region,  the 
appcopriate  form  of  the  Maxwell's  equations  is 

^  X  1,  «  i  -  i^(x)^(a-fh)  Ty 

^  X  «  -  i  oo 

(11.2) 

*  0 

where  ^  is  the  fictitious  "magneto  charge  density^  related  to  the  magnetic 
ourrwat  density  by  the  continuity  equation 

(11*3) 

Sinoe  X  is  nen^singidary  the  second  equation  of  (11.2)  can  be  written 


^ vx  •  i  6,  F, 


(11.4) 
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Nov  peffondag  a  oorl  operation  on  the  last  equation  and  using  the  first 
equation  of  (It *2)1  ve  obtain  the  vatre  equation 


^  X  ^  X  H,  +  k*  H,  s  ^  S  (x)  (s  +  h)  V 


Using  the  inverse  of  the  pernittitrity  tensor  of  equation  (2*9)  and  carx^ng 
out  the  neoessaxy  slgebraio  operation  where  vi  note  that  0,  ve  obtain  a 
set  of  simultaneous  equattons  as  foUovst 


Xki+3»  -bKdkdi 

H*. 

0 

-CKhi,* 

Ny. 

1 

0 

(11.6) 


iMoh  ve  shall,  leave  in  this  fore  for  the  ties  being.  From  equation  (11. 4)  we 
note  that  the  electric  field  oomponente  oan  be  expressed  in  tezns  of  the  mag^ 
netio  field  components  as  foUowit 


0 

0 

”  H,..  * 

*y. 

-  1 

iK^ 

ak 

Hy, 

i<0€«% 

-iKdx 

i»<% 

H*.  ^ 

(11.7) 


the  air.-region  the  apF^priate  foms  of  the  Maxwell's  equations  is 


r  *  fo  ® 

V  J5  if,  »  *  i 
'ff,  BO 
lo  ^0 


(11.8) 


•nd  the  Cartesian  oomponaats  of  the  magnetio 
Tsetor  vafs  eqqaUoa  (  ^ 


.4)  j 


and  eleotrio  fields  satisfy^  the 

} 


»  Q  . 


(11,9) 
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11.3  FOURIER  INTEGRAL  REPRESlNZAlION  IN  CARTESIAN  COOROINATES 


The  formulation  of  the  present  boundary  ralue  problem  can  be  simplified 
a  great  deal  by  expressing  the  field  oomponsate  ia  the  magnetoplasnia  and  in 
the  edr  in  terms  of  their  double  Fourier  Integral  representation  in  Cartesian 
coordinates  in  the  transform  space  as  vail  as  in  the  configuration  space. 
ThttSf  as  beforst  ve  introduce  a  double  Fourier  transform  pair  defined  by 


F 


{ot,fe0 


f  dxds 


(in  10) 


and 


F  (x,t) 


ifl. 


(«.  *  *3)  •  dot.  daf3 


(mil) 


in  idiat  follow  ve  shall  also  need  the  transforms  of  the  derivatives. 
These  can  be  obtained  by  integrating  by  parts.  As  vas  show  in  Chapter  2  for 
the  triple  transfonif  the  vanishing  of  the  integrated  part  at  the  upper  and 
lover  limits  is  assured  providing  the  radiation  condition  is  satisfied.  Thus, 
ve  oan  establlsb  the  folloving  oorrespondencest 


«===^  i  ot.  ? 

9===»  iotj?. 


(11.12) 


11.3  a  The  Hot^ii»iBf>ndtTU7  to  thu  aniira^— Tn  tranS- 

fbn  the  inhomogeneous  ayhtm  of  the  simultaneous  differential  equations  in 
(11.6),  one  multiplies  both  sides  by  aaqp'£  1  (  w,z  i-hiss)}  and  integrates 
vith  respect  to  the  real  variablss  x  and  s  betvitn  «  ^  and  *  ,  The 

rlghtHhand  side  of  (11.6)  yields  at  onps 

S  (9  *h)  r^(^**^**>  m»  ■  e<«aH. 


(11.13) 


MO 


imd  the  leftwhaad  aide  trihafons  aeodrdihg  to  (11.12)*  Coaaequently,  oaa 
obtaina 
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«,Ka(,Q(3 
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H,. 


Hy, 
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(11*14) 


The  aboira  ogrhten  of  algebraic  equatioha  can  be  reduced  to  a  almpler  one  aa 


- 

fk*  1 

-  - 

0 

H,. 

0 

C«^-S‘Xai**sD 

0 

tf 

^  ki-«  *  *«!?•) 

Hy, 

2lr 

Sf., 

0 

(11*15) 


lAiera  a,  and  a^  are  given  bjr  (8.15)*  Nov,  ve  can  vrite  eadb  tranafortted  oo»> 
ponant  of  the  aagnetio  field  aa  foUovsi 

Jl(»  ^  ^  .*«,k  (11.16a) 

2T  I  (9^-af)(a.Vaf)  J 


B  ^  r  .1  (11.16b) 

2Tr  •  [(o<^-a.M(<x^-.a|)J 


S  ,  1  (11.160) 

2Tr  {^3^0?  )  ( of '•of )  J 


The  inversion  with  respect  to  the  ot,  transform  variable  can  be  perforBad 
lanediateiar.  According  to  (2.53)  and  (2.54) >  ve  obtain 


gdj  s  tuc.C  K  r  XkS«»otr*»8?  ,ts.l»*kl 
y,  2  /ZTT  I  8,  (a»  -a*  ) 


8»  <  0?  -8|  3r 


j,  (11,17a) 
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2  /W  ‘  la,(8»-8*) 


s<v  = » f#ii^ 

z4>HVd  2  ^/2Tr  8,*  “8^ 


»x(8.*  -"St  ) 


B,  —Bj 


(11.17b) 


(11.17c) 


Finally^  v8  Inyeirt  td.th  respect  to  the  «,  trahsfoim  variable  and  obtain  the 
desired  repreaentatlon 


s 

“y 


4Tr  Jt  8.  (a*  -a|  ) 

«B4d8 


qcoi.x 

«i(8.*  *®X  )  ■* 

(11.18a) 
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3*  f  8  . .  *  ^  V  .  le*^***  dor.  (11*18b) 

8.(8»-8‘)  8,(8.‘-8,M} 
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<’«.  I  b«.l 


doc,  • 


(ll.lSo) 


11.3  b  The  complTOwntary  field  in  the  daaMS'M^ln  the  preceding  sabseotion 
ve  obtain  the  "particular  integrals”  of  the  ^stem  of  equations  in  (11.15) 
vhich  represent  the  primary  eKoitation  due  to  the  source.  IC  satisfy 
boundazy  conditions  of  the  problemf  ue  shall  need  an  appropriate  ooatplementazy 
solution  of  corresponding  honogeneous  system. 

Pzcsi  (11.15)  it  is  clear  thatf[^>  satisfies 

(  )(  8»  )  »  Q  (11.19) 

vhsrs  in  the  above  equation  o(§is  loolesd  upon  as  a  differential  operator  with 


m 


mspeot  to  the  »•  oodrdinate.  tho  so^tlon  to  (11.19)  ^  be  vrltteft 
iBBediately  la  oomponent  form 


tjfc) 

“Kl 

Cm 

-a.* 

+  C.2 

ffc) 

Cai 

+  Cii 

e“‘»»* 

tfj 

C3. 

+  c,. 

(11.20) 


t&ere  ve  discarded  Solutions  %Q.th  positi^re  exponentials  Since  ve  can  have  only 
vaves  going  avey  from  the  interface  upon  reflection.  The  ooeffldients 
dan  be  found  by  a  method  identical  to  the  one  of  Section  2.3.  Ue  obtain 


i  s  oc.  3. 
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(11.21a} 


i  KOC.^ 


’9l 
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(11.21b) 


C.,  B 


Ik  PC,  a 


^  Xki  -  oc5  -  aj 


zt 


(11.21c) 
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•zz 


(11.21d) 


For  convenienoe  in  Wbat  foUoiiSf  ve  shall  nonBaliie  the  coefficients  C2, 
and  C22  88  foUovst 

0.,  •  *.  («•»•) 


*zz 


2 


(11.22b) 


Nov  ve  invert  vith  rsapect  to  the  oi(,  tienafom  yariable  and  obtain  the  deiAred 
representation 
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dot. 


(11.23a) 


H|,>  «  {{ -  :  t  hs:i-^ - I  e*^*.**  d«. 


(11.23b) 


g^)  »  ^  fr  a.A.e-*^**^ 

'  4lr  XkJ*etJ^8^  J 


dot..  (11. 23d) 


1 1 .3-6  The  field  J.a_^thd  air^or  the  alj*»regldn  we  choose  the  following 
Bolutida  to  the  wara  equation  (11.9)t 
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yo 


I 

4T 


Js,  .i 


dot, 


(in24a) 


«*•**-  J  e  dtx. 


(11.24b) 


doi, 


4T 
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(11.240) 


where  b^  has  the  saae  weaning  as  before.  One  can  readl^  show  that  the  shore 
field  oaqionants  also  satisfy  the  dlrergenoe  equation  Q. 


11 .4  1SI  BQONOUa  OQNQITIONS 


In  the  preoedlng  seotlons  we  fowd  field  oonponents  in  the  aagnetoplasina 
and  alr<-reglons  tddeh  are  solutions  to  Maawell's  eqiuatlons.  Moreorary  in  sol" 
Ting  HsKweH's  squatLens  we  dbpse  sutdi  solutions  for  the  field  representation 
that  hare  proper  bebaTlpr  at  Infinitqr  bj  requiring  that  the  laaglnary  part  of 
^e  pertinent  nponents  be  non-negatlTs.  la  addition,  the  field  ooaqipnents 
eontala  oertf^  %tm  far  undsterwined  ooeffisl«ats  which,  upon  ii^sitlon  of 
tbs  hoandtt^  oonditlpos,  will  be  dstendasd  and  thos  render  the  solution 
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unique. 

11  «4  a  Statement  of  the  boundaiy  coiadition^The  bouadaiy  conditions  to 
be  satisfiod  by  the  Cartesian  components  of  the  field  vectors  require  continue 
ity  of  the  tangential  components  of  ^e  electric  and  magnetic  fields  at  the 
interface  s  s  o.  In  terms  of  the  magnetic  field  oniy  these  can  be  vritten: 


Hx,  (11.25a) 

Hyo-Hy,  (11.25b) 

(11.25c) 

X(  .  (11.25d) 

The  above  equations  can  be  recast  in  a  more  eonvenlent  form  as  follovst 

0  *  ^  s^  Hy,  +  i  D,  tty,  (11.26a) 

0  a  +12)^  tt„,  ♦  iKO(,Hy,  +  ot.ili,  (11.26b) 


idwn  V  denotes  the  field  representations  with  the  integral  sign  removed. 

11. A  b  Applloation  Of  the  boundary  eonditiona~Performing  the  necessary 
algebraic  operations  to  satisfy  equations  (11.27a)  and  (11.26b)  gives  a  set  of 
two  simultaneous  equations  in  two  ucJcoowns  as  follows  t 


**■*"  $s* 

A. 

*s.  (s.  -  ^  s.)  R.  >v  ,  S.  J  R,  Cs*)  e 

s,sa(:s,»*5^; 

s.  +  s. 

S%  4?  im 

Sa(s.^s.Ye‘**>^  «  s.  s.Ye^«»*v 

PhC*0 

s.s*(rs.»,s^) 

J 

U  J 

m  m 

where 

(8,,i)  =  %k^  *  ot?  •  8,^ 


He  solve  the  above  aysten  of  equations  uslgg  Gramer's  ruls  and  obtain 


(11.28) 
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a.  (a.^  *8*  )  Hh 


vhtro 

a  (a,  +  ^ao)(ai+a^)i‘H(a,H8^+^fl|,)(d,'f0^)P^(8i).  (11.30) 


fb  flad  B,  and  um  (11.25a)  and  (I1.25b)  aa  vail  aa  the  abotra  raaulta 
and  obtain 


(11.31a) 


(11.31b) 


11.5  rUU)  OQMFOBBilS  XH  A39  III  CZLZHnaCAL  GCIOBDIlUlXiS 

Having  found  rlgoroua  axpraaalona  for  tha  flald  oonponenta  In  tha 
Pa^iian  ooordlnataa  In  boUi  raglona*  va  oan  now  tranafont  than  to  oyllpdrl.- 
oal  ooordlnataf  Ia  idddl  ttaaaa  flald  oonponanta  vUl  ha^  aoamhat  aii^ilar 
fbna. 

Ha  tnnafoni  to  agrlindrioal  opordlnataa  In  both  poi^ginatlon  and  trwia>» 
fona  Qaaaa  mdag 


. . 
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6t.  a  ko  ain/a 

s  a  ^  oos  Q  (11*32) 

X  n  p  sin  4 

As  a  ooDseq^oa  ve  obtain 

^  /‘(ia-»ao)g.  (a.  y6^^^s.-^a6)^(aa)^^^  ,  ck.fco«<p.e)^^ 

2T  J  Ih 

(11.33a) 

•*  2T  4^  >N 

(11 .33b) 

(11. 33c) 
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vhere  C,  is  the  integration  path  in  the  oonpleK  (3  'plane  as  Shoim  in  Fig.  3.1  * 
In  the  abore  equations  it  is  understood  that  s,  »  s^  »  and  s^  sure  all  func¬ 
tions  of /3  according  to  the  transforoation  (11.32). 


11.6  CLOSUIffi 

In  the  foregoing  chapter  ve  forBulated  rigorously  the  problem  of  a  mag- 
natic  current  line  source  in  a  Bagnetpplaaoa  vhen  the  steady  magnetic  field  is 
noznal  to  it.  In  particular,  ve  found  the  integral  representation  for  the 
field  components  in  the  Cartesian  and  cylindrical  coordinates  is  both  regions* 
Unlike  that  in  the  corresponding  isotropio  case,  aH  of  the  electric  and  mag' 
netie  field  components  are  present  is  both  regions. 


CHAPTER  12 

RESULTS  FOR  THE  AIR  REaiOH  'iHEti  'HiE  STEAOT  MAGNETIC  FIELD  IS 
HORtiAL  TO  THE  LBil  SOURCE 

Having  fox^ulated  the  pi^blem  In  the  previous  chapter^  we  are  now  ready 
to  evaluate  the  definite  integrals  to  obtain  explicit  foriDs  for  the  electric 
and  magnetic  field  components. 

Since  the  present  problem  is  very  similar  to  the  corresponding  problem  of 
an  electric  line  source  of  Chapters  9  and  10,  we  shall  use  many  of  those  r»> 
salts  that  are  applicable  to  the  present  problem. 

12.1  E7ALUAI10N  OF  THE  FIELD  INTEGRALS 


The  field  integrals  to  be  evaluated  have  the  following  form: 


ojs.^  K  r 
2T  J 

Gy^coSj8e^*‘»fc«»»(p-«)  dS 

(12.1a) 

<oe«^  K  I  .. 
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(12.1c) 
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(12.2a) 
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tiui 

(or  +4r60S|ft  )(es  +  cos/S)P„(er,  )*(esrj  +  S‘G68|S)(a7  +  eos/a )P„{fl|)  (12.3a) 


K  ( o-.,x  )  -  sin’^/3  * 


(12.3b) 


and  tf,  and  are  given  by  (9.2a)*  We  also  note  that  P^(<^)  and  P^(Oi)  can 
be  uiltten  In  the  fora 


(12^4) 


12.1  a  Braluatloa  of  the  field  lntegrala'.at„  the  saddle  point^Apalving 
the  formuLLatlon  (9*12)  to  the  field  integrals  in  (12.1)  ve  obtain 
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(12.5a) 
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H<^)  s  0.„(G)OOS  G  2l^j^±^iaa 


(12.5b) 


(12.5c) 


12«1  b  Evaluation  of  the  field  integrals  along  the  branch  cuta^^Slnce  e» 
and  ax^  the  same  as  in  the  corresponding  problem  of  the  eieotric  line 
source^  the  branch  points  are  the  same  as  in  Chapter  9  and  the  fomulatipn  of 
(9.21)  applies.  Applying  this  foimulatipn  to  the  field  integzals  of  (12.1 )« 
ve  obtain  for  the  components  of  the  lateral  field 


Vrr  '  sinOGI  .  G,.  )1^ 
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yo  72^  [k,f  sin  (IGI  « ^ ) 3 ^ 
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g n  / 1 -C6-7;  1  (la.To) 

/&•»•€ b*i  I  J 
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*  \/i» C6 it7){[»* ?0*6^7)]y^ (C^ 6 ±7)  - r[i-?-C6±7)]Vi* (r4“*7)  }  (12*7d) 
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/•)  *  eoa  y/t  3in20,. 
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.  2  f  /-k(6+»j+-? )  +  C  /l^(a  +  7)  ]  (12.9a) 
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12.2  fovm.  FLOW 


In  the  radiation  field  all  components  of  the  electric  and  magnetic  field 
are  transtrerse  to  the  radial  direction.  Therefore^  the  Poynting  vector  baa 
only  a  radial  component  gltren  by 

Sf  *-T-{  ivl*  ♦  l»«l*  1  • 

Nov,  using  the  results  of  (12.5)  the  above  result  oaa  be  written  more 
explicitly  as  foUovs: 

s^,  - j  |o  (e)l'  +  kWb  |o,,(e)|*|  (w.ii) 

4  no  c  '  J 

\diere  Gy^  (0)  and  G*p  («)  are  given  by  (12.2). 

12.3  CLOSURE 

In  this  chapter  ve  found  the  field  components  in  the  alr-reglon  due  to  a 
magnetic  current  line  source  in  a  magnetoplaaDBa  idien  the  steady  nagnetio  field 
is  normal  to  the  line  source.  As  in  the  corresponding  case  of  an  electric 

current  line  source  in  Chapter  9,  ^9  field  in  the  air^reglon  oonaists  of  two 
main  oomponents,  the<  radiation  field  and  the  lateral  field. 
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the  radiation  field  contains,  in  addition  to  the  components  present  In 
the  oorrespehdlng  Isotropic  case,  additional  transverse  Components  of  the 
electric  and  magnetic  fields  vdiich  render  the  resulting  fields  to  be  elliptl- 
oally  polarlied  evexyidiere  emoept  right  above  the  seitfce« 


CHAPTER 


THE  STIAOT  }iAGN£TiC  PZm)  PARALLEL  TO  THE  LINE  SO  jRCE 

In  this  eha|>tef  we  shaLl  be  eoncexiied  id.tb  finding  the  electilc  and  mag^ 
netie  fields  due  to  a  magnetic  line  source  situated  in  a  magnetopilasma  when 
the  steadjr  magnetic  field  is  along  t^e  direction  of  the  line  source^ 

13.1  FOB»a^TlON  OF  THE  PBOBLM 

The  geometzy  of  the  problen  la  shown  in  Fig.  13.1*  The  horizontal  plana 
z  s  0  colnoides  with  the  Interface  between  the  anisotropic  homogeneous  plasma 
and  air.  for  Convenience  we  shall  call  the  plasma  medium  (1)  and  the  air 
medium  (O).  As  beforoi  we  assume  that. both  media  have  the  same  magnetic  in¬ 
ductive  capacity  of  free  space,  •  The  steady  magnetic  field  will,  in  this 
case,  be  oriented  along  the  line  source  idiereas  the  line  source  is  parallel  to 
the  interface. 


13i.l  fr 


solution  to  the  present  boundary  value 


prohliD  entails  solution  to  the  Hsmuell^s  equations  in  the  foUowing  foxmt 


^  i  ^  (y)  5  (i  ♦  h) 

^xIb,  * -icje«Tli 
«  0 


(13.1) 


where  Et  as  before,  denotea  the  mugnetio  ourr^at  end  ,  the  nagnetie  oharge. 
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llgnrf  13.1  «  OuaamiiBj  of  the  Problai  of  4  Hegnotlo  Current  I4ne  Sooroe  Uhen 
The  9te^  Megpetle  if  Fexellol 


In  the  alz^reglon  the  fields  must  satisfy 

^  X  =  i  CaJ^A^'S’^ 

^  X  H  p  =  *  i 
^•D*  =  0 
^•1*  =0 


(13.2) 


Moredver,  oiur  solution  nust  Inswe  pzx>per  behaiH.or  of  the  fields  at  infinity 
08  well  as  the  oontlnuity  of  the  tangential  oomponents  of  the  fields  at  the 
interface. 


13.1b  the  integral  representation  of  the  field  Components  in  the 
BLLft<OTa~Ua  form  the  wa/e  equation  by  performing  a  curl  operation  on  the  second 
equation  of  (13.1)  and  then  substitating  the  first.  This  giires 

*  X  r-‘  ^  X  f,  +  4?,  =  -  i  cut*.K  S(y)  ^  (s  +  h)  (13.3) 

Noting  that  0  and  using  the  components  Of  the  in/erse  of  the  permit* 
tivlty  tensor  from  equation  (2.9) »  we  obtain 


XkS+dy+h^  0 
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Hy, 

=  ^(y)  s  (»*h) 
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0  ^  3x 

Ck^+3/ 

0 

(13.4) 


From  (13.4)  it  appears  that  a  single  oon^nent  of  the  magnetic  field  may  be 
sufficient  to  solve  the  present  boondsrjr  veloe  prpblsii.  The  oorrespondlng 
differential  equation  is 


(xk|  +  ^  i&>«,x*^(y)^(i^) 


(13.5) 


FOr  opovenienoe  we  Jyntrpdupe  the  Fourier  transform  pair 
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a;, .  H.,  .--.X  ^ 


(13.6) 


idiiereupon  \te  tr^sform  eqiiatldn  (13*  5)  as  follovs: 


(  +  sf 


i  cac  o  X 

\/W 


{\h1) 


Si^ace  the  Irani ahlng  of  the  Integrated  pairt  is  assured  vhen  the  radiation  con¬ 
dition  is  satisfied.  The  parameter  s,  is  gi\ren  hy 


S,  =  \/  XkJ  *a|^' 
in  {  a,}  >  0  . 


(13.8) 


The  aolutlon  to  the  differential  equation  (13.7)  iS  velX-knovn  and  ve  urlte 


it  at  onoe 


h‘7  =- 


co&iXK  e 


tSiU-f-  K.I 


2  \/2T  s, 


(13.9) 


Inverting  vith  respect  to  the  transfom  variable,  ve  obtain  the  particular 


uW  - 
“ki  “ 


L  (axy  S,|«-»vl  ) 


(13.10) 


The  font  of  the  popplementaxy  field  in  the  plasna  iaiaediately  suggests  itself. 


He  vrite 


gCe)  _  _  ^ 


^A,  e  dajt 


(13.11) 


He  also  note  that  the  eleetxio  field 
^  s  0.  He  get 


om  be  found  fron  (2.27)  vith 
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£  S  0 


-•f  ®/i  r» 


3^«lKdt  0 


(13.12) 


In  th#  air  the  mag^tie  field  e<»ponents  aatlsf)r  the  homogeneous  wave  equation 


{  +  )  H^*0. 


IliUiy  ve  can  vrlte 


4^ 

f‘-‘ 


(13.13) 


(13.14) 


13.1-  d  The  boundary  Qonditiona»»th6  raqulrenent  of  the  continuity  of  the 
tangential  field  eomponents  at  the  boundary  a  s  0  amounts  to  satisfying  the 


following  equations} 


*  (iKa^  +  &«) 


(13.15) 


Upon  performing  the  indicated  differentiationy  «e  obtain  tvo  simultaneous 


equations  in  two  unknowns 


S,-i»^i 


s,  +  iKai 


(13.16) 


fyoiB  vliiab  ws  obtain 


Ao  » 


I  CS« 


S,  ♦  (Xa, •  iKof^) 


(13.17a) 


Z  a. 

s,  +  (Xa«  iK  ) 


(13.17b) 


Substituting  abpTa  raptilta  in  tba  ponaaponding  l^tegralsy  me  obtain  ^e 
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Hx,  =- 


cuegicK  (  u(i) 


I  (  /STlCo  ?,  )  *  (  v/Tk,  ft  )  ]  +  (13.21a) 


Where 


(13.21b) 


gj  4-  cocoyK  r  003 e» e d 3  (13.2 

**'  2T  J  GGs/i  +  X^^^"3in*/s  -  iK-sinfi 

the  Integrand  in  (13.21b)  contains  the  radical  (•^  «  ain^/&  )'^  as  a  re« 


suit  of  which  the  point  /%  -  ^ere 
9^  -  Ss  arc  sin  (‘"w) 


(13.22) 


Wdll  be  a  branch  point.  At  each  point  (A  the  integrand  in  (13.21b)  can  take 
on  two  values  depending  on  which  Sign  we  choose  for  the  radical.  As  before, 
it  is  convenient  here  to  talk  about  two  sheets  of  the  fi  -plane  (fonned  by  a 
two-sheeted  Rienann  surface)  on  which  the  integrand  is  Slnfpie  valoed.  The  two 


sheets  will  be  joined  along  the  lines 


Im  ^  -  sin*|§  ^  «*  0 


(13.23) 


starting  at  the  branch  point. 

As  before,  the  result  will  consist  of  the  sculdle-point  contribution  axui 
the  contribution  from  the  branch  cut  representing  the  radiation  field  and  the 
lateral  field  respective!^'.  Thus,  we  write 


HxT^  -Ki7^+Hl7>u  (  9-9e) 

and  evaluate  each  one  of  these  contributions  separately. 

To  find  the  saddle-point  contribution  we  use  (9*12)  and  obtain 


(13.24) 


\/zT  cos  9  +  -iKsin  9 


0  C 


(13.25) 


The  braneb-out  contribution  wae  fonNlated  in  (9.21a) .  Using  this  fonulation 
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ve  can  \alte 


(13.26) 


^ere 


R  (ai/i  )  =  (a:^  )  -  F^.  (^/i  ) 


(13.27) 


(3  =  flg  ♦ 


Nov  ve  find 


R  (i(b  )  = 


^  J£0  < 


003  fl  \/4.  *•  aln^g  ^ 

(eos^  ?  Ik  3ln(8  )*  -  X*(i|^  *  sin*/S  ) 


(13.23) 


Ve  oan  set  ^  Og  eveiyVbere  eKoept  in  the  radical  which  becoiaes  zero  at 
that  point.  We  ooapute  It  then  to  the  next  higher  approximation.  We  obtain 


-  3in^/5  ••  ^  ^ 


(13.29) 


where  of  the  two  possible  signs  of  the  radical  ve  have  chosen  the  one  that 


ookes  lm{  /jj-  -  8in*/5  }  >  0.  We  find  for  the  leading  term 

TT  (  /  TC  ^  1  +  !••<)*■ 


(13.30) 


Std)stitating  the  above  result  into  (3.26)  and  perfoxning  the  integration  t^iere 
ve  use  the  result  of  (4.31b) »  we  obtain  finally 


g<rsj  _  -<oeoX*K(7c^l^exp^  i  [  v5rk,ft.  oosdOl  ^  Wt)  *  Va 

*'  ( v?c^^iK)^fVx  ko  f*  ^  (i»i  -  »s  yj^ 


(13.31) 


We  observe  that  the  lateral  wave  in  this  case  is  not  syxpetric  as  before  but 
it  is  still  of  second  order. 

fo  s^snaariae  the  above  resulta  we  write  for  the  first  order  magnetic 


field  in  the  plasma 
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vlU  now  be  the  branoh  point*  Pollouing  the  method  of  the  preTioue  seetioii8> 
ve  write 

Hx6=  +  hW  tt  (  0  -  08  )  (13.37) 

Where  and  ere  the  eaddle«point  and  braneb*cut  eontrlbutiona  re<^ 
apeotlvelj.  The  eaddleopoint  oontribution  can  be  found  innedlately  using 
(9.12).  we  obtain 


uCs) 

“xo 


_  . 

^C^S?0  XoOS  0  a  i< Sin  0 


(13.38) 


To  evaluate  the  branoh.-eut  oontribution  we  use  the  formulation  of  equation 
(9.21a)  and  write 

(13.39) 


Hrt)  a  ie^‘‘.^.ooa^ei..e.)  j 


where 


and 


R  (±/5)  »  f*(±/a)  *F^(M) 

(5s  0 

2  • 


(13.40) 


Nov  we  find 

-«js,yZooe6[\/x-  sirfil'oo8(k,h  0p..eijnF<S)..i(XOQap*iicsinfi)3in(lc,h  v^-.Biij)9)] 


T  [(x  ooait  7  iK aia/a )  ..(x 


(13.41) 


We  set  A  ^  0ii  eveiTvhere  emoept  in  the  radical  idiich  becomes  aerp  at  that 
point.  We  oomputa  it  then  to  the  next  falser  approxiinatien  and  obtain 

/X-oin*(S  ^  e'"^^  Vx(1  -  X)  x  (13.42) 

where  of  the  two  poaaible  oboicee  of  the  sign  for  the  radioal  we  have  chosen 
^e  one  that  maloss  Im{  ain-fi  )'^  j  >  Q.  Suba'U.tuting  the  above  expreaidoa 
into  (13*41 )f  we  obtain  for  the  leading  tent 
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gtep  ) ^  [  1-1  >gM( 


(13.43) 


Sabstitutlag  the  last  result  into  the  Integral  of  (13.39)  end  perfoxioing  the 
Integration  uhere  we  use  the  results  of  (4.3lb}y  we  obtain  flnemy 


[<•)  ss  cmSoK  yx(1*x)^  ri*i\i^k«h(  v^%(1»>c)$  IK 

eio”  \/5lr  [  /X(1-x)  ^1  <  r  D^o  ?  sin(  l«l  -  ig)] ‘  * 


We  note  that  the  lateral  wave  in  the  air,  when  It  exists^  is  not  syimnetrio  but 
It  is  still  of  second  order. 

The  eleotrlo  field  conponents  in  the  air  can  be  found  from  the  appropriate 
curl  equation.  The  only  non^Vanishing  eoaponent  is  £«o  given  by 

Eeo*  *  K  flxo  (13»45) 

to  the  first  order. 

The  Panting  vector  has  only  a  radical  component  and  is  given  by 

^  4T^  I 


cos  »  e"***^^  ^ 

v'X*3in*fii  +  Xoos  0  ^  Ik  sin  0 


(13.46) 


Examination  of  (13.46)  reveals  some  interesting  possibilities.  For  when 
y,  >  sln^O  one  can  vrlte  for  a  lossless  plasma. 


S,  . 


cu€oX^K* 


and  lAen  %<  sii^O 


_ 

(  />c-8in*9  +  X  BOB  9)*  +  K^sln^  9 


(13.47a) 


o  _  we.x’E^ 
- - 


00^®  0  -  ^  koK  y/sCrx^0  —  % 


X^cos*®  +  (  /sin*'®  »-  x  ksln*0) 


,»Ql  t 


(13.47b) 


Then  It  follows  that  for  X>  sin^®  the  power  pattern  will  be  syn>ostrie 
about  the  axis  ®  s  o  and  unaffected  by  the  depth  of  the  spuroe  burial,  h. 
Vben  X<  sin^®  the  opposite  is  true;  the  pattern,  in  general,  will  not  be 
syneetrie  and  its  magnitude  will  deorease  with  the  Inoreased  depth  of  the 
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W.imr«  13.4  -  Pbwer  Pattern  in  Air  of  a  Magnetic  Cuurrent  Une  Source  in  Mag- 
Mtulum  H«-aXi»g  th«  J4»  So«ro».  II  =  750  P«  ®»M,«  ow^tw, 

hVSTw  0,  .*  gmM.  t .  10^  5  »  10»  «5rae.  pw  wi>4 


196 


198 


source  burial,  h.  Moreover,  ue  note  that  the  reversal  of  the  direction  of  the 
Stea^  magnetic  field  which  chasoges  the  sign  of  "k”  has  the  same  effect  as 
the  reversing  of  the  sign  of  the  angle  G.  Thus,  the  pattern  corresponding  to 
has  a  minor  symnetiy  with  respect  to  the  one  oonespending  to  Hp^  . 
lb  lilustrate  the  above  conelosions  a  nufflexd.eal  example  was  considered* 

It  ws  assumed  that  the  soturoe  Is  situated  In  the  lower  edge  of  the  Ionosphere 
which,  for  the  purpose  of  this  example,  was  assumed  to  be  homogeneous  and 
sharply  bounded,  having  an  electron  density  N  ~  750  electrons  per  cubic 
oentimeter  and  the  earth's  magaetic  field  *4  gauss.  The  povier  patterns 

for  the  depth  of  burial  of  100  meters  and  sero  meters  for  several  frequencies 
are  showa  In  Figures  13*3  *  13*8  idiere,  for  comparison  also,  the  power  pat^ 
terns  corresponding  to  0  were  superimposed. 

The  physical  Interpretation  of  these  patterns  Is  as  follows.  For  polar 
angles  smaller  than  the  critical  euigle,  the  transmitted  rays  undergo  the  usual 
refraction  fixed  by  Snell's  lav  for  real  angles 


^sin  G,ss  sin  0  (13*48} 

idiere  0,  Is  the  Incidence  angle  in  the  magnetoplasma.  Since  the  radiation 
from  the  line  source  Itself  Is  symmetric,  the  resultant  pattern  Is  also  sya^ 
metric  within  the  wedge  region  Gp<  G  <  0^  .  Beyond  the  critical  angle  the 
transmitted  energy  is  produced  by  the  Ix^mogeneous  waves  which  are  also  being 
ananated  from  the  souroe.  A  pombined  aoUon  of  tbe  interface  and  the  steady 
magnetic  field  (see  equations  (13.15)  and  (13*17)  destroys  the  symmetry  of  the 
inhomogeneous  wave  isplitudes  In  the  magnetoplasma  ^ch  then  results  In  lack 
of  symmetiy  of  the  power  patterns  In  the  air  In  the  regions  beyond  the  critl<^ 
eal  angle*  An  excellent  oorrelation  between  the  values  of  the  orltlcal  angle 
in  Fig*  13*2  cuid  the  power  patterns  in  Figures  13*3  ^  13*5  can  be  noted. 

The  oritical  uigls  Gg  does  opt  hnrp  a  ri«l  ui^^  frequency  reaches 
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about  52  z  10^  ogrbles  par  seoond.  Thusi  the  patterns  In  Fig.  13*3  and  cor¬ 
responding  to  frequeneles  10^  and  5  z  10^  ^cles  per  seoond  are  not  synmetrlo 
alnoe  they  are  solely  produced  by  the  lnhoiK>geneous  wa7es  emanating  from  the 
souree.  At  a  frequency  of  10^  cycles  per  second^  the  critical  an^e  9^-61 
degrees  and  It  Is  apparent  from  the  upper  pattern  of  Fig*  13*4  that  this  Is 
the  point  idiere  the  non^iymmetiy  of  the  pattern  begins. 

Finally,  ue  mention  about  the  effect  of  the  depth  of  the  soiree  burial* 
It  is  appenrent,  especially  at  lev  frequencies  that  the  transmission  through 
the  ionosphere  Is  Improred  considerably  due  to  the  action  of  the  steady  mag- 
netlo  field.  Whereas  the  magnlt'ide  of  the  transmitted  power  diminishes  with 
the  Increased  depth  of  the  source's  btulal  If  no  steady  magnetic  field  Is 
present,  it  is  almost  independent  of  it  when  the  steady  magnetic  field  Is 
present*  This  Is  due  to  the  changes  produced  by  the  steady  magnetic  field  In 
the  effective  dielectric  constant  of  the  plasma. 


13*3  CLosass 


In  this  chapter  we  have  formulated  and  solved  the  problon  of  a  magnetic 
oiureat  line  source  situated  In  a  magpetoplasma  vben  the  steady  magnetic  field 
is  along  the  line  soiu'ce  and  the  line  source  Is  parallel  to  the  boundary.  In 
partleulari  we  found  the  radiation  field  In  the  air  and  in  the  magnetoplasma 
as  well  as  the  latend  field  in  each  region  when  the  oonditlons  are  favoraUe 
for  its  eKLstence* 

The  resultant  radiation  pattern  In  air  are  not  symmetric  at  lower  fren- 
quenoies  but  they  become  GQnnDetrio  with  increasing  frequency.  Moreover,  we 
found  the  presence  of  the  steady  magnetic  fidld,  eren  of  small  magnitude, 
lapcove  the  transmission  tlnough  the  plasma  a  great  deal  in  oonparison  with 
the  ease  of  no  steady  magnetio  field* 


CHAPTER  14 
GONCLCSIORS 


In  the  ^ceding  chapters  ue  hare  ahom  the  results  of  sueeessful 
analyses  of  sereral  eases  of  radiation  from  aburoes  In  the  presence  of  an 
anlsotrople  plaaEoa  half-space. 

In  Peurt  1  uhleh  includes  Chapters  2  -  7>  we  were  ooneemed  with  the  pro- 
blflB  of  a  horizontal  magnetic  dipole  in  and  out  of  mag^toplasma  \dien  the 
steady  magnetic  field  is  parallel  to  the  axis  of  the  dipole.  The  rigorous 
fomulation  of  the  problen  was  carried  out  to  the  point  tdiere  the  detennina- 
tion  of  the  pertinent  boundary  coefficients  remained  to  be  a  straight-forward* 
but  not  simple*  algebralo  process.  Due  to  the  inohibitlre  algebraic  oompleK- 
i^  involved  in  the  expUclt  finding  of  tne  necessary  boundary  aoeffleients* 
the  high  frequency  approximation  was  introduced  and  the  boundary  ooeffloients 
pertinent  to  the  air-region  subsequently  found.  In  this  prooess  it  was  found 
possible  to  separate  the  oontributions  of  the  plaama's  aaLsotzopy  soEplioitly* 
The  radiation  field  in  the  air^region  was  found  wLtbia  the  restrictions  of  the 
high  frequency  approximation  but  the  lateral  field  could  not  be  evaluated  in 
this  fomulation  due  to  the  fact  that  the  main  contribution  to  the  lateral 
field  cosies  from  the  point  idiere  the  high  frequency  approximation  is  not 
valid.  The  result  of  this  analysis  was  than  applied  to  the  proUan  of  a 
horizontal  magnetic  dipole  situated  near  the  lower  edge  in  the  Ionosphere.  It 
appeared  that  in  the  region  of  wave  frequencies  at  iddoh  the  high  frequency 
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appradaatlon  is  valid,  the  eorreotlon  due  to  the  earth's  magnetic  field 
aaiilfests  itself  moat  strongly  in  the  polar  plane  through  the  aids  of  the 
dipole  and  lr/4  radians  from  it.  Furthermore,  this  oorreetion  appears  to  be 
signlfloant  only  for  large  polar  angles,  i,e.,  in  the  regions  close  to  the 
interface. 

Part  II  idiloh  oonsists  of  Chapters  3  10  is  devoted  to  the  problems  of 

eleetrio  current  line  sovuroes  in  the  magnetoplasma  half-space.  The  case  of 
the  steady  magnetic  field  normal  to  the  line  soviree  is  rigorously  fozmulated 
In  Chapter  3  and  the  field  integrals  pertaining  to  the  airBregion  are  evalu¬ 
ated  in  Chapter  9  vlth  no  restriction  on  the  vave  freqiiendy.  ^e  radiation 
field  contains,  in  addition  to  the  components  present  in  the  corresponding 
Isotropic  case,  additional  transverse  components  of  the  electric  and  magnetic 
field  ehieh  render  the  resulting  fields  to  be  eliptioally  polarised  everyidiere 
«eoept  right  above  the  source,  n^e  lateral  field  is  of  second  order  as  in  the 
corresponding  isotropio  case,  however,  unlike  in  that  ease,  it  may  consist  of 
two  waves  rather  than  one.  These  waves  propagate  along  the  interface  with 
phase  velocities  characteristic  of  the  rights  and  left-hand  circularly  polar¬ 
ised  plane  waves  in  the  magnetoplaaaa.  From  the  numerical  eatainple  considered 
it  is  ooneluded  that  the  transmission  of  energy  txon  the  electric  current  line 
source  throu|di  the  pi, asms  and  into  the  free  space  is  generally  enhanced  by  the 
presence  of  the  steady  magnetic  field  especially  at  wave  frequencies  below 
plasma  frequency.  At  wave  frequencies  well  above  the  plasma  frequency  the 
aetLon  of  the  steady  imgnetio  field  becomes  insignificant  which  is  not  at  all 
surprising.  Ibrcover,  the  radiation  patterns  omhibit  peaks  at  points  cor- 
resppnding  to  the  (oitical  angles  which  in  turn  correspond  to  the  indloes  of 
rsfraotlon  of  the  ri^t-  end  loft-hand  oiroularly  polariaod  waves  in  the  mag- 
natopLasma.  Tha  remaining  obapier  of  Fart  II,  Gbaptar  10,  is  devoted  to  the 
fbrwaatioa  and  solntion  of  tte  prohlsm  of  an  alsotrie  eorrent  line  souroe  in 


the  magnetoplasna  tdiea  the  steady  nagnetie  field  is  peurallel  to  the  line 
Source.  In  this  case  it  is  found  that  the  steady  magnetic  field  has  no  effect 
on  the  radiation  field.  This  is  consistent  idth  \duit  one  would  expect  since 
in  this  Oase  the  electric  field  is  always  along  the  steady  magnetic  fields 
thuS}  the  cross-product  H  x  vanishes  and  no  interaction  with  the  steady 
aagnetic  field  takes  place. 

Part  Ill  consisting  of  Chapters  11  «  13  is  dstroted  to  the  problems  of 
magnetic  current  line  sources,  the  case  of  the  steady  magsetie  field  normal 
to  the  line  source  is  formulated  rigorously  in  Chapter  11  and  the  integrals 
are  etialuated  in  Chapter  12  with  no  restriction  on  frequen^.  This  ease  is 
analogous  to  the  corresponding  case  of  Part  11  and  the  conclusions  drawn  in 
that  psurb  apply  here  also.  Chapter  13  is  devoted  to  the  formulation  and 
solution  of  the  ease  when  the  steady  magnetic  field  is  along  the  line  soturee. 
Here  it  is  found  that  the  magnetoplasma  acts  as  a  single-refracting  medium 
i^oh  is  consistent  with  physical  reasoning  since  in  this  case  the  steadjr 
magnetic  field  is  always  normal  to  the  alternating  electric  field)  both  field 
vectors  being  in  a  plane  normal  to  the  direction  of  propagation.  The  action 
of  the  interface  produces)  however)  a  very  interesting  phenomenon  in  the  form 
of  the  transmission  coefficient  idiich  is  not  symoetrlo  with  respect  to  the 
vertical  axis  thxou^  the  source.  Moreover)  the  magnitude  of  the  tranamlssipn 
opeffLeient  and  the  phase  nay)  ooder  pertain  conditions)  depend  on  the  orien¬ 
tation  of  the  steedjr  magnetle  field.  This  produces  some  rather  interesting 
phenqiaena  in  the  radiatipn  pattern  in  the  alr-regipn.  When  the  effective 
indax  of  refraction  in  the  msgnetoplaama  is  Imagiiuuy  which  occurs  at  very  low 
frequencies)  the  resulting  radiation  pattern  is  not  synobtrie  at  all.  More¬ 
over)  in  this  case  the  steady  magnetic  field  reverses  its  orientation) 
the  radiation  pattern  tehes  on  a  form  that  is  esaetly  the  mirror  image  about 
the  vertical  aadla  of  Vsa  radiation  pattern  before  the  steady  magnetio  field 
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Has  reimrsed.  When  the  effective  Ixidex  of  refraction  In  the  raagnetopleana  le 
between  sera  and  one,  the  radiation  pattern  is  symmetric  in  the  wedge  region 
about  the  vertical  axis  bounded  by  the  lines  representing  the  critical  angles 
but  it  is  not  s^omatrlc  between  these  lines  and  the  interface.  Finally,  vdien 
the  effective  index  of  refraction  in  magnetoplasma  is  greater  than  one,  there 
is  no  critical  an^  and  the  pattern  is  perfectly  symmetric#  The  physical 
Intes^sretation  of  these  phenomena  is  as  follows.  For  polar  angles  smaller 
than  the  critical  angle,  the  transmitted  rays  undergo  the  usual  refraction 
fixed  by  dnell’e  law  for  real  an^Lea.  Now  ainee  the  radiation  from  the  source 
itself  le  SyBnetrio,  the  resultant  portion  of  the  pattera  in  the  air  is  also 
synastrio  within  the  wedge  region  limited  by  the  lines  representing  the 
erlticsl  angles.  For  the  trananisaion  angles  beyond  the  critical  angles,  the 
Snell's  law  can  be  satisfied  only  by  oomplex  angles  in  the  magnetoplasma  idilch 
that  the  transmitted  ene:^  in  the  air  is  produced  by  the  inhomogeneous 
waves  in  the  magnetoplasma  which  are  also  being  emanated  from  the  source.  A 
oomhlned  action  of  the  interface  and  the  steady  magnetic  field  destroys  the 
symmetry  of  the  Inhamogeneoua  wave  amplitudes  in  the  magnetoplaama  idilch  then 
results  in  lack  of  symmetry  of  the  radiation  pattern  in  the  air  in  the  regions 
beyond  the  orltlcal  angles. 

Altitough  the  prphlems  treated  in  Part  II  and  111  are  two»»dimeasional, 
they  give  a  great  deal  of  insight  into  the  corresponding  three>i^ynensipnal 
problems*  lbs  fpllpwing  oorrospondsnoes  can  be  easily  deduced  by  physical 
rsasoning; 

(a)  The  radiation  pattern  from  an  electric  current  line  aouroe  with 
Mnml  to  the  line  sowoe  is  the  same  as  it  would  result  from  a 
eleotrle  dipole  in  the  jdane  normal  to  the  dipole  axis 
when  Hp^  is  normal  to  the  dipols  eeds* 


Nov  vs  eonslder 


J  (a  + 


g 


idlers  r  is  the  path  in  the  com^^  a^plane  as  shown  in  Fig.  A.1.  Sinee  the 
Oontour  oontains  no  singiilazltieSf  then  ty  Cauehy' s 


Figure  A.1 


theojrem  (15*  p*  30) 
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Furtb.morei  on  r:  v.  put  #  s  ^  1  +  <  ,  '^2  > «  >  -  ^2  wd  Qbtain 
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n,  % 


On  Cj  we  put  z  =  1  +  ee"'*-®  ,  -  ^/2  >  0  >  -  3  ^2  and  obtain 

r . if d 0  ^ ml  . 

J  J  (2+  €e-‘»)  ee-®  - - .^0  2 

Thus,  by  (A. 10),  (A.11),  and  (A. 12)  we  obtain 
I  s  TTi  [  1  H)*"  ]  . 

How  substituting  (A.  13)  ixito  (A.  3)  gives 


(A.11) 


(A.12) 


(A. 13) 


I*  TTi^ 

?or  integer  values  of  n  the  following  relation  is  Valid  (22,  p.  357) 

J.,^-(-irJ^.  (A.15) 

Thus,  we  obtain  finally 

I  a  4Tri2~  (-ir  sin(2n+  1  )<f  (Af )  .  (A.16) 

n«0 


[i  -  (-ir]  j^(Af) 


(A.  14) 
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